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Abstract 


The  role  of  composite  materials  in  modem  hclicoptei  biade  design  has  become  most 
important  during  recent  years.  By  exploiting  the  direction.  ]  stiffness  properties  of  these 
composites,  favorable  torsion  modes  or  "twisting"  can  be  achieved.  The  capability  to 
apply  this  potential  as  a  design  parameter  is  generally  known  as  aeroelastic  tailoring. 

The  bending" torsional  coupling  of  static,  hingeks;^  c  omposite  rotor  blades  is  investi¬ 
gated  using  finite  element  theory.  The  hingeless  bladi  i-,  treated  as  a  single  cell  lam¬ 
inated  shell  beam.  Each  laminate  is  composed  of  duT.rent  lay-ups  of  graphite-epoxy 
composite  plies  and  is  categorized  as  isotropic  or  anisot.xpic  based  upon  this  lay-up. 


A  systematic  study  is  made  tc  identify  the  effects  of  ply  orientation  and  lamina 
thickness  on  blade  section  properties.  The  results  of  this  study  are  used  to  solve  the  beam 
equations  for  composite  materials.  /First,  the  beam  is  modeled  and  the  fiber  orientation 
and  thickness  variations  are  selected.  This  is  done  using  PATRAN®  by  PDA  Engineer¬ 
ing  as  both  a  solid  modeler  and  as  a  pre-pi’ocessor.  The  model  then  undergoes  finite  ele¬ 
ment  analysis  and  post-processing  using  ANSYS®  version  4.3A  by  Swanson  Analysis 
Systems  Inc.  The  results  arc  displayed  both  graphically  and  in  numerical  form.  It  is  con¬ 


cluded  that  this  is  an  effective  method  of  tailoring  and  analysis. 
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Chapter  One 

Introduction 

Aeroeiastic  tailoring  is  defined  as  "the  incorporation  of  directional  stiffness  into  an 
aircraft  structural  design  to  control  aeroeiastic  deformation,  static  or  dynamic,  in  such  a 
fashion  as  to  affect  the  aerodynamic  and  structural  performance  of  that  aircraft  in  a  bene¬ 
ficial  way"  [Ij.  In  recent  years,  numerous  efforts  have  dealt  with  the  effects  of  anisotro¬ 
pic  design  on  the  aeroelasticity  or  deformation  coupling  of  fixed  and  swept  wing  aircraft, 
but  very  little  has  been  done  io  research  its  effects  on  rotorcraft  and  specifically  rotorsys- 
tems.  Helicopter  rotor  blades  r)perate  in  an  aeroeiastic  environment  consisting  of  inertial, 
aerodynamic  and  elastic  loadings.  This  environm  ,  n  is  ideal  for  anisotropic  design  and 
improvement  and  provides  for  the  creative  use  of  directional  stiffness  characteristics. 

Composite  material  sys'ems  are  now  the  primary  materials  for  helicopter  rotor  sys¬ 
tem  applications.  Currently,  attention  is  being  focused  on  the  use  of  composites  for 
designing  rotor  blades  in  which  the  elastic  coupling  associated  with  unbalanced  ply  layup 
is  employed  to  enhance  the  dynamic  and  aerodynamic  characteristics.  The  type  of  rotor 
and  its  control  determine  largely  the  aeroeiastic  behavior  of  a  helicopter.  Of  special 
interest  nowadays  are  hingeless  and  bearingless  rotor  systems.  The  LHX  helicopter  pro¬ 
ject,  with  its  bearingless  rotor  designs,  is  a  prime  example. 

There  are  a  number  v)f  aeroeiastic  phenomena  associated  with  the  design  of  hel¬ 
icopters.  The  dynamic  stability  and  response  characteristics  of  rotary-wing  aircraft  are 
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dependent  on  parameters  which  have  to  be  defined  in  the  preliminary  design  phase.  In 
order  to  achieve  maximum  stmctural  efficiency  through  aeroelastic  tailoring,  it  is  neces¬ 
sary  to  use  composite  laminates  with  different  ply  angles.  As  a  result,  composite  beams 
will  exhibit  coupling  among  extensional  stiffness,  bending  stiffness  and  torsional  stiff¬ 
ness.  For  example,  bending-torsional  stiffness  coupling  is  necessary  for  pitch-flap  stabil¬ 
ity  of  rotor  blades,  and  extensional-torsional  stiffness  coupling  is  desirable  to  change  the 
linear  twist  distribution.  Even  isotropic  beams  with  pretwist  will  exhibit  coupling 
between  extension  and  torsion  [2], 

Aeroelastic  tailoring  of  a  composite  structure  involves  a  design  process  in  which  the 
materials  and  dimensions  are  selected  to  yield  specific  torsional  characteristics.  Due  to 
the  directional  nature  of  the  composite  rnaierials,  it  is  possible  to  construct  blades  with 
different  ply  orientations  and  laminae  thickne.'jses.  These  laminations  can  be  either  sym¬ 
metric  or  asymmetric,  where  "symmetric"  means  that  a  layer  of  material  at  some  distance 
above  a  structural  midsurface  reference  location  has  the  iderdcal  ply  thickness,  angular 
orientation,  and  material  properties  as  that  of  a  lamina  at  an  identical  distance  below  the 
midsurface  [1].  If  the  fibers  are  placed  off-axis  in  the  upper  and  lower  portions  of  Lhe- 
blade,  a  twist  nhould  be  induced  when  exerting  a  non-torsional  force.  This  provides  a 
potential  for  improving  the  perforaiance  of  a  lifting  surface  through  aeroelastic  tailoring 
of  the  primary  load-bearing  structure.  The  design  of  such  advanced  structures  requires 
simple  and  reliable  analytical  tools  which  can  take  into  consideration  the  directional 
nature  of  these  materials. 

The  analysis  of  helicopter  rotor  blades  has  traditionally  depended  on  classical  beam 
or  shell  theory.  This  is  inadequate  for  the  dynamic  characteristics  of  orthotropic  and 
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anisoaropically  laminated  shells  which  can  be  thought  of  as  rudimentary  rotor  blade 
structures.  While  analyses  based  on  beam-like  representations  are  generallj  feasible 
methods  of  structural  analyses  for  many  composite  rotor  blade  structures,  the  design  of 
advanced  rotor  blades  will  often  require  a  detailed  shell- finite-element  representation  [3]. 

Finite  element  analysis  has  become  an  established  and  widely  used  tool  in  the 
design  of  rotorcraft  airframes,  kecendy,  however,  research  effons  have  indicated  a 
growing  interest  in  applying  these  analysis  techniques  in  composite  rotor  blade  structural 
analysis. 

Most  of  the  previous  work  in  this  area  involves  the  feasibility  of  modeling  compo¬ 
site  helicopter  blades  using  finite  element  analysis  and  the  treatment  of  the  flexbeam  or 
spar  as  the  model  base  [2-5].  From  this  point  investigators  have  attempted  to  determine 
the  structural  properties  by  theoretical  and  experimental  methods.  Although  some  work 
has  been  done  on  composite  rotor  design  parameters  that  affect  aeroe’astic  stability  and 
response,  it  mainly  centers  around  arbitrary  cross-sectional  warping  and  vibration 
response.  It  is  uncertain  if  anyone  has  attempted  to  vary  the  lamination  sequence  and 
angular  orientation  to  maninnlate  the  classical  beam  equations. 

In  the  present  research,  a  composite  rotor  blade  is  modeled  as  a  hollow  beam  with 
end  caps.  It  can  be  shown  that  the  composite  rotor  blade  concept  offers  the  chance  to 
accomodate  aeroelastic  coupling  effects  which  may  improve  or  deteriorate  the  response 
behavior  of  the  helicopter.  Using  the  commercial  package  PATRAN®,  a  finite  element 
model  of  the  constructed  blade  is  created  and  an  evaluation  of  the  blade  sectional  proper¬ 
ties  is  conducted.  Four  individual  loading  conditions  are  applied  to  eight  distinct  blade 
laminations  to  achieve  the  desired  reactions.  The  results  of  the  analysis  using  ANSYS® 
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are  applied  to  solve  the  classical  beam  equations. 

Following  a  background  chapter  on  the  history  of  aeroelastic  tailoring  and  a  sum¬ 
mary  of  the  work  already  accomplished  in  this  area,  a  concise  but  complete  chapter  on 
composite  materials  which  includes  lamination  theory  is  provided.  This  has  the  advan¬ 
tage  of  building  a  firm  basis  for  the  later  development  of  the  model  and  the  significance 
of  the  numerical  results.  The  classical  beam  equations  are  then  developed  in  Chapter  4 
and  their  validity  for  use  with  composite  material  beams  is  demonstrated  through  lam¬ 
inate  analysis  theory.  Chapter  5  involves  the  modeling  process  with  the  ensuing  chapters 
reporting  the  results  and  conclusions.  The  readers  may  find  the  appendices  useful  for 
similar  work  involving  PATRAN®  and  ANSYS®.  The  FORTRAN  code  for  laminate 

Id  aidu  IiiCiUuc>Q  aS  au  appciiuiA.* 

One  aspect  of  this  research  which  the  author  found  both  frustrating  and  enriching 
was  the  interfacing  of  the  pre-  and  post-processors  and  the  lessons  learned  when  dealing 
with  large  commercial  packages  whose  companies  are  clearly  competitors. 
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Chapter  Two 

Background 

M.H.  Shirk,  TJ.  Hertz  and  T.A.  Weisshaar  describe  the  historical  evolution  of 
aeroelastic  tailoring  and  cite  significant  examples  of  the  early  use  of  tailoring  principles 
[6],  These  consider  primarily  the  well-known  effects  of  the  bending  of  a  sweptback  wing 
and  the  effective  angles  of  attack  which  it  induces.  In  his  review  paper,  "The  First  Fifty 
Years  of  Aeroelasticity"  f?],  A.  R.  Collar  mentions  the  application  of  this  bending 
induced  load  relief  principle  to  an  unswept  wing,  where  a  downward  bending  would 
cause  increased  lift. 

In  a  recent  review  article  [8],  H.  Ashley  and  his  co-authors  note  that,  at  the  incep¬ 
tion  of  manned  powered  flight,  lifting  surface  flexibility  was  used  to  generate  favorable 
aerodynamic  loads.  The  biplane  structural  truss  of  the  original  Wright  Flyer  exhibited 
wing  torsional  flexibility  which  allowed  a  control  cradle,  operated  by  the  prone  pilot’s 
hip  movements,  to  provide  differential  twist  to  the  biplane  wings.  The  result  was  a  rol¬ 
ling  moment.  While  the  use  of  stnictural  deformability  was  quickly  replaced  by  mechan¬ 
ical  devices,  the  Wright  Flyer  ^lovided  what  might  be  regarded  as  the  world’s  first 
aeroservoelastic  system. 

Closer  to  home.  Professor  Max  Munk,  while  at  the  Catholic  University  of  America 
in  Washington,  D.C.,  used  aeroelastic  deformation  to  provide  favorable  propeller  perfor¬ 
mance.  This  design,  patented  in  1949,  was  entitled  "Propeller  Containing  Diagonally 
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Disposed  Fibr  'Us  Material".  The  fibrous  material  was  wood.  The  grain  of  the  wood  was 
oriented  to  provide  spanwise  anisotropy  such  that  bending  due  to  lift  caused  the  propeller 
to  twist  nose-down  in  such  a  way  that  some  lift,  in  this  case  propeller  thrust,  was 
decreased  at  high  speeds.  As  a  result,  an  automatic,  passive,  favorable  pitch  change  in 
the  propeller  was  achieved  without  the  aid  of  a  mechanical  device.  Wind  tunnel  tests 
verified  predictions  of  increased  propeller  performance  at  high  speed  as  a  result  of  this 
tailored  grain  orientation  [9]. 

Although  ingenious,  these  first  designs  had  minor  repercussions  in  the  aircraft 
industry.  This  was  primarily  due  to  a  lack  of  elastic  materials  other  than  wood  and  a  lack 
of  understanding  of  the  true  need  to  incorporate  this  flexibility  in  the  designs.  However, 
this  was  a  step  toward  its  use  for  helicopter  blade  design  and  these  problems  would 
disappear  with  the  introduction  of  composite  materials. 

In  1947,  Cornell  Aero  Lab  built  what  were  probably  the  first  composite  rotor  blades 
which  used  some  fiberglass  reinforced  plastic  (FRP)  construction.  These  blades  had 
wood  spars  with  FRP  skins.  They  flew  on  a  Sikorsky  R-5.  In  1953,  Glenview  Metal 
Products  built  a  GMP-2  Flyride  helicopter  with  a  main  rotor  blade  of  laminated  spruce 
forward  of  the  30%  chord  and  FRP  skins  on  a  balsa  wood  core  aft  of  the  30%  chord.  In 
1956,  Prewitt  Aircraft  Company  built  three  sets  of  researcti  rotors  for  the  Piasecki  HUP-2 
aircraft.  One  set  was  made  of  stainless  steel,  one  of  titanium,  and  one  of  FRP.  The  com¬ 
posite  blade  had  a  wood  spar.  The  Gyrodine  XRON-1  introduced  in  1956  had  all¬ 
fiberglass  blades  on  some  later  versions.  Parsons  also  built  two  different  experimental 
blades  for  the  H-21  from  stainless  steel  and  FRP  in  the  same  period.  One  of  the  compo¬ 
site  blades  failed  in  the  whirl  test  stand  and  was  never  flight  tested  [10]. 
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The  early  1960’s  saw  the  introduction  of  several  composite  rotor  blade  designs. 
Among  these  were  the  Kaman  H43-B  blade  of  1960  which  had  a  wood  core  with  FRP 
skins  on  the  afterbodies.  The  BoUcow  BO-103  helicopter  which  first  flew  in  1961  had  a 
single  main  rotor  blade  made  of  FRP.  The  Boeing  Vertol  CH-47  blade  of  1961  consisted 
of  a  steel  D  spar  with  aluminum  ribs  and  FRP  skin.*!.  In  1962,  Kaman  flew  an  all-FRP 
blade  on  the  HH-43B  helicopter  [10]. 

It  was  the  year  1968  that  marked  the  beginning  of  the  all-boron  Advanced 
Geometry  Blade  program  for  the  Boeing  Vertol  CH-47.  An  all-FRP  version  of  this  blade 
was  designed  and  built  in  the  same  time  period.  Both  of  these  blades  were  flight  tested  in 
the  early  1970’s.  By  this  time,  the  future  of  helicopter  rotor  blades  was  clearly  pointed  in 
the  direction  of  all-composite  construction  and  would  progress  to  the  advanced  graphite- 
epoxy  blades  of  the  1980’s  [10]. 

By  1982,  at  least  50  different  all-composite  main  rotor  blades  had  been  designed. 
A’  out  25  of  those  had  been  flight  tested  and  over  15  had  progressed  to  limited  or  full 
production  status.  The  US  Army  had  approximately  25  different  composite  rotor  blade 
programs  underway  at  this  time  of  which  six  had  reached  flight  status.  These  six  were 
[10]: 

1)  the  CH-47  Advanced  Geometry  Blade  (1972) 

2)  the  AH-1  Multi  Tubular  Spar  Blade  (1975) 

3)  the  CH-47D  Composite  Blade  (1978) 

4)  the  AH-lS  Improved  Main  Rotor  Blade  (1980) 

5)  the  AH-64  Prototype  Blade  (1981) 

6)  the  BO-105  Bearingless  Main  Rotor  (1981) 
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Helicopters  are  usually  classified  according  to  the  mechanical  arrangement  of  the 
hub  in  order  to  accommodate  the  blade  flap  and  lead-lag  motion  ant.  .ccording  to  the 
blade  and  hub  bending  stiffness.  The  most  widespread  helicopter  configuration  uses  a 
single  main  rotor  and  a  small  tail  rotor.  During  the  last  two  decades,  the  hingeless  rotor 
concept  and  its  successor,  the  bearingless  rotor,  have  found  continuosly  growing  interest 
among  helicopter  manufacturers  and  research  organization. 

With  the  evolution  of  advanced  composites,  the  feasibility  of  designing  rotor  sys¬ 
tems  for  high  speed,  dema’iding  maneuver  envelopes,  and  high  aircraft  gross  weights  had 
become  a  reality.  The  desire  to  tailor  the  elastic  qualities  of  these  blades  was  soon  to  fol¬ 
low. 
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Chapter  Three 

Composites 

Advanced  fiber-reinforced  composite  materials  combine  vastly  superior  specific 
stiffness  and  strength  characterisdes  (i.e.  stiffness/unit  weight  and  strength/unit  length) 
with  reduced  weight  and  increased  fatigue  life.  These  improved  fatigue  characteristics 
are  very  significant  to  helicopter  design  because  fatigue  is  probably  the  primary  cause  of 
rotor  blade  failure.  Composites  make  possible  designs  with  fewer  parts  and  improved 
maintainability.  In  terms  of  manufacturing,  it  is  possible  to  achieve  more  general  aero¬ 
dynamic  shapes  such  as  flapwise  variation  in  planform,  section  and  thickness  [11].  Com¬ 
posite  materials  are  ideal  for  structural  applications  where  high  strength-to-weight  and 
stiffness-to-weight  ratios  arc  required.  The  adjective  "advanced"  in  advanced  fiber- 
reinforced  composite  materials  is  used  to  distinguish  composites  with  new  ultraliigh 
strength  and  stiffness  fibers  such  as  boron  and  graphite  from  some  of  the  more  familiar 
fibers  such  as  glass. 

To  understand  fully  the  significance  of  composite  materials,  it  is  necessary  to  start 
with  the  very  basic  element  or  building  block  which  is  the  lamina.  A  lamina  is  an  assem¬ 
blage  of  fibers  oriented  in  the  same  direction  in  a  supporting  matrix.  The  fibers  supply 
most  of  the  strength  and  stiffness  characteristics  and  the  matrix  holds  these  fibers  in  the 
desued  position.  Advanced  composites  are  made  by  laminating  individual  laminae  in 
specific  directions  in  order  to  obtain  the  desired  properties  for  the  given  structural  com- 
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poiicnL  Therefore,  the  knowledge  of  the  mechanical  behavf.or  of  a  lamina  is  essential  to 
the  understanding  of  laminated  composites. 

There  are  several  assumptions  which  are  made  about  a  lamina  and  its  components: 

(a)  The  fibers  are  homogeneous,  linearly  elastic  and  isotropic.  They  are  equally 
spaced  and  perfectly  aligned. 

(b)  The  matrix  is  homogeneous,  linearly  elastic  and  isotropic. 

These  assumptions  then  allow  us  to  say  that  a  lamina  is  macroscopically  homogene¬ 
ous,  linearly  elastic,  macroscopically  orthotropic  or  transversely  isotropic  and  initially 
stress  free. 

The  directional  properties  of  the  lamina  are  the  key  to  the  limitless  uses  of  compo¬ 
site  materials.  In  the  fiber  direction,  the  lamina  is  very  stiff  ajid  strong  while  in  the  direc¬ 
tion  which  is  transverse  to  the  fibers,  the  lamina  is  compliant  and  weak.  The  difference 
in  the  strength  in  these  directions  is  usually  at  least  one  order  of  magnitude;  it  is  not 
uncommon  for  the  lamina  to  be  100  times  stronger  in  the  fiber  direction.  Figure  3.1  [12] 
helps  to  illustrate  the  directional  properties  of  a  lamina.  The  longimdinal  direction  is 
represented  by  0  =  0°  and  the  transverse  direction  by  0  =  90°. 
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WGURH3.1.  Directional  Properties  of  Lamina 

The  description  of  the  mechanical  behavior  of  a  composite  lamina  is  based  on  stress- 
strain  relations  referred  to  as  Hooke’s  Law.  To  model  the  stress-strain  response  of  a  lam¬ 
ina,  we  use  the  generalized  version  for  an  anisotropic  body.  These  reladons  have  two 
commonly  accepted  manners  of  expression;  compliances  and  stiffnesserj  as  coefficients 
of  the  stress-strain  relations.  Jones  [12]  provides  a  description  of  these  relations.  The 
generalized  Hooke’s  law  relating  stresses  to  strains  can  be  written  in  contracted  notation 
as 


Oj  —  CfjF.j  i,j  —  1,...,6 

where  are  the  stress  components,  Cij  is  the  stiffness  matrix,  and  are  the  strain  com¬ 
ponents.  The  stiffness  and  compliance  components  are  often  referred  to  as  elastic  con¬ 
stants.  Both  sets  of  components  can  be  represented  by  the  engineering  constant.' 
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Ei,  Vij,  Gij  which  are  Young’s  moduli,  Poisson’s  ratio,  and  the  shear  moduli  respec¬ 
tively.  The  engineering  constants  are  particularly  helpful  in  describing  lamina  behavior. 
These  relations  reduce  to  21  independent  constants  and  can  be  represented  by  Equation 
(3.1).  These  relations  represent  an  anisotropic  material  since  there  are  no  planes  of  sym¬ 
metry  for  tiic  material  properties.  This  is  the  case  in  which  none  of  the  axes  of  the  lam¬ 
ina  line  up  with  the  principal  axes  of  materia*  symmetry. 
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If  there  is  one  plane  of  material  property  symmetry,  the  stress-strain  relations  reduce 
to  Equation  (3.2)  which  is  termed  monoclinic.  This  relationship  represents  a  lamina  in  a 
coordinate  system  which  is  rotated  in  its  plane. 
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If  there  are  two  orthogonal  planes  of  material  property  symmetry,  the  stress-strain 
relations  are  shown  in  Equation  (3.3)  and  arc  said  to  define  an  orthotropic  material.  Note 
that  there  is  no  interaction  between  nomial  stresses  O]  ,  02,  o?  and  shearing  strains 
Tto.  Yi2  as  in  anisotropic  materials. 
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If  at  every  point  of  a  material,  there  is  one  plane  in  which  the  mechanical  properties 
are  equal  in  all  directions,  then  the  material  is  termed  transversely  isotropic  and  has  the 
relations  shown  in  Equation  (3.4).  This  relationship  expresses  the  response  of  a  lamina 
in  the  principal  material  coordinate  system  (PMCS).  The  PMCS  is  the  coordinate  system 
in  which  the  axes  are  aligned  with  respect  to  the  structural  material  symmetry. 
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if  there  are  an  infinite  number  of  planes  of  material  property  symmetry,  then  the 
material  is  isotropic  with  only  two  independent  constants  and  Equation  (3.4)  becomes: 


^  • 

Cii 

C12 

Ci3  0 

0 

0 

El 

02 

C12 

^22 

C  -Q  0 

0 

0 

f2 

Cfs 

Cr, 

C23 

C3.,  0 

0 

0 

E3 

T23 

*  ■= 

0 

0 

0  2(Ci,-C,2) 

0 

0 

4 

Ya 

X31 

0 

0 

0  0 

2(C„-C,2) 

0 

Ysi 

'h2 

0 

0 

0  0 

0 

2(C  11  -  C 12 ) 

Yi2 

u 


14 


For  plane  stress  in  the  plane  of  the  lamina.  Equation  (3.4)  can  be  rewritten  in  terms 
of  the  reduced  stiffness  matrix  [Q\.  Setting  03,  T23,  and  X31  equal  to  zero,  \  'hich  implies 
that  £3,  Tto*  zero  also,  leads  to  the  following  relationship  between  the  inplane 

stress  components  and  the  inplane  strain  components: 


Qvl  Qn  0 

*  ■ 

£t 

02 

’  = 

(2 12  G22  0 

£2 

ti2 

0 

0 

to 

Y12 

»  J 

(3.6) 


Where  the  elements  of  the  reduced  stiffness  matrix  [Q]  are  given  in  terms  of  tlie  elements 
of  the  3— D  stiffness  matrix  [C]. 


Equation  3.6  defines  inplane  stresses  and  strains  in  the  principal  material  coordinate 
system  of  an  orthotropic  or  transvercsly  isotropic  material.  However,  the  principal  direc¬ 
tions  of  orthotropy  often  do  not  coincide  with  coordinate  directions  that  are  geometri¬ 
cally  natural  to  the  structure  elements.  Therefore,  a  coordinate  transfomiation  is  neces¬ 
sary  in  order  to  relate  stresses  and  strains  in  a  coordinate  system  not  aligned  with  the 
PMCS.  Suppose  that  a  layer  is  rotated  through  an  angle  0  about  the  z  axis  as  in  Figure 
3.2.  We  can  express  its  constitutive  equations  in  the  rotated,  that  is  x-y  coordinate  sys¬ 
tem  by  proceeding  as  follows. 
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FIGURE  3.2.  Rotation  of  Coordinate  System 


If  we  let 


c  =  COS0,  s  =  sin0 


then  the  inplane  stress  and  strain  components  in  the  principal  material  coordinate  system 
can  be  related  to  the  inplane  stress  and  strain  components  in  the  rotated  coordinate  sys¬ 
tem  through  the  transformation  matrices  [Ti]  and  [Ti]  defined  below: 


Ics 

[T’l]^  -2c5 

2  2 

-cs  cs 
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[72]  = 


2  2 

C  S  CS 

2  2 
^  C 

2cj 


(3.8) 


Using  the  above  relations  in  Equation  (3.6)  allows  us  to  relate  the  inplane  stresses  to  the 
strains  in  the  rotated  coordinate  system  in  the  following  fashion: 


cr'  =  [G]e'  or  e'  =  [S]CT' 

where 


(3.9) 


[I2]  =  [7’ir’[<2][r2]  (3.10) 

and 

[S]  =  [T2r^[S][Ti]  (3.11) 

With  the  other  benefits  of  composite  materials,  comes  the  ability  to  design  or 
prescribe  laminate  geometry.  As  mentioned  previously,  laminated  composites  consist  of 
differently  oriented  layers  or  laminae  that  are  bonded  together. 

It  is  important  to  understand  the  motivation  for  laminating  the  plies  discussed 

above.  Because  each  individual  ply  has  very  directional  properties  (Figure  3.1)  we  lam¬ 

inate  to  combine  the  best  aspects  of  the  constituent  layers  in  order  to  achieve  a  more  use¬ 
ful  material.  This  is  done  by  tailoring  the  directional  dependence  of  strength  and  stiff¬ 
ness  of  a  material  to  match  the  loading  environment  of  the  structural  element.  For  exam¬ 
ple,  a  structure  that  has  a  dominant  uniaxial  loading  in  the  fiber  direction,  such  as  the 
centrifugal  force  on  a  rotor  blade,  would  need  a  predominantly  unidirectional  composite. 
For  more  general  lof-  Uuations,  the  principal  material  directions  are  oriented  to 
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prodiicot  a  Ssructural  element  capable  of  resisting  load  in  several  directions.  Figure  3.3 
prcx-iv  5S  an  e:f.ample  of  layering  with  different  ply  orientations. 


FIGURE  3.3.  Laminate  Construction 


Laminated  composites  are  generally  strong  and  stiff  in  the  plane  of  lamination  and 
weak  and  flexible  in  the  tiansverse  direction.  This  stiffness  is  obtained  Irom  the  proper- 
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ties  of  the  conistitucnt  laminae.  Laminated  plates  are  one  of  the  simplest  and  most 
widespread  practical  applications  of  composite  laminates.  They  will  be  used  in  this  work 
as  they  are  uniquely  suited  to  its  objective.  Therefore,  a  summary  of  lamination  theory  is 
provided. 

Lamination  theory  allows  us  to  predict  the  ensuing  properties  of  a  laminated  plate  or 
straemre  when  the  properties,  dimensions  and  stacking  sequence  of  individual  plies  are 
specified.  Figure  3.4  [12]  displays  the  geometry  of  a  layered  laminate. 


^  LAYER  NUMBER 


FIGURE  3.4.  Geometry  of  an  N-layered  Laminate 


There  arc  two  static  parameters  which  are  used  to  describe  the  response  of  a  lam- 


inate.  The  first  is  in-plane  force  resultants  defined  as: 
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••  n 

H 

=  j  '  Oy  ’  dz 
-H 


(3.12) 


The  second  is  the  in-plane  moment  resultants  defined  as: 


n  o* 

My  '  =  f '  Oj,  -  zdz 
-H  _ 


(3.13) 


The  kinematic  parameters  which  are  used  to  describe  the  response  of  a  laminate  are 
based  on  the  following  assum])tions: 

(a)  The  laminate  consists  of  perfectly  bonded  laminae,  each  of  which  are  in  a  state  of 
plane  stress. 


fill  fii2  fii6  c, 
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(b)  The  normals  to  the  midplane  remain  straight  and  normal  to  the  deformed  midplane. 

(c)  The  normals  do  not  change  in  length. 

The  last  two  assumptions  are  a  consequence  of  the  Kirchoff-Love  hypothesis.  Tliis 


yields: 


du  _  dw 
dz  dx 


(3.15) 


dv  _  dw 
dz  dy 


(3.16) 
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Also,  by  expressing  th^  inplane  displacement  con^xjnents  in  the  following  fashion: 


u(x,y,z)=^u\x,y)  +  rf(x,y)  (3.17) 

and 


v(x,y,z)  =  v°(x,y)  +  zg(x,y)  (3.18) 

which  is  consistent  with  the  Kirchoff-Love  hypothesis,  it  can  be  shown  that  the  in-p)ane 
strains  are  a  function  of  the  in-plane  strains  of  the  midplane  and  the  curvatures  (  k  )  of 
the  midplane: 


e* 

•  ^ 
^  X 

-  " 

Kx 

% 

• 

e^y 

•  +z- 

K, 

Yxy 

7°jcy 

^  J 

»  - 

where 


Combining  this  equation  with  Equation  (3.14)  yields: 

k  r, 

Ox-y  =  Q  e  +  zK 


dx^ 

dxdy 


(3.19) 


n  9m 


(3.21) 


By  integrating  the  equation  for  in-plane  forces  with  this  equation,  we  arrive  at  the  follow¬ 
ing  result: 
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N  =Ae“  +  5K 


where  is  the  cxtcnsional  stiffness  laminate  matrix  which  is  written  as: 


(3.22) 


~  2  Qij  ~  ^k-l  ) 
Jk=l 


and  where  is  the  extensional  bending  coupling  matrix  which  is  expressed  as: 


(3.23) 


^  *=i 


(3.24) 


This  particular  matrix  will  become  zero  when  a  symmetric  laminate  (see  Chapter  1,  page 
2)  is  used.  By  integrating  the  equation  for  in-plane  moments  with  equation  (1.20)  we 


M  =  Be°  +  DK 


where  is  the  bending  stiffness  matrix  which  is  written  as; 


(3.25) 


=  Qij^n-n-i ) 

^k=l 


(3.26) 


The  bending  stiffness  matrix  is  important  to  this  research  as  it  is  the  factor  which  causes 
the  torsional  motion  of  the  blade.  The  D  j,.;  term  is  the  key  element  of  the  matrix. 

The  tlnal  results  of  this  procedure  axe  the  constitutive  equations  for  laminate 
analysis  theory.  They  aie  expressed  for  all  composite  material  laminates  as: 


i¥r[isji« 


(3.27) 
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Laminates  are  not  without  their  problems,  however,  as  potentially  adverse  shearing 
and  normal  stresses  are  introduced  between  adjacent  layers  in  the  vicinity  of  the  free 
edge.  Such  stresses  may  cause  delamination.  These  stresses  are  out-of-plane  or  inter¬ 
laminar  stresses  which  arise  to  satisfy  equilibrium  requirements  at  the  free  edge. 
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Chapter  Four 


Governing  Equations 

A  thorough  understanding  of  the  mechanical  behavior  of  a  structure  is  essential  for 
its  effective  design.  This  includes  helicopter  rotor  blades.  The  use  of  physical  properties 
of  composite  materials,  as  well  as  numerous  theoretical  laws  and  concepts,  becomes  tlie 
basis  for  this  understanding.  Mechanics  of  materials  is  a  branch  of  applied  mechanics 
that  deals  with  tiie  behavior  of  bodies  subjected  to  various  types  of  loading.  The  objec¬ 
tive  of  this  particular  analysis  will  be  the  determinadon  of  the  deflections  and  displace¬ 
ments  of  a  simple  blade  model  produced  by  these  loadings.  This  data  will  then  be 
applied  to  the  Euler-Bemoulli  beam  equations  to  determine  the  off-diagonal  term 
response  to  the  loading  of  composite  material  lamination  sequences. 

This  research  will  center  around  only  one  finite  element  model,  however  there  will 
be  eight  individual  blade  cases  examined  These  eight  cases  are  concerned  with  seven 
anisotropic  composite  material  lay-ups  and  one  isotropic  lay-up.  Each  of  these  cases  will 
be  subjected  to  identical  loading  conditions  consisting  of  four  separate  loads.  They  con¬ 
sist  of  two  bending  moments,  a  torsional  force  and  an  axial  force,  and  will  be  described 
in  derail  at  a  later  point. 

There  are  certain  properties  of  area  which  are  needed  for  this  beam  theory,  such  as 
the  location  of  the  cc.ntroid,  the  fust  moment  of  the  area  and  the  second  moment  of  the 
area  or  the  moment  of  ,rtia  with  respect  to  each  axis. 
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The  bending  axis  of  the  beam  imder  investigation  is  the  longimdinal  axis  through 
which  the  transverse  bending  loads  must  pass  in  order  that  the  bending  of  the  beam  shall 
not  be  accompanied  by  twisting  (for  isotropic  materials).  The  shear  center  or  center  of 
twist  for  any  transverse  section  of  the  beam  is  the  point  of  intersection  of  the  bendi  ig 
axis  and  the  plane  of  the  transverse  section  [13]. 

The  bending  axis  and  the  shear  center  are  of  special  importance  in  beams  having 
cross  sections  composed  of  thm  parts  which  offer  large  resistance  to  bending  such  as  in 
airplane  or  rotor  blade  construction. 

The  bending  axis  for  a  beam  whose  cross  section  has  two  axes  of  symmetry  is  the 
longitudinal  centroidal  axis  of  the  beam  and  hence  the  shear  center  for  such  a  section  is 
also  the  centroid  of  the  section.  When  conducting  experiments,  for  example,  the  axial 
force  and  constraints  must  act  through  this  centroid  of  the  cross-sectional  area. 

Calculaticn  of  the  moments  of  inertia  (second  moments  of  the  area)  requires  a  sub- 
.stantial  effort  for  complex  geometry  of  this  type.  Using  the  equation 

Iy=\]z'^dA  (4.1) 

A 

and  the  fact  that  the  blade  is  symetric  about  the  y  axis,  we  obtain 

/y=2ffz^dA  (4.2) 

Aa 


where  z  =  rco^p.  See  Figure  4.1.  Substituting  for  z  we  obtain 
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h~\  \  (4.3) 

An 

Since  dA~dst  and  ds  =  rtd^,  then 

n 

/j,  =  2 1  (4.4) 


FIGURE  4;  1 ,  Blade  Moment  of  Inertia 


A  similar  derivation  can  be  completed  to  determine  the  moment  of  inertia  about  the 
z  axis.  Since 

I,  =  jy^dA  (4.5) 

/I 


then 
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Jt 

/j  =  2  |/r^sin^pjp  (4.6) 

Sines  the  radius  r  is  an  arbitrary  function  of  P.  a  numerical  solution  to  this  equation 
is  preferred.  A  jhort  code  was  written  to  divide  one  half  of  the  blade  into  24  elements 
and  to  calculate  the  average  radius  for  each  element  to  the  centroid.  Then  a  simple  sum¬ 
mation  for  each  moment  of  inertia  was  completed. 
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Iy=2YjiSimi  (4.7) 

»=i 

and 
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i=l 


(4.8) 


where  s  is  the  element  length  and  where  m  —  rcosP  and  n  =  /"sinp. 

The  four  leading  conditions  will  now  be  examined.  In  this  problem,  a  rotor  blade 
represented  by  a  hollow  shell  beam  of  length  L  with  end  caps  as  shown  in  Figure  4.2,  will 

iicnH  V  pam  ic  ct/mmptTnr'ol  tKn  ^  ra^c* A  nv«ci  oMy-l  rtv«r« 

ww^-w.  •  -  •  V  w  JT  WK  SUll?  TV*  «JV  lA/XXa  t**I\A  tliV 

runs  along  tlic  length  of  the  blade. 
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FIGURE  4.3.  Axially  Loaded  Beam 


If  the  force  P  acts  at  tliC  centroid  of  the  cross  section,  the  stress  in  the  beam  at  sections 
away  from  the  ends  is  given  by: 


o  = 


A 


(4.9) 


where  A  is  the  cross-sectionai  area.  Additionally,  the  strain  can  be  expressed  as 


e  =  |-  (4.10) 

where  8  is  the  displacement  produced  by  the  axial  force.  Utilizing  Hooke’s  law 
( a  =  £e ),  the  equation  for  beam  displacement  is  given  as: 


8  = 


PL 

EA 


(4.11) 


Since  the  axial  force  varies  along  the  axis  of  the  beam,  an  expression  for  displacement 
along  the  entire  length  of  the  beam  is  obtained.  Also,  we  will  denote  the  displacement  of 
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the  beam  in  the  x  direction  by  u. 

(4.12, 

The  second  force  to  be  examined  is  that  of  torsion.  Torsion  refers  to  the  twisting  of 
a  structural  member  when  loaded  by  couples  that  produce  rotation  about  its  longitudinal 
axis.  In  this  case,  the  x-axis  is  the  longitudinal  axis.  This  type  of  loading  is  pictured  in 
Figure  4.4  which  shows  the  same  beam  now  constrained  at  one  end  and  loaded  at  the 
other  end  with  a  moment  about  the  x  axis. 


z 
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1 
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FIGURE  4.4.  Torsional  Loading 


During  twisting,  there  will  be  a  rotation  about  the  x-axis  of  one  end  of  the  blade 
with  respect  to  the  other.  The  free  end  will  rotate  through  a  small  angle  p  with  respect  to 
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the  fixed  end.  The  angle  P  is  referred  to  as  the  angle  of  twist. 
This  angle  of  twist  can  be  expressed  in  the  form: 


(4.13) 


where  T  is  the  torque  and  GIp  is  the  torsional  rigidity.  A  combination  of  the  polar 
moment  of  inertia,  Ip,  the  shear  modulus,  G,  and  the  length  yields  the  torsional  flexibility 

or  This  expression  is  analagous  to  the  axial  flexibility,  given  in  Equation 
GIp  EA 

(4.11)  [14],  The  shear  modulus  relates  to  the  modulus  of  elasticity  by  the  equation 


E 

2(1 +v) 


(4.14) 


The  term  Ip,  or  the  polar  moment  of  inertia,  is  not  an  ideal  term  for  this  model.  It  is 
best  suited  for  solid  circular  beams.  Thin  walled  beam  theory  describes  a  similar  term,  J, 
or  torsion  constant,  which  will  be  used  where: 


line. 


FIGURE  4.5.  Median  Line  of  Cross  Section 


The  next  loading  condition  will  be  acting  transversely  to  the  longitudinal  axis.  The 
created  load  is  that  of  a  bending  moment.  These  lateral  loads  cause  the  beam  to  bend  or 
flex,  thereby  deforming  the  axis  of  the  beam  into  a  curved  line.  An  illustration  is  given 
in  Figure  4.6  which  shows  the  same  beam  which  is  subjected  to  a  bending  moment  about 
the  transverse  ory  axis. 


/ 

/ 


A 


FIGURE  4.6.  Cantilever  Beam  Under  Bending  Moment 
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After  loading,  the  axis  is  bent  into  a  curve  that  is  known  as  the  deflection  curve  of 
the  beam.  Now  consider  two  points  m  i  and  m2  on  the  deflection  curve  (Figtire  4.7)  [14], 
where  point  mi  is  at  distance  x  from  the  z  axis  and  point  m2  is  situated  a  small  distance 
ds  further  along  the  curve. 


FIGURE  4,7.  Curvature  of  a  Bent  Beam 


At  each  point,  if  we  draw  a  line  normal  to  the  tangent  to  the  deflection  curve,  a 
center  of  curvature  will  occur  at  point  O'.  TIic  length  of  this  normal  is  called  the  radius 
of  curvature  and  is  defined  as  p.  As  defined  in  calculus  and  analytic  geometry,  the  cur- 
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vature  k,  is  the  reciprocal  of  the  radius  of  curvature: 


1 

K=  — 

P 


Also,  from  the  geometry  of  the  figure,  we  obtain 


(4.17) 


pd(^  =  ds  (4.18) 

where  is  the  small  angle  between  the  normals  and  ds  is  the  distance  along  the  curve 
between  the  normals.  If  the  deflections  of  the  beam  are  small,  which  is  the  case  in  this 
research,  then  the  deflection  curve  is  very  flat  and  the  distance  ds  along  the  curve  may  be 
set  equal  to  its  horizontal  projection  dx  and  we  obtain 


1  d(b 

(4.19) 

Using  the  fact  that  the  integral  of  all  elemental  moments  over  the  entire  cross-sectional 
area  must  result  in  the  total  applied  moment 


Mo=-jc;cydA  (4.20) 

and  Hooke’s  law  for  uniaxial  stress 

ax  =  -EKy  (4.21) 

we  obtain 


M=~kEI 


(4.22) 


Rearranging  and  combining  with  equation  (4.17)  we  arive  at  [14] 


34 


and  consequently 


K  = 


p  El 


(4.23) 


3<|)  My 
dx 


(4.24) 


A  similar  bending  moment  about  the  vertical  or  z  axis  yields  an  analagous  equation: 


^\^f  _  Mj 

dx  ~~~Ef 


(4.25) 


The  result  of  each  of  these  loadings  and  the  equations  which  will  be  used  to  analyze 
the  data  can  be  expressed  in  maiix  form  as  follows: 


fie.' 

dx 

0  0  o' 

il 

3x 

^  ^  °  0 

iHL 

3x 

du 

»  0  i, 7  0 

0  0  0  -^ 

a,. 

""J 

(4.26) 


The  off-diagonal  terms  of  Equation  (4.26)  in  general  will  not  be  zero  for  anisotropic 
materials.  This  is  due  to  coupling,  llic  response  of  these  terms  for  different  laminations 
sequences  will  be  the  focus  of  the  present  lesearch.  As  anisotropic  composite  materials 
are  loaded,  these  terms  should  take  on  a  non-zero  identity  and  should  have  varied  reac¬ 
tions  to  different  lay-ups. 


It  is  important  to  note  at  this  point  that  one  can  use  the  displacement  information  for 
this  model  which  is  taken  at  a  sufficient  distance  from  the  ends  as  a  result  of  St.  Venant’s 
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principle.  We  can  say  that  the  effects  of  surface  tractions  over  a  part  of  the  boundary  thai 
are  felt  relatively  far  into  the  interior  of  an  elastic  solid  are  dependent  only  on  the  rigid- 
body  resultant  of  the  applied  tractions  over  this  part  of  the  boundary.  By  this  principle 
we  can  replace  the  complex  supporting  force  distribution  exerted  by  the  constraints  on 
the  beam  by  a  single  force  and  couple.  Although  we  shall  not  present  them  here,  it  is 
pointed  out  that  mathematical  justifications  have  been  advanced  for  St.  Venant’s  princi¬ 
ple  (Goodier,  1937;  Hoff,  1945;  Fung,  1965)  [13]. 

Before  moving  on  to  another  area,  some  mention  of  the  relationship  of  equation 
(4.26)  to  composite  materials  must  be  made.  Since  beam  type  specimens  under  loading 
conditions  are  utilized  in  composite  material  characterization,  a  theory  for  laminated, 
anisotropic  beams  is  also  desirable. 

Beam  bending  is  often  based  on  homogeneous  isotropic  beam  theory.  For  lam¬ 
inated  materials  such  as  the  composite  blade  model  in  this  research,  the  classical  beam 
formulas  must  be  modified  to  account  for  the  stacking  sequence  of  individual  plies. 
Whimey  provides  the  best  description  of  this  modification  [15]. 
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Consider  the  laminated  beam  shown  in  Figure  4.8.  It  has  been  shown  by  Hoff  and 
Pagano  that  layered  beams  of  this  type  in  which  the  plies  arc  oriented  symmetrically 
about  the  midplane  and  the  orthotEopic  axes  of  material  symmetry  in  each  ply  arc  parallel 
to  the  beam  edges  can  be  analyzed  by  the  classical  beam  theory  previously  explained  if 
the  bending  stiffness  El  is  replaced  by  the  equivalent  stiffness  E^I  defined  in  the  follow¬ 
ing  manner: 

£»/  =  £  £f/‘  (4,27) 

k=\ 

where  E^  is  the  effective  bending  modulus  of  the  beam,  £/  is  the  modulus  of  the  /nh 
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layer  relative  to  the  beam  axis  and  N  is  the  number  of  layers  in  the  laminate. 

Equations  which  are  applicable  to  a  general  class  of  symmetric  laminates  can  be 
derived  by  considering  a  beam  as  a  special  case  of  a  laminated  plate  with  a  length  much 
larger  than  the  width,  i.e.,  L:>b. 

For  bending  of  symmetric  laminates,  the  constitutive  relations  (Chapter  3)  reduce  to 
the  form: 


'll  ^12 


Di2  D22  D26  '  iCy 
^16  ^26  ^66 


(4.28) 


whe.n*. 


3x^  dxdy 


(4.29) 


For  this  derivation,  it  is  useful  to  consider  Equation  (4,28)  in  the  inverted  form 


1C_  D\^  D\2  ^16 

'~J¥  -  • 

Ky  ►  =  0*12  ^22  ■^26  " 

^16  ^26  ^66_ 


(4.30) 


where  Dij  are  elements  of  the  inverse  matrix  of  Dij. 


In  order  to  derive  a  beam  theory  the  following  assumptions  are  made: 


My  ~  Mxy  =  0 


(4.31) 


Using  Equations  (4.29)  and  (4.30)  in  conjuction  with  Equation  (4.31),  we  find 
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a^w 


=DnMx 


(4.32) 


Since  this  beam  is  assumed  to  have  a  high  length-to-width  ratio,  we  can  state  that 


w  =  w(x) 


(4.33) 


Combining  Equations  (4.32)  and  (4.33),  we  obtain  the  following  result: 


d^w  _  M 
dx^  ~  E^I 

where 


(4.34) 


12 

h^Dn 


,  M  =  hA4,  /  = 


. ,  1 
on' 

12 


(4.35) 


and  b  is  the  width  of  the  beam.  Equation  (4.34)  is  in  the  same  form  as  classical  beam 
theory  with  the  homogeneous,  isotropic  modulus  E  replaced  by  the  effective  bending 
modulus  of  the  laminated  beam,  Ex  [13].  This  composite  beam  analysis  allows  an 
elastically-coupled  composite  structure  to  be  described  in  terms  of  its  equivalent 
engineering  properties  {EA,  El,  GJ). 
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Chapter  Five 


The  Model 


5.1.  Finite  Element  Modeling 

The  first  step  toward  the  complete  analysis  of  an  aeroelastic  problem  is  the  develop¬ 
ment  of  an  accurate  geometric  and  structural  model.  This  chapter  summarizes  that 
development  process.  A  theoretical  rotor  blade  model  based  upon  composite  construc¬ 
tion  is  used  in  this  research.  Certain  stmctural  elements,  such  as  spars  and  helicopter 
rotor  blades,  can  be  approximated  as  beams.  Rotor  blades  almost  always  involve  compli¬ 
cated  cross-sectional  geometries,  and  combined  axial,  bending  and  torsional  loadings.  A 
major  obstacle  in  applying  this  model  toward  understanding  the  effects  of  composite 
materials  on  aeroelastic  tailoring  is  that  analytical  tools  for  such  a  purpose  are  limited. 

The  finite  element  method  was  chosen  for  this  work  as  it  provides  a  basis  of  modem 
structural  analysis.  It  is  a  powerful  method  to  treat  nonuniformities  and  complex 
geometries  such  as  a  twisting  helicopter  blade.  Finite  element  analysis  has  been  used 
previously  (Hong,  1985)  to  examine  the  aeroelastic  stability  of  helicopter  blades  in  both 
hovering  and  forward  flight  [16]. 

As  stated  earlier,  finite  element  analysis  has  become  an  established  and  widely  used 
tool  in  the  design  of  rotorcraft  and  more  recently  in  composite  rotor  blade  analysis.  In 
addition  to  multi-purpose  finite  element  codes  which  have  composite  material  capability, 


composiie  plate  and  shell  finite  elements  have  been  formulated.  Numerous  research 
efforts  have  resulted  in  the  development  of  shell  elements  with  anisotropic  laminate 
capabilities.  The  result  of  these  efforts  are  found  in  numerous  large  commercially  avail¬ 
able  packages  such  as  PATRAN*  which  can  be  used  for  comprehensive  analytical 
representations  using  detailed  shell  finite  element  models. 

The  specific  advantages  of  a  finite  element  solution  procedure  are  as  follows  [4]: 

(a)  The  formal  derivation  of  the  complex  nonlinear  equations  of  motion  of  the 
problem  is  not  required. 

(b)  Each  of  the  nonlinear  terms  is  dealt  with  in  a  rational  fashion  while  bypassing 
the  need  for  an  ordering  scheme. 

(c)  The  complex  structural  behavior  of  the  blade  is  accurately  modeled. 

(d)  Both  the  undefomied  and  deformed  geometry  of  the  blade  are  taken  into 
account  in  a  natural  fashion  and  ate  available  for  graphics  viewing. 

For  the  analysis  of  a  slender  rotor  blaxlc,  a  beam  element  seems  most  appropriate 
since  only  a  small  number  of  degrees  of  freedom  are  required.  However,  the  accuracy  of 
a  simple  Eultr-Bemoulli  beam  model  for  this  losearch  has  been  questioned  in  Chapter  3. 
This  problem  J  “'’’sfore  requires  a  two-dirnensionai  laminated  sheii  element  to  model 
clastic  couplings  of  anisotropic  composite  materials. 

In  the  present  work,  the  composite  static  blade  is  treated  as  a  hollrv  semi- 
cantilevcrcd  laminated  shell  beam  with  end  caps  as  shown  in  Figure  5.1.1.  The  model  is 
representative  of  the  size  of  a  main  rotor  blade  in  order  to  assess  the  importance  of  ply 
angle  orientation,  'fhe  beam  section  is  assumed  to  be  a  single  airfoil  shaped  cell  consist¬ 
ing  of  six  sections  which  represent  laminates;  two  on  top,  two  on  bottom  and  two  as  end 
caps,  rigures  5.1.2  and  5.2  provide  a  cross-section  J  view  of  the  end  cap  and  laminate 
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section  view  respectively. 

The  solid  modeling  is  done  on  PATRAN*  by  using  a  series  of  grids,  lines  and 
patches.  PATRAN*  is  an  interactive  engineering  program  with  the  ability  to  construct, 
view,  verify,  analyze,  manipulate,  and  demonstrate  [17].  The  four  surface  laminates  and 
two  end  caps  were  constructed  as  six  individual  patches.  Each  patch  was  assigned  a 
material  property  identification.  This  property  was  based  on  the  construction  of  a 
number  of  composite  layers  with  specified  angle  orientations  and  thicknesses.  The 
material  property  option  of  LAM  or  LAMS  in  PATRAN*  was  used  to  specify  the 
number  of  plies,  their  individual  thicknesses  and  their  orientations.  A  desaiption  of  all 
of  the  commands  for  model  generation  is  given  as  a  session  file  in  Appendix  B. 

The  constitutive  relations  of  an  orthotiopic  graphite  fiber  (T50  Fiber)  and  isotropic 
matrix  (Fiberite  934  Epoxy  Resin)  were  used.  The  volume  mix  of  fiber  to  matrix  used 
for  the  composite  was  60%  to  40%  respectively.  The  properties  can  be  found  in  Table 
5.1  below. 


1  MATERIAL  PROPELTIES 

Property 

(Units) 

T50  Fiber 

Epoxy  Matrix 

El 

(Psi) 

56.3 

0.674 

Ef 

(Psi) 

1.10 

0.674 

'^LT 

— 

0.41 

0.363 

Vtt 

— 

0.45 

0.363 

Glt 

(Psi) 

2.20 

0.247 

Gtt 

(Psi) 

0.50 

0.247 

wl 

{\xinlin  F) 

0.50 

25.0 

Ctj 

i\un  /in  F) 

6,50 

25.0 

TABLE  5.1 
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The  resultant  graphite/cpoxy  composite  properties  are  shown  in  Table  5.2. 


Graphite/Epoxy 

Property 

(Units) 

Gr/Ep 

Ex 

(Psi) 

19.2  X  10® 

El 

(Psi) 

1.56  X  10® 

El 

(Psi) 

1.56  X  10® 

Vi2 

— 

1.95  a:  10“^ 

V23 

— 

.21 

Vi3 

.801 X  10“2 

0x1 

(Psi) 

8.2  X  10® 

TABLE  5.2 


The  analysis  code  used  was  ANSYS*  by  Swanson  Analysis  Systems  Inc.  The 
ANSYS*  program  is  a  large-scale,  general  purpose  computer  program  for  the  solution  of 
several  classes  of  engineering  analyses.  It  uses  the  wave-front  direct  solution  method  for 
the  system  of  simultaneous  linear  equations  developed  by  the  matrix  displacement 
method.  The  program  has  the  capability  of  solving  large  structures  with  no  limit  on  the 
number  of  elements  used  in  an  analysis  [18].  There  is  no  "band  width"  limitation  in  the 
analysis  definition;  however,  there  is  a  "wave-front"  restriction  of  2000  degrees  of  free¬ 
dom.  There  are  some  peculiarities  in  ANSYS*  associated  with  the  elements  used  for  this 
laminated  shell  model,  'j^hc  element  is  an  elastic  flat  triangular  clement  with  both  bend¬ 
ing  and  membrane  capabilities  (See  Figure  5.3)  [18].  Up  to  15  different  layers  of 
material  ate  permitted  per  element  and  both  in-plane  and  normal  loads  may  be  applied. 
The  element  has  six  degrees  of  freedom  at  each  node.  There  is  no  significant  stiffness 
associated  with  rotation  about  the  element  z  axis.  A  nominal  value  of  stiffness  may  be 
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present  to  prevent  free  rotation  at  the  node.  The  triangular  2-D  elements  are  less  accurate 
than  equivalent-sized  quadrilateral  elements  but  must  be  used  for  laminated  shell 
analysis.  Wlien  using  triangular  elements  in  a  rectangular  array  of  nodal  points,  the  best 
results  are  obtained  from  an  element  pattern  having  alternating  diagonal  directions. 
Also,  since  the  element  coordinate  system  is  relative  to  the  I-J  line,  the  stress  results  are 
most  easily  interpreted  if  the  I-J  lines  of  the  elements  are  all  parallel. 


z 


FIGURE  5.3.  Laminated  Shell  Element  (ST1F53)  Geometry 


The  model  was  set  up  for  finite  element  analysis  using  flat  two-dimensional  triangu¬ 
lar  elements  constructed  over  tluee  nodes,  TRH3/53,  where  53  is  the  stiffness  coefficient 
corresponding  to  the  laminated  shell  configuration  of  ANSYS*.  A  nine  node  element. 
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TRI/9,  was  picfcncd  but  was  impossible  under  tlie  ANSYS*  restriction  of  three  nodes 
per  element  when  using  laminated  shells.  The  result  was  760  elements  and  310  nodes 
constructed  over  the  entire  blade  region.  In  addition,  the  end  caps  were  assumed  rigid 
and  were  therefore  modeled  with  much  thicker  and  stiffer  elements.  Once  these  finite 
elements  were  defined  using  the  connectivity  sub-program  of  PATRAN*,  the  model  was 
subjected  to  automated  geometric  equivalencing  over  each  of  the  elements.  This  is  sim¬ 
ply  the  process  of  selecting  a  single  ID  for  all  node  points  which  coexist  at  a  given  point 
(and  as  a  result  have  different  ID  numbers)  and  propagating  that  change  through  any 
existing  connectivity  detinitions.  The  next  step  was  the  optimization  of  the  model  by 
selection  of  the  maximum  wavefront  criterion  for  both  the  CuthilUMckee  and  Gibbs- 
Poole-Stockmeyer  methods.  TJic  result  is  element  compaction  based  on  the  assignment 
of  sequential  ID  numbers  to  remove  any  "holes"  or  spaces  in  the  original  numberini, 
sequence. 

Forces,  displacements,  pressures  and  constraints  are  initiated  at  this  point.  It 
required  several  attempts  to  determine  the  best  way  to  constrain  this  model.  The  accu¬ 
racy  of  the  displacement  results  was  very  much  dependent  on  the  constraints  chosen.  An 
example  of  a  poorly  constrained  model  can  be  seen  in  Figure  5.4.  This  particular  model 
was  constrained  (six  degrees  of  freedom)  at  one  node  on  the  leading  edge  of  the  center  of 
the  blade.  An  axial  force  was  applied  at  each  end  in  opposite  directions.  The  result 
should  have  looked  something  like  Figm-e  5.5. 
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nCURE  5.4 

Improperly  Constrained  Model 
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FIGURE  5.5.  Beam  Response  to  Axial  Load 


Finally,  for  simplicity  in  arriving  at  a  solution,  one  end  of  the  blade  model  was  con¬ 
strained.  This  constraint  occurred  at  each  node  on  the  end  cap  in  each  of  its  six  degrees 
of  freedom  (three  translational  and  three  rotational).  This  was  used  to  simulate  the  hub 
attachment  for  the  blade. 

The  other  end  underwent  forces  in  varying  amounts  and  directions.  Each  of  the 
eight  ply  configurations  were  subjected  to  identical  loadings.  Three  of  the  four  loading 
conditions  were  applied  at  a  node  on  the  end  cap  which  represented  the  centroid  (See 
Chapter  4)  of  the  blade  and  consisted  of  an  individually  applied  torsional  load  Tx,  and 
moments  about  the  y  and  z  axes.  My  and  M^.  The  fourth  load  was  an  axial  force  P  and 
was  applied  at  the  neutral  axis  of  the  blade.  Since  there  was  not  a  node  located  on  the 
neutral  axis,  the  force  was  distributed  over  the  nodes  of  the  element  in  which  the  neutral 
axis  was  located. 

All  of  this  infonnation  is  pre-processed  by  PATRAN^  and  compiled  in  a  neutral 
file  in  preparation  for  analysis  by  ANSYS*.  Following  the  analysis,  the  information 
would  normally  go  throi  ,  ■  a  reverse  translator  and  be  post-processed  by  PATRAN*. 
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However,  we  were  unable  to  post-process  in  iPATRAN*  due  to  the  two  commercial  pro¬ 
grams  being  incompatible  on  the  University  of  Virginia’s  current  version  of  PRIMOS 
(Version  21.0.2).  Instead  the  post-processing  was  done  using  POSTl  in  ANSYS43.A*. 
The  results  are  displayed  graphically  for  deformed  and  imdeformed  geometries,  fringe, 
carpet  and  contour  plots,  or  stress  and  displacement  plots.  The  graphics  capability  of  this 
new  version  are  much  improved  over  the  older  version  of  ANSYS43*. 

S2.  Coordinate  System 

The  rotor  blade  is  treated  as  a  nonrouting  (Q  =  0)  elastic  beam.  Therefore,  a  rec¬ 
tangular  coordinate  i-ystem,  x,  y,  z,  is  used  and  attached  to  the  undeformed  blade.  The 
longitudinal  direction  of  the  blade  aligns  with  the  x  axis  and  the  y  axis  corresponds  to  the 
transverse  direction.  The  z  axis  correlates  to  a  right-handed  coordinate  system.  The 
angles  of  roution  about  each  of  the  axes,  (P,  <j>  and  y  for  x,y,z  respectively)  are  shown  in 
Figure  5.6. 
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It  is  assumed  that  the  end  caps  provide  the  required  blade  structural  stiffnesses  in  each 
direction.  Each  laminate  consists  of  a  number  of  laminae  with  ply  orientations  which  are 
different  than  tlie  beam  axis  orientation,  A  point  p  on  the  undeformeu  axis,  undergoes 
displacements  w,  v,  w  in  the  x,  y,  z  directions  respectively  and  occupies  point  p '  on  the 
newly  deformed  axis. 
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5J.  Pty  Configuration 

Two  criteria  were  used  in  the  selection  of  ply  lay-ups  and  orientations  for  each  con- 
riguration.  The  first  required  that  the  blade  did  not  exhibit  static  instability  and  the 
second  was  that  the  result  would  have  coupling  characteristics  which  were  large  enough 
to  have  identifiable  effects.  The  blades  are  eight-ply  graphite-epoxy  laminate  beams 
which  have  strong  bending-twisting  coupling.  Favorable  and  unfavorable  angle-of-attack 
changes  can  be  obtained  depending  on  the  bending-twisdng  coupling  (D 15)  terms.  There 
are  eight  different  lay-ups  used  in  this  research.  Seven  are  anisotropic  configurations 
with  the  angle  orientation  of  two  layers  changing  systematically  to  determine  its  result  as 
a  parameter  of  the  off  diagonal  terms  in  the  beam  equations.  The  eighth  lay-up  is  an  iso¬ 
tropic  configuration  and  is  used  to  verify  the  accuracy  of  the  blade  model.  These  ply 
configtuations  can  be  seen  in  Table  5.3. 


LAMINA 

SYMMETRIC 
nON  CONFIGURATIONS 

PLY 

nf  A  r\c  1 

OT  A  r\c  ^ 

*  rvT"  J 

nuixuc  ** 

DLJ<Uti  b 

3LADE  D 

tyiMjt,! 

BLAues 

1 

0“ 

0“ 

0“ 

0“ 

■IIIIQIHI 

0“ 

0“ 

2 

30“ 

i^jn^niiiii 

30“ 

30“ 

30* 

45“ 

3 

0" 

15“ 

30“ 

45“ 

60“ 

75* 

90“ 

135“ 

90“ 

90“ 

90“ 

90* 

90“ 

90“ 

90“ 

90“ 

MIDPLANE 

s 

iHSiSUH 

HKiSMH 

HKtiifll 

IDB 

HHQSH 

30“ 

45“ 

75“ 

90“ 

135“ 

mm 

30“ 

30“ 

30“ 

IHESiHi 

30“ 

HjjHlSH 

45“ 

Si 

0“ 

0" 

0“ 

0" 

0 

0 

0“ 

0“ 

TABLE  5.3 


Table  5.4  shows  the  laminate  properties  for  each  blade.  They  were  calculated  using  the 
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laminate  analysis  program  in  Appendix  A. 


LAMINATE  PRt 

IPERTIES 

B 

BLADE  1 

BLADE 2 

BLADE3 

BLADE4 

BLADE5 

BLADE6 

BLADE7 

BLADES 

1. 

11.63x10* 

6.11  X  10* 

1.43  X  10* 
.117 
.062 

9.19*  10* 
6.11*  10* 
1.38*  10* 
.105 
.070 

7.63*  10* 
6.13*  10* 
1.63*10* 
.123 
.099 

7.23  *  10* 
6.44*10* 
1.85  *  10* 
.106 
.094 

7.29*  10* 
7.29*  10* 
1.71*10* 
.057 
.057 

7.33  *  10* 
8.96*  10* 
1.44*  10* 
.023 
.028 

7.24*  10* 
10.51*  10* 
1.36*  10* 
.068 
.099 

7.58*  10* 
7.58*  10* 
2.92*  10* 
.298 
.298 

TABLE  5.4 


5.4.  Blades  I*Vn  (Anisotropic) 

The  first  seven  blades  were  designed  to  exhibit  bending-coupling  characteristics. 
Using  the  LAMS  option  for  laminated  materials,  an  eight  layer  composite  was  built. 
Each  layer  had  the  same  thickness  of  0.015  inches.  The  four  layers  on  the  top  were  sym¬ 
metric  with  the  bottom  four  plies.  The  angles  of  orientation  can  be  seen  in  Table  5.3 
above.  The  60140  mix  of  T50  Fiber  (Graphite)  and  Fiberite  934  Epoxy  Resin  was  used. 
The  angle  of  ply  orientation,  0,  was  changed  systematically  for  the  third  and  sixth  layers. 
The  first  blade  has  0  =  0*’  and  blade  seven  has  0  =  90®  with  a  15°  increment  for  blades 
two  through  six.  Four  forces,  M^,  My,  and  P,  were  applied  individually  to  each 
blade. 


5.5,  Blade  Yin  (Isotropic) 


'ITiis  blade  model  was  designed  to  exhibit  isotropic  characteristics.  It  was  built 
using  the  LAMS  option  of  PATRAN®.  There  are  eight  layers  of  equal  thickness  (0.015 
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inches).  The  angles  of  ply  orientation  are  shown  in  Table  5.3  and  should  produce  isotro¬ 
pic  behavior  when  the  forces  are  applied.  Once  again,  four  forces,  M^,  My,  M^  and  t, 
were  applied  individually. 
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Chapter  Six 


Results 

Numerical  calculations  are  carried  oai  through  finite  element  analysis  for  a  compo¬ 
site  material  rotor  blade.  The  blade  is  discretized  into  616  elements  and  the  data  is  col¬ 
lected  ac  twelve  nodes;  six  along  the  leading  edge  and  six  along  the  trailing  edge.  The 
behavior  at  two  blade  sections  is  evaluated.  One  section  (Section  1)  is  located  at  the 
center  of  the  blade  length  and  the  other  section  (Section  2)  is  located  in  the  last  quarter  of 
the  blade  but  at  a  sufficient  distance  from  the  end  cap.  The  model  is  analyzed  for  eight 
blade  laminations  (seven  anisotropic  and  one  isotropic)  which  undergo  four  seperate 
loading  conditions. 

The  output  from  ANSYS*  (See  Appendix  E),  consists  of  displacement’s  in  the  three 
principal  directions;  UX,  UY,  and  UZ  for  each  node.  The  selected  nodal  displacements 
are  entered  into  a  data  file  for  further  evaluation  by  the  FORTRAN  code  in  Appendix  C. 
This  program  calculates  the  rotations  about  each  axis  (P,  <j),  y  for  x,y,z  respectively)  and 
the  displacements  along  the  longitudinal  axis.  The  data  is  also  manipr.lated  to  produce 
the  first  derivative  of  each  of  these  values  with  respect  to  the  blade  span.  The  stiffness 
matrix,  [C  ]  ,  is  formed  for  both  sections  of  each  blade.  A  modified  finite  differencing 
method  is  used  to  calculate  the  rotations  from  the  given  displacements  and  is  explained 


below. 
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The  beam  equadoii  matrix  which  is  represented  by  Equation  (4.26)  can  be  expressed 


agam  as: 


an  ai2  ai3  au 

021  (122  ^23  ^24 

^31  <*32  <*33  <*34  | 

^a4i  a42  a43  044  P 


Since  the  forces  are  known,  one  must  calculate  and  in  order  to  solve 

Bx  dx  dx  dx 

for  the  stiffness  terras  in  the  matrix  above.  The  use  of 


93  _  Pi-t-i-Pi-i 
dx  ~  21 

where 


Pi+i  - 


-[uzhi-uzLi] 


and  similarly 


-\UZ 


P.-i  = 


ti  -  uzu 

¥ 


and  where  i-1,  i,  and  /+1  arc  consecutive  increments  of  length  (/)  along  the  blade  length 
allows  for  the  computation  of  p  and  its  first  derivative  and  leads  to  the  solution  of  the 
first  column  of  the  stiffness  matrix. 


If 
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UZU 

I 


(6.5) 


where 


t;zf  =  (/2f'[i-c] 

and 

uzu^  -=uz\^x  [i-c]  +t^J’+iC 

then 


a4>  -2UZf-UZU 

dx 

And  similarly, 


(6.6) 


(6.7) 


(6.8) 


^  UYUi-2UYf-UYU 
dx  ~  /2 


Finally,  the  derivative  of  u  is  defined  as: 


(6.9) 


Bu  _  UXUi-UXU 

dx  21 


(6.10) 


Equations  6.8,  6.9  and  6.10  are  used  to  solve  for  the  second,  third  and  founh  columns  of 
the  stiffness  matrix  respectively.  The  matrices  are  then  assembled  for  each  blade.  They 


are  shown  below. 


MATRIX  FOR  SECTION:  1,  BLADE:  1 

.582198E^  -,648599E-08  .110397E-09  -.372616E-08 
-.646000E-08  .259820E-07  -.59083  lE-09  -.966120E-08 
.202181E-09  -.59'i695E-09  .579233E-08  .210487E.09 
-.357498E-G8  -.969?97E-08  .195244E-09  .501077E-07 

MATRIX  FOR  SECTION:  2,  BLADE:  1 

.576136E-07  -.6776^  1.E-08  .129622E-09  -.360332E-08 
-.63654'/E-08  .258174E-07  -.68:..216E-09  -.968831E-08 
-.138691E-09  -.536155E-09  .565750E-08  .274261E-10 
-.4I2470E-08  -.972916E-08  -.237509E-09  .495678E-07 


MATRIX  FOR  SECTION:  1,  BLADE:  2 

.597071E-07  -.998768E-08  .828957E-10  -.101376E-07 
-.996914E-08  .2870I7E-07  -.391494E-09  -.691174E-08 
.210321E-09  -.396723E-G>  .669389E-08  .7958122-10 
-.993I47E-08  -.693385E-08  .585059E-10  .555784E-07 

MATRIX  FOR  SECTION:  2,  BLADE:  2 

.592261E-07  -.104780E-07  .649079E-10  -.100703E-07 
-.981416E-08  .283538E-07  -.585838E-09  -.694357E-08 
-.130799E-09  -.419017E-09  .641803EO8  -.478981E-10 
-.104133E-07  -.722706E-08  -.769398E-09  .550485E-07 


MATRIX  FOR  SECTION:  1.  BLADE:  3 

.500480E-07  -.111859E-07  .221594E-10  -.971455E-08 
-.111937E-07  .308369E-07  -.194795E-09  -.400552E-08 
.138787E-09  -.190835E-09  .737674E-08  .145800E-11 
-.952498E-08  -.399955E-08  .483645E-li  .604304E-07 

MATRIX  FOR  SECTION;  2,  BLADE;  3 

.497724E-07  -.118054E-07  -.164852E-10  -.969432E-08 
-.110411E-07  ,305075E-a7  -.370755E-09  -.405130E-08 
-.306253E-09  -.238398E-09  .703087E-08  -.109083E-09 
-.100369E-07  -.433010E-08  -.902971E-09  .598740E-07 


MATRIX  FOR  SECTION:  1,  BLADE:  4 

.438927E-07  -.114499E-07  -.2959I7E-li  -.422881E-08 
..114735E-07  .307596E-07  -.104341E  09  -.212956E-08 
.902146E-10  -.995612E-10  .742222E-08  -.171720E-10 
..408496E-08  -.212293E-08  -.112840E-10  .604906E-07 
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MATRIX  FOR  SECTION:  2.  BLADE:  4 

.437602E-07  -.121026E-07  .872368E-11  -.419622E-08 
-.113331E-07  .304896E-07  -.217743E-09  -.216331E-08 
-.404176E-09  -.585793E-10  .704398E-08  -.153504E-09 
-.465988E-08  -.225900E-08  -.887374E-09  .599169E-07 


MATRIX  FOR  SECTION:  1.  BLADE:  5 

.473131E-07  -.118224E-07  -.357026E-10  -.135087E-10 
-.118509E-07  .290284E-07  -.993379E-12  .184891E-11 
.404266E-10  -.140037E-12  .701748E-08  -.167640E-10 
.927771E-10  -.I32395E-11  -.189042E-10  .565798E-07 

MATRIX  FOR  SECTION:  2.  BLADE:  5 

.472095E-07  -.124535E-07  .223954E-10  .I39275E-10 
-.116972E-07  .287571E-07  -.564925E-10  -.283551E-10 
-.365083E-09  .118199E-09  .662338E-08  -.184275E-09 
-.466199E-09  .330131E-10  -.905433E-09  .560019E-07 


MATRIX  FOR  SECTION:  1,  BLADE:  6 

.566060E-a7  ..112748E-07  -.112823E-09  -.139777E-08 
-.112952E-07  .264781E-07  .201259E-09  .31521  lE-08 
-.452902E-10  .200743E-09  .633755E-08  .155130E-1C 
-.129999E-08  .314820E-08  ,374775E-11  .510225E-07 

MATRIX  FOR  SECTION:  2,  BLADE:  6 

.119927E*05  -,117927E-07  -.825691E-10  -.145466E-08 
.410326E-05  .261994E-07  .182398E-09  .309478E-08 
-.257063E-09  .250628E-09  .600472E-08  -.165144E-09 
-.170257E-08  .312056E-08  -.807286E-09  .504383E-07 


MATRIX  FOR  SECTION:  1,  BLADE:  7 

.607422E-07  -.802202E-08  -,136835E-09  -.840365E-08 
-.803070E-08  .247008E-07  .377322E-09  .546028E-08 
-.693500E-10  .37772 lE-09  .578679E-08  .74385 IE- 10 
-.827150E-08  .546064E-08  .•505670E-’0  .476835E-07 

MATRIX  FOR  SECTION:  2.  BLADE;  7 

.603513E-07  -.834859E-08  -.168020E-05  -.855583E-08 
-.784810E'08  .245449E-07  ,337192E-09  .53824ih-08 
-.123595E-09  .262l84E-)9  .564403E-08  -.122913E-09 
-.848671E-08  .517938E-08  -.345758E-09  .470697E-07 
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MATRIX  FOR  SECTION:  1.  BLADE:  8 

.287997E-07  -.150000E-11  .563859E-n  -.935840E-11 
-.259201E-11  .259144E-07  -.178833E-il  -.248651E-12 
.184126E-10  ~.232749E-12  .619388E-08  -.486900E-10 
.453763E-10  .968932E-12  -.228322E-10  .543169E-07 

MATRIX  FOR  SECTION:  2,  BLADE:  8 

.285857E  07  .166050E-09  .185933E-10  -.479692E-11 
.275293E-09  .259307E-07  .158014E-10  -.262454E.il 
-.276093E-09  .804174E-12  .632200E-08  -.187553E-09 
-.347824E-09  .253299E-10  .373140E-09  .538126E-07 

To  examine  the  accuracy  of  the  resultant  matrices,  it  is  important  to  begin  with  the 
isotropic  blade.  Examination  of  Blade  8  reveals  that  the  diagonal  tenns  are  all  positive 
and  representative  of  the  correct  responses.  Theoretically,  the  off-diagonal  terms  of  this 
matrix  should  be  zero.  Since  the  finite  clement  method  was  used  for  analysis  and  numer¬ 
ical  methods  are  not  exact,  these  terms  will  never  quite  be  zero.  The  values  of  these 
terms  range  from  two  to  five  orders  of  magnitude  smaller  than  the  diagonal  terms  and 
can  be  considered  zero.  The  only  term  which  is  marginally  acceptable  is  the  043  term 
(Equation  6.1).  Tlic  value  is  only  one  or  two  orders  of  magnitude  smaller  than  the  diago¬ 
nal  term  in  that  column  and  may  represent  some  unwanted  coupling  in  the  isotropic 
model.  This  term  represents  an  extensional-bending  coupling  caused  by  a  moment  about 
the  z-axis.  This  is  probably  due  to  the  fact  tiiat  the  bending  moment  is  applied  about  tlie 
usymmetric  axis  of  the  model  and  that  the  point  of  application  is  in  eitur.  Because  of  this 
discrepancy  in  the  isotropic  case,  the  a  43  term  in  each  of  the  other  blades  must  be  con¬ 
sidered  to  be  slightly  inaccurate.  There  are  distinct  coupling  terms  present  in  each  of  the 
anisotropic  blades.  The  value  of  a  coupling  term,  say  a  13,  of  one  lamination  can  be  com¬ 
pared  to  the  same  temi  of  another  anisotropic  lamination  to  determine  the  impact  that  the 
ply  orientation  angle  has  on  the  coupling  effect.  This  information  could  be  stored  in  a 
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file  which  contained  similar  data  for  numerous  blades  and  laminations  and  might  be  use¬ 
ful  for  an  aeroclasdcian  who  is  attempting  to  construct  a  composite  blade  which  exhibits 
elasdc  coupling  responses  desirable  for  given  flight  conditions.  It  might  even  be  possible 
with  enough  research  and  data  collection  to  construct  a  catalog  of  composite  blades 
which  exhibit  unique  aeroelastic  properties. 

After  all  of  the  calculations  are  completed,  the  response  of  P',  <j>\  and  u'  are 
plotted  with  respect  to  the  ply  angle  0.  The  results  are  shown  in  Figures  6.1  through  6.8. 

In  the  plot  of  p'  versus  the  ply  angle,  (Figure  6.1),  the  results  are  as  expected.  The 
rate  of  twist  of  the  blade  (for  a  fixed  torque)  should  be  inversely  proportional  to  the  in¬ 
plane  shear  modulus,  Pxy»  of  the  laminate.  The  value  of  the  laminate  in-plane  shear 
modulus  increases  from  0  =  0®  to  0  =  45°  where  it  reaches  a  maximum.  As  the  ply  angle 
continues  to  increase  to  0  =  90°,  the  in-planc  shear  modulus  decreases.  If  the  values  of 
for  this  model  (See  Table  5.4)  ate  plotted  with  respect  to  the  ply  angle,  the  inverse  of 
Figure  6. 1  is  the  result.  Both  of  these  plots  are  offset  to  the  right  so  that  the  minimum 
angle  of  twist  and  the  naaximum  in-plane  shear  stress  arc  located  near  49°,  This  is  due  to 
the  geometry  of  the  blade  which  employs  a  curved  laminate  and  therefore  the  effect  of 
Gjcy  is  altered. 

Figiues  6.2,  6.3  and  6.4  demonstrate  the  effect  of  the  ply  angle  on  the  curvature, 
rotation  and  rate  of  displacement  of  the  blade.  These  plots  should  follow  the  functional 
dependence  of  the  laminate  Ex  and  Ey  on  the  ply  angle.  The  results  correlate  closely 
with  the  moduli  from  0  =  0°  to  0  =  45°.  However,  one  would  expect  the  curvature  and 
rotation  to  continue  to  increase  as  the  the  ply  angle  increases  up  to  0  =  90°.  This  is  the 
expected  behavior  since  the  strength  in  the  longitudinal  direction  is  about  /7  times 
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greater  than  the  strrrgth  in  the  transverse  direction.  The  UZ  and  UY  displacements  were 
examined  for  each  blade  and  were  determined  Pj  he  correct.  The  laminate  properties 
(Table  5.4)  do  not  behave  as  expected  but  have  been  verified  by  three  laminate  analysis 
programs.  It  is  also  assumed  that  the  geometry  of  the  blade,  modeled  as  a  hollow  unsym- 
metric  laminated  beam,  causes  the  deviation  from  the  plot  of  £jc  and  Ey  which  actually 
correspond  to  a  fiat  laminated  plate. 

The  responses  of  P,  4>,  y  and  u  arc  plotted  along  the  blade  length  and  are  shown  by 
Figures  6.5  through  6.8.  The  results  are  as  expected. 


In  Figure  6.5,  which  is  the  plot  of  the  angle  of  twist  along  the  blade  length,  it  is 
expected  that  the  twist  of  Blade  4  (See  Chapter  6  for  blade  designation  and  description) 
would  exhibit  the  minimum  values  for  the  anisotropic  cases.  It  is  also  expected  that 
Blade  1  and  Blade  7  would  be  similar  to  each  other  and  reach  the  highest  values  fr»r  twist. 
This  was  tiic  case  in  Figure  6.5.  The  isotropic  blade  (Blade  8)  results  are  plotted  on  this 
graph  as  a  reference  or  base  line  case  and  are  depicted  by  a  straight  line  which  reflects  a 
minimum  value  linear  twist  Figures  6.5  through  6.8  also  accurately  refiect  the  increase 


in  twist,  CuTratuie,  rctaticti  and  displawcuiciU  along  the  blade  span  as  the  distance  from 


the  constrained  end  increases.  There  do  not  appear  to  be  any  negative  end  effects  and 
one  can  assume  that  the  data  was  taken  at  a  sufficient  distance  from  the  end  of  the  beam. 


The  nodal  displacements  and  blade  deformations  provide  valuable  information  in 
determining  the  elastic  coupling  responses  of  the  anisotropic  blade  muJels.  They  also 
provide  an  acceptable  vehicle  for  comparing  the  effects  of  ply  orientation  as  a  design 
parameter.  In  reference  to  the  blade  displacements,  two  points  need  to  be  made.  First, 
the  end  caps  were  designed  to  be  extremely  stiff.  To  veiify  the  accuracy  of  the  force 
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application  on  the  end  cap,  two  procedures  were  attempted.  A  bending  moment  of  1000 
units  was  applied  at  one  node  on  the  end  cap  and  the  displacements  were  recorded.  The 
bending  moment  was  then  distributed  about  each  node  on  the  outside  edge  of  the  endcap 
and  the  displacements  were  again  recorded.  The  displacements  were  identical  in  both 
cases  and  verified  that  the  end  cap  was  not  deforming  undesirably  and  that  the  moment 
application  was  correct.  The  second  point  of  interest  is  that  the  outside  surface  of  tlie 
blade  had  a  tendency  to  deform  unfavorably  under  stress  and  as  a  result  the  displacement 
data  was  only  taken  at  the  leading  and  trailing  edges.  To  correct  for  this,  the  model  needs 
to  be  reinforced  with  ribs  along  the  entire  length  of  the  blade.  These  ribs  could  be  con¬ 
structed  in  a  similar  fashion  as  the  end  caps. 

Figure  6.9  represents  a  theoretical  view  of  what  might  happen  when  a  single  force  is 
applied  to  an  anisotropic  blade.  The  actual  displacements  can  be  seen  in  Figures  6.10 
through  6.41  which  are  located  at  the  end  of  this  chapter.  As  a  point  of  reference,  the 
isotropic  blade  (Blade  8)  needs  to  be  examined.  Figures  6.38  through  6.41  present  a 
visual  representation  of  the  blade  response  of  an  isotropic  lamination  sequence  under  tor¬ 
sion,  bending  and  extension.  Figure  6.38  reveals  that  as  a  constant  torque  is  applied  at 
the  shear  center,  the  entire  blade  twists  around  that  point.  When  comparing  that  with  an 
anisotropic  lay-up.  Figure  6.18  for  example,  one  can  see  that  the  blade  twists  about  a  dif¬ 
ferent  axis  which  is  located  toward  the  trailing  edge  and  that  there  is  some  flapping  or 
bending  about  the  y-ar.is.  This  is  the  coupling  which  is  a  result  of  the  anisotropic  lamina¬ 
tion  and  which  is  predicted  iu  the  bending-stiffness  matrix,  [D],  of  lamination  theory.  A 
similar  investigation  of  the  application  of  an  axial  force  along  the  neutral  axis,  (Figure 
6.41),  reveals  thai  as  the  isotropic  blade  is  pulled,  the  result  is  extension  in  the  x  direction 
and  nothing  else.  When  this  is  compared  to  the  response  "  an  anisotropic  blade  under 
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identical  loading  conditions,  Figure  6.21  for  example,  one  can  see  that  in  addition  to 


P 


FIGURE  6.9.  Elastic  Bending-Torsion  Coupling  of  Rotor  Blade 

extension,  the  blade  underwent  bending  about  the  y  and  z  axes  and  twisting  about  the  x- 
axis.  This  can  be  predicted  by  the  extensional-stiffness  matrix,  [A],  and  the  bending- 
stiffness  matrix  [D].  Since  the  laminates  used  in  this  research  were  symmetric  (Sec  page 
2),  the  cxtensional-bending  coupling  matrix,  [B],  is  zero  and  would  not  nonoaDy  have 
any  effect  on  the  response  of  the  model.  However,  since  the  model  is  not  symmetric  with 
respect  to  tlic  x-z  plane,  there  are  some  minor  effects  from  this  type  of  coupling.  Tliese 


can  be  seen  in  Figures  6.19  and  6.20.  All  of  the  blade  deformations  are  depicted  in  Fig¬ 
ures  6.10  through  6.41  and  can  be  used  to  examine  the  varying  effects  of  lamination 
sequence  on  displacement  and  coupling. 

The  results  provided  by  each  of  these  media  allow  us  to  examine  the  effects  of  ply 
orientation  and  to  compute  the  coupling  terms  of  the  accompanying  stiffness  matrices. 
This  work  should  lead  to  a  much  improved  capability  for  evaluating  blade  properties  for 
use  in  rotorsystems. 
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FIGURE  6.10 


Blade  One  (Anisotropic)  Under  Vorsionai  Load 
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FIGURE  6.12 


Blade  One  (Anisotropic)  Under  Bending  Moment  (Z) 
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FIGURE  6.13 

Blade  One  (Anisotropic)  Under  Axial  lx)ad 
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FIGURE  6.14 

Blade  Two  (Anisotropic)  Under  Torsional  Load 
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FIGURE  6.15 

Blade  Two  (Anisotropic)  Under  Bending  Moment  (Y) 
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FIGURE  6.16 

Blade  Two  (Anisotropic)  Under  Bending  Moment  (Z) 
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FIGURE  6.17 

Blade  Two  (Anisotropic)  Under  Axial  Load 
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FIGURE  6.  i  8 


Blade  Three  (Anisotropic)  Under  Torsional  Load 
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FIGURE  6.19 

Blade  Three  (Anisotropic)  Under  Bending  Moment  (Y) 


FIGURE  6.20 

Blade  Three  (Anisotropic)  Under  Bending  Moment  (Z) 
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FIGURE  6.21 

Blade  1  nn  e  (Anisotropic)  Under  Axial  Load 
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FIGURE  6.22 


Blade  Four  (Anisotropic)  Under  Torsional  Load 
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FIGURE  6.23 


Blade  Four  (Anisotropic)  Under  Bending  Moment  (Y) 
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FIGURE  6.24 

Blad(-  Four  (Anisotropic)  Under  Bending  Moment  (Z) 
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FIGURE  6.25 

Blade  Four  (Anisotiopic)  Under  Axial  Load 
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HGURE  6.26 

Bladt  Five  (Anisotropic)  Under  Torsional  Load 
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FIGURE  6.27 

Blade  Five  (Anisotropic)  Under  Bending  Moment  (Y) 
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FIGURE  6.28 

Blade  Five  (Anisotropic)  Under  Bending  Moment  (Z) 


ANSY3  4 . 3 
MAR  E8  1989 
10:55:56 
P0ST1  DISPL. 
DMX  =0.02256 


BLRDE5P  I  1  YF  =4.5 


X  CJ  N  X  C 
■ni-HC  <ICM 

:5  s: 

i>-<xox£: 

-C  X  ON 

O  M  CO  t)  3 

M 

CM 

m 

-n  *0  M  c 

2 

H  H  CO  IT)  3>  3^ 

CO  3 

:z 

O  CO  3 

CO 

n:  01  CO  ••  ^  tz 

H  o 

H 

o 

r~  HO 

H 

:o  "0  CO  CO 

II  II  11  II  II  II 

II  II  II 

II  II 

II  II  II  l)  II  II 

II  II  II  II 

II 

II  II  -V  nj  X 

(Jj  1  ->■ 

-Ck 

(jj 

1  H  js.  ro 

4>. 

CO 

— »  — ^  •’  H  CO 

m 

•  CS3  4>> 

CO  •  o  . 

•  O 

. 

O  Jk 

.813 

cn 

csjcn  CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

00  ^ 

CO  CO  4i- 

TO  00  • 

r~  CO  CjJ 

• 

FIGURE  6.29 


Blade  Five  (Anisotropic)  Under  Axial  Load 
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Blade  Six  (Anisotropic)  Under  Torsional  Load 
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Blade  Six  (Anisotropic)  Under  Bending  Moment  (Y) 
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FIGURE  6.32 

Blade  Six  (Anisotropic)  Under  Bending  Moment  (Z) 


FIGURE  6.33 

Blade  Six  (Anisotropic)  Under  Axial  Load 


A 


XF  =40 

BLADE7MX  YF  =4.5 


97 


r 


ON-<;xc 
M  C  C  C  tH 

c 

l-l 

n  ^  i-H  c 

o  M  cn 
3  H  H 

cn  3 

on 

zz. 

Ch  in 

:z 

on 

X  m  m 

H  o 

H 

o 

r-  H 

o 

— 1 

3C  "D 

II  II  II  II  II 

II  II  II 

II  II 

II  II  II  II  II 

II 

II  II  II  II 

II  II  II 

00-*  1 

Jx 

-k  00 

1  41  — 

ro 

41  41  41  -<- 

03  -1  -* 

•  03  4^ 

■  CD  • 

•  CO  41 

" 

on 

CD  on  (o 

on 

d 

on 

4i 

00  CO 


FIGURE  6.34 


Blade  Seven  (Anisotropic)  Under  Torsional  Load 
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Blade  Seven  (Anisotropic)  Under  Bending  Moment  (Y) 
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Blade  Seven  (Anisotropic)  Under  Bending  Moment  (Z) 
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Blade  Seven  (Anisotropic)  Undi;r  Axial  Load 
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Under  Torsional  Load 
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Blade  Eight  (Isotropic)  Under  Bending  Moment  (Y) 
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Blade  Eight  (Isotropic)  Under  Bending  Moment  (Z) 
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Blade  Eight  (Isotropic)  Under  Axial  Load 
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Chapter  Seven 


Summary  and  Concluding  Remarks 


The  advent  of  composite  material  construction  in  helicopter  rotor  blade  design  has 
opened  new  areas  of  research  and  development.  In  particular,  the  unique  coupling  pro¬ 
perties  afforded  by  anisotropic  composites  provide  a  potentially  powerful  new  set  of 
design  variables  for  designing  aeroelastic  stability.  The  development  of  shell  elements 
with  anisotropic  laminate  capabilities  which  are  incorporated  into  commercially  avail¬ 
able  finite  element  packages  makes  this  type  of  research  possible. 

Composites  have  better  fatigue  characteristics  than  metal  and  they  permit  the  flexi¬ 
bility  of  tailoring  the  structural  properties.  One  drawback  of  composite  materials  is  that 
their  structural  characteristics,  in  particular  the  couplings  of  different  bending  modes,  are 
not  very  well  understood  by  a  typical  design  engineer.  Due  to  the  complex  directional 
nature  of  composite  materials,  the  analysis  of  anisotropically  laminated  structures  is 
complicated  by  the  introduction  of  these  elastic  couplings.  This  is  especially  true  in 
dynamic  behavior  since  the  elastic  couplings  have  such  a  significant  influence  on  charac¬ 
teristics  like  mode  shape  and  frequency.  Recent  design  practice  has  been  to  treat  compo¬ 
sites  similarly  to  metals.  This  approach,  however,  does  not  permit  the  description  of  gen¬ 
eral  composite  lay-ups  and  its  accuracy  is  questionable  in  view  of  the  following  facts 
[19]: 
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(a)  The  use  of  laminated  composite  materials  for  the  blade  structure  results  in 
signincant  shearing  and  warping  deformations.  These  effects  are  far  more 
pronounced  for  highly  anisotropic  composites  than  for  metallic  materials. 

(b)  The  proper  tailoring  of  composite  lay-ups  will  result  in  elastic  couplings,  e.g., 
bending-twisting  or  extension-twisting  couplings  which  are  known  to  strongly 
influence  the  dynamic  behavior  and  stability  of  die  blade.  Such  couplings  are 
not  appropriately  modeled  with  standard  beam  theories. 

(c)  in  specific  applications,  the  blade  planform  becomes  wider,  resulting  again  in 
increased  shearing  ar>d  warping  deformations. 

(d)  Recent  research  shows  that  the  torsional  behavior,  and  the  effects  of  pre-twist 
of  the  blade  are  accurately  predicted  only  if  out-of-plane  warping  of  the  cross- 
section  is  included. 

In  this  research,  the  different  stiffness  coupling  terms  for  a  composite  rotor  blade 
were  identified.  These  included  the  terms  caused  by  bending-torsion  and  extension- 
torsion  couplings  of  angle  plies.  It  was  shown  that  the  ply  lay-up  and  orientation  has  a 
substantial  effect  on  the  blade  aeroelasdc  characteristics.  It  is  concluded  that  the  ply 
orientation  angle  and  the  ply  thickness  are  suitable  parametric  variables  for  future  studies 
aimed  at  assessing  the  influence  of  variations  in  structural  properties  which  can  in  turn  be 
used  in  the  prediction  of  natural  frequencies  and  aeroelasdc  stability. 

The  complex  structural  behavior  of  a  rotor  blade  was  modeled  accurately.  The  fol¬ 
lowing  features  were  included:  thin-walled  cross-sections,  proper  scaling  and  material 
anisotropy.  This  resulted  in  an  accurate  modeling  of  the  torsion,  and  coupled  extension- 
twist  behavior.  While  the  undeformed  blade  geometry  was  accurate,  the  deformed 
geometry  was  somewhat  deficient.  The  deformation  was  skewed  at  locations  on  the 
blade  surface.  This  can  be  corrected  by  adding  stiff  ribs  (similar  to  the  end  caps) 
throughout  the  length  of  the  structure.  Despite  the  fact  that  the  geometries  presented  in 
this  paper  were  rather  simple,  the  formulation  is  applicable  to  complicated  geometries. 


Although  we  cannot  increase  the  accuracy  of  prcdicred  trequencies  for  the  bending 
modes,  which  are  basically  uncoupled,  we  can  greatly  improve  the  prediction  of  the  fre¬ 
quency  of  the  torsional  mode.  This  mode  is  elastically  coupled  and  since  there  is  an 
absence  of  off-diagonal  coupling  terms  in  the  formulation  of  the  equivalent  engineering 
stiffnesses  using  beam  theory,  they  must  be  predicted  using  lamination  theory  or  deter¬ 
mined  experimentally. 

In  this  research,  the  extension-twist  terms  were  most  pronounced.  Potential  applica¬ 
tions  of  extension-twist  coupling  include  the  design  of  a  tilt-rotor  blade  that  would 
change  twist  as  a  function  of  rotor  speed.  Such  a  blade  would  be  extremely  useful  on  an 
aircraft  such  as  tlie  V-22  OSPREY.  Besides  the  potentially  positive  aeroelastic  benefits  of 
twist,  there  are  other  useful  applications  of  elastic  coupling.  These  include  the  redurtion 
or  modification  of  structural  vibration  in  rotor  blade  design.  This  is  accomplished  by 
building,  into  the  blades,  a  natural  aeroelastic  compliance  for  vibration  reduction  [20]. 

To  apply  this  model  toward  further  analysis  of  elastic  coupling,  one  would  need  a 
three-dimensional,  non-linear  beam  element  which  includes  shear  and  warping  effects  as 
well  as  a  correct  modeling  of  elastic  couplings.  This  type  of  element  is  unavailable  in 
ANSYS*  and  one  would  need  to  use  a  program  such  as  MSC/NASTRAN*. 
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Appendix  A 


Laminate  Analysis  Program 


PROGRAM  COMMAT  (WUTGinPL]T.TAPE5=INPm,TAPE6=OUTPUT) 

itr*  %  4ii|i  4i4i  4i  i(i  i|i  tliif  fi  <*i « <4*  4(  4*  >t(  4i  4>  4*  i|(  4<  >l>  4<  %  %  if>  %  4i 

coMPOsrrE  materials 

LAMINATE  ANALYSIS  PROGRAM 

PATRICK  GRAHAM  FORRESTER 
UNIVERSITY  OF  VIRGINIA  1989 

THESIS  RESEARCH 

I|i4(4i4t4i  4c  411(1 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  CALCUl  ATE 
THE  LAMINATE  PROPERTIES  OF  COMPOSITE  MATERIAL 
HELICOPTER  BLADE  MODELS  GIVEN  THE  LAMINA 
MATERIAL  PROPERTIES  OF  GRAPHITE/EPOXY 

I»c^i^i4i44i4i4i*4i*4c4i4i4>4i4c4<4<4i4>4i4>4'4>4>4<4i4<c^4i4'4i4i4i4i4i4i4i«4<*c»4ic(i4i4>4i4> 


INTEGER  lAM  JSf.IAINV.EER 

INTEGER  NOMAT.ANGLE.PLIESd.OAD.ILOAD.ITEMP.QUEST 
REAL  SI  1(9),S12(9),S13(9),S22(9).S23(9),S33(9).S44(9),S55(9) 

REAL  El  1(9)E22(9),E33(9),G12(9).G23(9).G31(9).S66{9) 

REAL  NU 1 2(9)  JW23(9)2W  1 3(9).G  1 3(9),RTHCTA(9) 

DIMENSION  D3(3,3).EPS3(3,9) 

INTEGER  THETA(9).ISYM 

REAL  Q1 1(9),Q12(9),Q22(9),Q66(9),QBAR1 1(9,9),QBAR12(9,9) 

REAI.  QBAR22(9.9),QBAR16(9.9),QBAR26(9.9),QBAR66(9,9) 
REALEXX(9,9)£YY(9.9)£XY(9,9)2WXY(9.9),GXY(9,9) 

REAL  ETAXYX(9,9)ETAXYY(9.9),KAP(3),H 

REAL  NX JTir.NX Y AIX  J4Y>IX Y,DELTT,TEMP.TEMPN,TEMPM, ALAM(3) 
DIMENSION  EPS  1T(9,3).EPS  1  B(9,3),XST(9),YST(9),SSS(9) 

DIMENSION  RMB(3),RNB(3),EPSB(3).EPSXB(9,3) 

DIMENSI0NALP1(9),ALP2(9),ALP3(9).ALPX(9.9).ALPY(9,9).ALPXY(9,9) 

DIMENSIONSIGXB(9,3),SIGXT(9,3),EPSXT(9,3),EPST(9,3),U(3,3) 

REAL  KAPT(3),KAPB(3) 


non  non  nnn  ooo 
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DIMENSION  EPSTH(3)X>EL(9,3)  jraC(9.3) 

REAL  F1,F2,F64"1 1 4^12  J^2,F66,BE  rAl(9,3).BETA2(9,3) 

DIMENSION  EPST1(9,3)3PSTH1(9,3) 

DIMENSION  RNT(3)jyviT(3).SIGlT;9,3).SIGlB(9, 3) 

REAL  MIDEPXMIDEPY>1IDE?XY,!LAPX.KAPY,KAPXY 

DW4ENSIONX(6.6).P(3.3).T(3.3)JR(6,6).PN(3),PM(3) 

REAI-PI,C.S.TOL,WKl(12).WK2(6).S1.(6).S2(3) 

REAL  ETABXYX,CTABXYY,ErABX:<YJETABYXY3BARXJEB  ARY 
DIMENSION  A(6,5.9).B(6,6.9)^(3.3.9.9).D(3.3,9.9) 

DIMF-NSION  MATN(9), ANGN(9),THICK(10),2(1 1).D4(3.3),A4(3.3).B4(3  3) 
DIMENSION  A1(3.3).B1(3.3).D1(3.3)^(3.3).B2(3,3).D2(3.3),B5(3,3) 
DIMENSION  EPS(3)  JIN(3).RM(3),P1  (3).P2(3),A3(3.3),B3(3.3) 

DIMENSION  SIG(3,9)^PS1(9, 3), SIGl(9.S).SIG2(3,9).SIG3(3,9j 
DIMENSION  SIGX(9,3).SIGY(3,9),SIGXY(3.9).EPS2(3, 9) 

DIMENSION  EPSX(9.3).EPSY(3,9)£PSXY(3, 9), AA(3,3)3B(3,3).DD(3, 3) 
DIMENSION  AT£MP(6.6.9).CTEMP(3.3.9.9).T1(3,3.9),T2(3,3,9) 

OPEN  (UNIT=6JFILE= ’RESULT’) 

. ASSIGN  INPUT  DATA . 

PI  =  3.14159 
STHICK  =  0.0 


. OBTAIN  LAMINATE  PARAMETERS . 

WRITE(5,*)  ENIER  THE  NUMBER  OF  DIFFERENT  MATERIALS  AND  ANGLES’ 
WRITE(5,*)  ’AND  THE  NUMBER  OF  PLIES  AND  LOADS  IN  LAMINATE’ 

READ (5*)  NOMAT^GLE^LIESJLOAD 

WRrrE(5 ,*)  ’IS  THIS  A  SYMETRIC  LAMINATE?  YES=0  NO=l  ’ 

READ(5,*)  ISYM 

WRITE(5 ,*)  ’IS  THERE  THERMAL  LOADING?  YES=0  NO=r 


WRITE(5,*)  ’DO  YOU  KNOW  NXX  FOR  FAILURE  CRITERIA?  YES=1  N0=0’ 
READ(5,’*)  QUEST 

. OBTAIN  VALUES  OF  ENGINEERING  PARAMETERS . 

DO  2  1=1, NOMAT 

WRITE(5,*)  ’ENTER  El  1,E22  FOR  MATERIAL  NUMBER  ’.I 
READ(5,*)  Ell(I),E22a) 

C 

WRnTCS,*)  ’ENTER  G12  FOR  MATERIAL  NUMBER 
READ(5,*)  G12(I) 

C 

WRITE(5  *)  ’ENTER  NU12,NU23  FOR  MATERIAL  NUMBER  ’,  I 
READ(5,*)  NU12(I),NU23(I) 

C 

WRITE(5,*)  ’ENTER  ALPHA  1,ALPHA2  AND  ALPHAS  FOR  MATERIAL  ’.I 
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READ(5  ,*)  ALP1(I).ALP2(D,ALP3(I) 

C 

WRITE(5  *)  ’ENTER  X  SUB  T.  Y  SUB  T,  S  AND  F12  FOR  MATERIAL  ’,I 
READ(5  ,*)  XST(I),YST(I),SSSa) 

C 

FI  =0.0 
F2  =  0.0 
F6  =  0.0 

Fll  =-1.0/(XSTa)'*'*2) 

F22  =  -1.0/(YSTa)**2) 

F66  =  1.0/(SSS(I)**2) 

E33a)  =  E22(I) 

G13(D  =  G12(I) 

G23(D  =  E22a)/(2.0*(L0  +  NU23a))) 

NU13(I)  =  NU12(D 

2  CONTINUE 

DO  3  I  =  1.ANGLE 

WRITEfS  *)  ENTER  ANGLE  NUMBER  I 
READ(5,*)THETA(I) 

RTHETA(I)  =  PI/180.0  ■*  THETA(I) 

3  CC'MTE-IUE 
D04I  =  1,PLIES 

V/KrrE(5.*)  ’ENTER  MATN.ANGN.THICK  FOR  PLY  NUMBER  ’  ,I 
;iEAD(5,*)  MATN(I).ANGN(I),THICKa) 

C 

STHICK  =  STHICK  +  THICK(I) 

4  CONTINUE 
C 

D05I  =  1U0AD 

WRnUCS,*)  ’DO  YOU  KNOW  [N]  AND  [M]?  YES=0  NO=r 
READ(5,*)  ILOAD 
IpgLOAD  .EQ.0)THEN 

WIGTECS,*)  ’ENTER  NX,  NY  AND  NXY  FOR  LOADING  CONDHION 
READ(5,’'')  PJ^(l),RN(2)jy^(3) 

WRrrE(5.*)  ’ENTER  IvIX,  mV AND  MXY  FOR  LOADING  CONOmON 
READ(5,*)  RM(1)JIM(2),RM(3) 

ELSE 

WRniE(5,*)  ’THAT  MEANS  YOU  KNOW  MIPLANE  S'fRAINS  AND  CURVATURES’ 
WRITEC5.*)  ’ENTER  EPSX,  EPSY  AND  EPSXY  FOR  LOADING  CONDITION  ’,I 
READ(5.*)  EPS(1),EPS(2),EPS(3) 

WRITE(5,*)  ENTER  KAPX,  KAPY  AND  KAPXY  FOR  LOADING  CONDITION 
READ{5.*)  K-VP(n,KAP(2),KAP(3) 

ENDIF 

5  CONTINUE 
C 

IF  OTEMP  .EQ.  1)  CO  TO  739 

WRrrE(5,*)  ’ENTER  CHANGE  IN  TEMPERATURE’ 

READ(5,*)  DELTT 
739  CONTINL’E 
C 
C 


uu  uuuuu  uuu 


. PERFORM  COMPUTATIONS . 

D0  6I=1JM0MAT 
SI  10)=  1.0/El  1(1) 

S12a)  =  -NUl2(I)/Elia) 

S22(I)  =  1.0/E22(I) 

S330)=1.O/E33(I) 

8440)  =  1.0/023(1) 

3130)  =  -NU130)/E110) 

S230)  =  -NU230)/E110) 

S550)  =  1.0/013(1) 

S66(I)  =  1.0/012(1) 

. CONSTRUCT  THE  COMPLIANCE  MATRIX, 

. SET  THE  INITIAL  VALUES  TO  ZERO . 


DO  1K=-.  1.6 
D01J=1,6 
A0C,J.I)  =  0.0 
1  CONTINUE 

. ESTABLISH  NON-ZERO  VALu'ES. 

A(L1,D  =  S110) 

A(1.2,.I)  =  S120) 

A(1,3.I)  =  S130) 

A(2.1,I)  =  S120) 

A(2.2.I)  =  S220) 

A(2.3,I)  =  S23(I) 

A(3,1,I)  =  S130) 

A(3,2,I)  =  S230) 

A(3,3.I)  =  S330) 

A(44,I)  =  S44(n 
A(5!5!i)  =  S550) 

At6,6,I)  =  S660) 

C 

D0  21  K=  1,6 
D0  21  J=  1.6 

21  ATEMP(K.J,r)  =  A(KJ.l) 

C 

DO  1000  K=  1.6 
DO  10001=1,6 
X0C,J)  =  0.0 
1000  CONTINUE 
X(1,1)  =  S11(I) 

X(1.2)  =  S12(I) 

X(1,3)  =  S13(I) 

X(2,1)  =  S12(I) 

X(2.2)  =  S22(I) 

X(2,3)  =  S23(I) 


115 


X(3,l)  =  S13(I) 

X(3^)  =  S23(I) 

X(3,3)  =  S33(D 
X(4.4)  =  S44(l) 

X(54)  =  S55(I) 

X(6.6)  *  S66(D 
C 

C  . CAI.CULATE  THE  STIFFNESS  MATRIX . 

C 

IA  =  6 
IAINV=:6 
M  =  6 
N=6 
TOL  =  0.0 
C 

CALL  LGINF  (X.IAJ^J4.T0L,R,IAINV,S1.WK1,IER) 

C 

XY1=R(1.1) 

XY2=:R(U) 

XY3  =  R(1.3) 

XY4=:R(2.1) 

XY5  =  R(2^) 

XY6  =  R(2.3) 

XY7  =  R(3,1) 

XY8aR(3;j) 

XY9  =  R(3.3) 

XY10  =  R(4.4) 

XYllaR(5^) 

XY12  =  R(6,6) 

XY13  =  0.0 
C 

DO  1001  K=  1,6 
DO  1001  J=  1,6 
B(KJ,I)  =  0.0 
1001  CONTINUE 
B(1,1,I)  =  XYl 
B(1.2,I)  =  XY2 
B(1,3,I)  =  XY3 
B(2,1,I)  =  XY4 
B(2,2,I)  =  XY5 
B(2,3,I)  =  XY6 
B(3,1,I)  =  XY7 
B(3.2,I)  =  XY8 
B(3,3,I)  =  XY9 
B(4,4,I)  =  XYIO 
B(5,5,I)  =  XY11 
B(6,6,I)  =  XY12 
C 

C  . CALCULATE  THE  REDUCED  STIFFNESS  MATRIX  WITH  ROTATION, 

C 

Qll(I)  =  S22(I)/((S11(I)*S22(I))  -  (S12(I)*S12(I))) 


Q12(I)  =  -S12(D/((S1  l(I)*S22a))  -  (S12(I)*SI2a))) 

Q22(D  «  Slia)/((SU(I)*S22(I))  -  (S12(I)*S12(I))) 

Q66(I)  =  L0/S66a) 

C 

C  . CONSTRUCTTHEREDUCEUSTIFFNESSMATRIXWITHROTATION, 

C 

D0  68K=  l^GLE 
C 
C 

C  =  COS(RTHETA(K)) 

S  =  SIN(RTHETA(K)) 

C 

QBARlia.K)  =  Q1  ia)*(C*C*C*Q  +  2.0*(Q12(I)  + 

A  2,0*Q66a))*(S*S)*(C*Q  +  Q22(I)*(S*S*S’*S) 

QB.\R12a.K)  =  (QlUO  +  Q22a)  -4.0*Q66(I))*(S'»S)*(C*C)  + 

A  Q12(I)*((S*S*S*S)+(C*C*C*Q) 

QBAR22a,K)  =  Qlia)*(S*3*S*S)  +  2.0*(Q12(D  +  2.0*Q66(I))*(S*S) 

A  *(C*C)  +  Q22a)*(C*C*'C*C) 

QBAR16a.K)  =  (Qll(I)  -  Q12a)  -  2.0*Q66(I))^S*(C*C*C)  +  (Q12(I)  - 
A  Q22(D  +  2.0*Q66a))  *(S*S*S)*C 

QBAR26a.K)  =  (Qll(I)  -  012®  -  2.0*Q66(I))*(S*S*S)*C  +  (Q12(I)  - 
A  Q22®  +  2.0*066®)  *S*(C*C*C) 

QBAR66a.K)  s  (Qll®  +  Q22®  -  2.0*Q12®  -  2.0*Q66®)*(S*S) 

A  *(C*C)  +  Q66®*((S*S*S*S)  +  (C*C*C*C)) 

C 

C  . CONSTRUCT  THE  REDUCED  STIFFNESS  MA'FRDC  WITH  ROTATION. 

C 

E(1,1.I.K)  =  QBAR11(I.K) 

E(U,IJC)  =  QBAR12aK) 

E(1,3,I,K)  =  0BAR16(I,K) 

E(2,1,I.K)  =  0BAR12(I.K) 

E(2;2,I.K)  =  QBAR22(I.K) 

E(2,3.I,K)  =  QBAR26(IJC) 

E(3,1,I.K)  =  0BAR16(I.K) 

E(3^,I4C)  =  QBAR26(I.K) 

E(3,3,I,K)  =  QBAR66(I,K) 

C 

DO  31  M=  1,3 
D031  J=l,3 

31  CTEMP(M,J,I,K)  =  E(MJ.I,K) 

C 

C  . ...CALCULATE  THE  REDUCED  COMPLIANCE  MATRIX . 

C 

P(l,l)=.-OBARliaK) 

P(1,2)  =  QBAR12(I,K) 

P(1.3)  =  QBAR16a,K) 

P(2,1)  =  QBAR12(I,K) 

P(2,2)  =  QBAR22(I,K) 

P(2,3)  =  QBAR26a.K) 

P(3,l)  =  OBAR16aK) 

P(3.2)  =  QBAR26(I,K) 


P(3.3)  =  QBAR66a,K) 


C 

c 

IA=:3 
IAINV  =  3 
M  =  3 
N=3 
TOL  =  0.0 
C 

CALL  LGINF(P,IAAIJ^,TOL,T,IAINV,S2.WK2,IER) 

C 

XYZ1  =  T(1,1) 

XYZ2  =  T(U) 

XYZ3  =  T(1.3) 

XYZ4  =  T(2,1) 

XYZ5  =  T(2.2) 

XYZ6  =  T(2,3) 

XYZ7  =  T(3,1) 

XYZ8  =  T(3^) 

XYZ9  =  T(3,3) 

C 

D(1.1.I.K)  =  XY7: 

D(1,2,!.K)  =  XYZ2 
D(1.3,I.K)  =  XYZ3 
D(2,1,I,K)  =  XYZ4 
D(2^.I,K;  =  XYZ5 
D(2.3,I,F0  =  XYZ6 
D(3.1,:,K)  =  XYZ7 
D(3^,I,K)  =  XYZ8 
D(3.3,I,K)  =  XYZ9 
C 

C  . CALCULATE  THE  LAMINA  ENGINEERING  PARAMETERS . 

C 

EXXa.K)  =  1.0/(((1.0/E11(I))*(C*C*C*Q)  +  ((1.0/G12a))  - 
A  (2.0*NU12(I)/Elia))) 

B  *(S*S)*(C*C)  +  ((1.0/E22(I))*(S*S*S*S))) 

EYYa,K)  =  1.0/(((1.0/EI  1(I)/(S*S*S*S))  +  ((1.0/G12(I))  - 
A  (2.0*NU12(I)/E11(I))) 

B  *(S*S)*(C*C)  +  ((1 .0/E22(I))*(C*C*C*C))) 

NUXYa.K)  =  EXX(I,K)*((NU12(I)/niia))n(S*S*S*S)  +  (C*C»C*C))  - 
A  ((1.0/£11(I))  + 

B  (1.0/E22(I))  -  (1.0/G12(I)))*(S*S)*(C*Q) 

GXYa,K)  =  1.0/((2.0*((2.0/Elia))  +  (2.0/E22(I))  + 

A  (^.0*NU12(r)/Ell(I))- 

B  (1.0/G12(I))))*((S*S)*(C*C))  +  ((I.0/GI2(I))*((S*S*S*S)  + 

C  (C*C*C*C))^) 

ETAXYXa.K)  =  EXXa.K)*(((2.0/El  1(1))  +  (2.0*NU12(I)/E1 1(1))  - 
A  (1.0/G12a)»*S*(C*C*C)- 

B  ((2.0/F22(I))  +  (2.0*NU12(I)/E1 1(1))  -  (1.0/G12(I)))*<S*S*S)*C) 
ETAXY'/(I.K)  =  EYYa.K)*(((2.0/El  KD)  +  (2.0*NU12(I)/E1 1(1))  - 
A  (1.0/G12(I)))*(S*S*S)^C 
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B  -  ((2.0/E22(I))  +  (2.0*NU12(I)/E11(I))  -  (1.0/G12(I)))*S*(C*C*C)) 

68  CONTINUE 
C 

6  CONTINUE 
C 

C  ....CALCUL  ATE  THE  [A],  [B],  AND  [D]  MATRICES . 

DO  691=  1,3 
D0  69J=  1.3 
A1CI,J)  =  0.0 
BiaJ)  =  0.0 
Dia,J)  =  0,0 

69  CONTINUE 

c 

Z(l)  =  -STHICK/2.0 
C 

DO  701=  1.3 
D0  79J=  1.3 
D0  72K=  IJPLIES 
M  =  MATN(K) 

N  =  ANGN(K) 

2:(K+  1)  =  2(K)  +  THICK(K) 

C 

AT  =  Ea.J.MJ^*(Z(K^•l)  -  Z(K)) 

BT  =  Ea.J>W*(Z^+l)**2  -  Z(K)**2)*(1.0/2.0) 

DT  =  Ea.JJ^J^*(Z(K+l)**3  -  Z(K)**3)*(1.0/3.0) 

C 

Al(I.J)  =  A1(I.J)  +  AT 
BiaJ)  =  Bl(IJ)  +  BT 
Dl(I.J)  =  Dia.J)  +  DT 
72  CONTINUE 
79  CONTINUE 

70  CONTINUE 
C 

DO  9091=  1,3 
D0909  J=  i,3 
DDa.J)  =  Dia.J) 

AAa.j) = Aia.D 
BBa.j)=Biaj) 

909  CONTINUE 

c 

DO  933  I  =  1. ANGLE 

Tl(l.l.i)  =  COS(RTHETA(I))**2 

T1(1.2,I)  =  SIN(RTHETA(I))**2 

T1(1,3,I)  =  COS(RTHETA(I))*SIN(RTHETA(D) 

T1(2.1.I)  =  T1(1.2,I) 

T1(2.2.I)  =  T1(1.1.I) 

Tl(2.3.r  =  -Tl(l,3.I) 

T1(3.;.I)  =T1(1,3,I)*2.0 
T1(3,2.I)=  -T1(3,1,I) 

T1(3!3!i)  =  COS(RTHETA(I))**2  -  SIN(RTHETA(I))**2 
■  933  CONTINUE 
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C 

C 

C  . CALCULATE  THE  THERMAL  COEPHCIENTS . 

C 

D0  332J=  1.NOMAT 
DO  333 1  =  LANGLE 

ALPXai)  =  T1(1.1.I)*ALP1(J)  +  Tl(i;2.I)*ALP2(J) 

ALPY(J,ri  =  T1(2,1.I)*ALP1(J)  +  T1  (2,2,1)*  ALP2(J) 

ALPXY(JJ)  =  T1(3,1,D*ALP1(J)  +  T1(3,2,I)*ALP2(J) 

333  CONTINUE 
DO  4000 1=1,4 

WRITE  (6,555)  ALPX(1,D 
WRITE  (6,555)  ALPY(1,I) 

WRITE  (6,555)  ALPXY(1,I) 

555  FORMAT  (5X,4E1 1.5) 

4000  CONTINUE 
332  CONTINUE 
C 

D0  334L=1J»LIES 
J=MATN(L) 

K=ANGN(L) 

EPST(L,1)  =  ALPX(J,K)*DELTr 
EPST(L,2)  =  ALPY(jr,K)*DELTr 
EPST(L.3)  =  ALPXY(J4C)*DELTT 

334  CONTINUE 

DO  9000 1  =  1,  PLIES 

9000  WRITE  (6,55)  EPSTa,l).EPSTa.2)JEPST(I,3) 

Z(l)  =  -STH!CK/2,0 
C 

DO  335  L  =  14>l  IBS 
J  =  MATN(L) 

K=  ANGN(L) 

DO  3361=  1,3 

TEMF . .  Ea.U,K)*EPST(L,l) 

A  +  E(i,2,J,K)*EPST(L,2)  +  Ea,3,J,K)*EPST(L,3) 

TEMPN  =  TEMP*(Z(L+1)  -  Z(L)) 

TEMPM  =  TEMP*(Z(L+l)**2-Z(L)**2)/2.0 
RNT(I)  =  RNT(I)  +  TEMPN 
RMTG)  =  RMTG)  +  TEMPM 
336  CONTINUE 

335  CONTINUE 
CALL  INVRS(A1,A2) 

C 

C 

C  ....CALCULATE  THE  LAMINATE  COEFFICIENT  OF  THERMAL  EXPANSION . 

C 

IF  (ISYM  .EQ.  1)  GO  TO  921 
DO  4441=  1,3 

ALAMG)  =  (A2a,l)*RNT(l)  +  A2(I,2)*RNT(2)  +  A2(I,3)*RNT(3)) 

A  *(1.0/DELTT) 

444  CONTINUE 


uuwuuuuuuu  uuuuu  uuu 
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921  CONTINUE 


...INVERT  THE  A1  MATRIX  TO  USE  FOR  LAMINATE  ENGINEERING  PARAMETERS.. 
CALL  INVRS(A1  A2) 

. INVERT  THE  B 1  MATRIX  AND  THE  D1  MATRIX  FOR  CALCULATIONS . 


....CALCULATE  THE  LAMINATE  ENGINEERING  PARAMETERS, 

IF  aSYM  .EQ.  1)  GO  TO  200 
EBARX  =  1.0/(STKICK  *  A2(l,l)) 

EBARY  =  1.0/(STHICK  *  A2(2,2)) 

GBARXY  =  1.0/(STHICK  ♦  A2(3.3)) 

NUBARXY  =  -A2(U)/A2(1.1) 

ETABXYX  =  A2(1.3)/A2fl.l) 

ETABXYY  =  A2(2.3)/A2(2;>) 

ETABXXY  =  A2(1.3)/A2(3,3) 

ETABYXY  =  A2(2.3)/A2(3.3) 

XICONTINLTE 

IF  (ILOAD  .EQ.  0)  GO  TO  600 


. COMPUTE  [M]  AND  m  GIVEN  MIDPLANE  STRAINS  AND  CURVATURES.... 

CALL  MPLY(A1.EPSF1.3,3,1) 

CALL  MPLY(B1,KAP.I»2.3,3,1) 

CALL  PLUS(P1,P2.PN.3,1) 

CALL  MPLY(Bl,EPSF1.3,3,l) 

CALL  ^LY{D1,KAPJ?2,3,3,1) 

CALL  PLU5i(Pi4^JPM,3.1) 

GO  TO  300 
10  CONTINUE 

...COMPUTE  THE  MIDPLANE  STRAINS  AND  CURVATURES  GIVEN  [M]  AND  [N].... 

IF  aSYM  .EQ.O)  THEN 
CALL  INVRS(A1,A3) 

CALL  INVRS(D1,D2) 

DO  4001=  1.3 
DO  401  J=  1,3 
B3(I4)  =  0 
B4(U)  =  0 
401  CONTINUE 
400  CONTINUE 
ELSE 

CALL  INVRS(A1,A4I 
CALL  MPLY(  A4,B1.B2,3,3,3) 
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DO  402 1  =  1,3 
DO  403  J=  1,3 

B2aJ)  =  -B2a.J) 

403  CONITNUE 
402  CONTINUE 

C 

CALL  MPLY(B1,A4,B5,3,3,3) 

CALL  MPLY(B1,B2,D3,3,3,3) 

CALL  PLUS(D1,D3,D4,3,3) 

CALL  INVRS(D4  J12) 

CALL  MPLY(B2,D2,B3.3,3,3) 

CALL  MPLY(B3,B5,D3,3,3,3) 

CALL  MINUS(A4,D3,A3,3,3) 

CALL  MPLY(D2,B5,B4,3,3,3) 

DO  4041=1,3 
DO405J=  1,3 
B4aJ)  =  -B4a.J) 

405  CONTINUE 

404  CONTINUE 
ENDBF 

C 

CALL  MPL  Y(  A3  JRNJ»1, 3,3,1) 

CALL  MPLY(B3,RMJ>2„3,3,1) 

CALL  PLUS(P1  J>2,EPS,3.1) 

CALL  MPLY(B4,RN,P1.3,3,1) 

CALL  MPLY(D2JIM,P2,3,3,1) 

CALL  PLUS(P1  J>2,KAP,3,1) 

C 

C  ...CALCULATE  MIDPLANE  STRAINS  AND  CURVATURES  WITH  THERMAL  LOADS... 
C 

CALL  MPLY(A2J1NTJ»1, 3,3,1) 

CALL  MPLY(B3,RMT4»2,3,3,1) 

CALL  PLUS(P1  J»2,EPSTH,3,1) 

Tvitl-  I  ynHyisj.'i  i  ,r  i ,  J,  J,  i) 

CALL  MPLY(D2RMT,P2,3,3,1) 

CALL  PLUS(P1,P2,KAPT,3,1) 

C 

C  ..CALCULATE  MIDPLANE  STRAINS  AND  CURVATURES  WHH 

C  MECHANICAL  AND  THERMAL  LOADING . 

C 

DO  408  I  =  1,3 
RNB(I)  =  RN(I)  +  RNT(I) 

RMB(I)  =  RM(I)  +  RMT(I) 

408  CONimJE 
C 

CALL  MPLY(A2RNB, PI, 3,3,1) 

CALL  MPLY(B3,RMB,P2,3,3,1) 

CALL  PLUS(Pl,P24iPSB,3,l) 

CALL  MPLY(B4,RNB,P1,3,3,1) 

CALL  MPLY(D2RMB,P2,3,3,]) 

CALL  PLUS(P1  J>2,KAPB,3,1) 
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C 

300  CONTINUE 
C 

C  ....CALCULATE  THE  SFRESSES  IN  EACH  PLY  IN  THE  X-Y  DIRECnON..,. 
C 

DO9G0  L=  LPLIES 
DO901  J=l,3 

H = (z(L+i)-zao)/2.o + za) 

WRITE  (645)  H 

EPSX(L,J)  =  EPS(J)  +  H*KAP(J) 

EPSXT(L.J)  =  EPSTH(J)  +  H-^KAPT(J) 

EPSXB(L  J)  =  EPf  B(J)  +  H*KAPB(J) 

EPSTIOLJ)  =  EP'.XT(L4)  -  EPST(L,J) 

EPSTH1(LJ)  =  EPSXB(L.J)  -  EPST(LJf) 

901  CONTINUE 
900  CONTINUE 
DO  9001  I=1J‘LIES 

WRITE  (645^  EPSXT(I,l)£PSXTa;i),EPSXT(I,3) 

9001  CONTINUE 

DO  900/ 1  =1^LIES 

WRITE  (6.5.5)  EPST1(I,1),EPST1(I.2)JEPST1(I,3) 

9002  CONTINUE 

WRITE  (645)  RNT(1)RNT(2).RNTIX3) 

WRITE  (645)  EPSTH(1)JEPSTH(2)EPSTH(3) 

WRITE  (645)  ICAPT(1).KAPT(2).KAP1(3) 

C 

DO  701 1  =  i  JULIES 
K  =  ANGN(I) 

J  =  MATN(I) 

D0  7C3L=1,3 

SIGXrT,L)  =  E(L.1.J,K)*EPSX(I,1)  +  E(L.2.I,K)*EPSX(I,2)  + 

A  E(L,3.J,K)*EPSXa.3) 

SIGXra-L)  =  ^(L,U,K)-'EPSTia.l)  +  E(L,24,K)*EPST1(I,2)  + 

A  E(L,3,I.K)-^EPST1(I.3) 

SIGXB(IL-)  =  E(l.,l.J,K)*EPSTHl(I,l)  +  E(L;’,J,K)*EPSTH1(I,2)  + 

A  E(L,34.K)*EPSTH1(I,3) 

703  CONTINUE 
701  CONTINUE 
C 

C  ....CALCULATE  THE  SI  RESSE.S  IN  THE  1  2  DIRECnON . 

C 

DO  932  I  =  1. ANGLE 

Tld.l.I)  =  COS(RTHETA(I))*'^2 

1 1  (1.2, 1)  =  SIN(RTHETA(I))**2 

Tl(]  ,3,1)  =  COS(PTHETA(I))*SlN(RTHETA(I)) 

T1(2,1.I)  =  T1(1,2,I) 

T1(2,2.I)- Tld.l.I) 

T1  (2,3,1)  =  -Tl  (1,3,1) 

T1(3,1,I)  =  T1(1.3,I)*2.0 
T1(3.2,1)  =  -T1(3.1.I) 

T1(3,3,I)  =  COS(RTHETA(I))**2  -  SIN(RTHETA(I))**2 


932  CONTINUE 
C 

DO  963  L=  1, PLIES 
K  =  ANGN(L) 

D0964  1=1,3 

EPS1(L4)  =  Tl(I.l,K)*EPSXa-,l)  +  Tia^.K)*EPSX(L.2) 

A  +Tia,3.K)*EPSX(L.3) 

EPS1T(LJ)  =  Tia.l.K)*EPSXT(L,l)  +  T1(I,2.K)*EPSXT(L.2) 

A  +T1(I,3,K)*EPSXT(L,3) 

EPS1B(L,0  =  Tl(I,l.K)*EPSXB(L.l)  +  T1(I.2,K)*EPSXB(L,2) 

A  +T1(I,3.K)''EPSXB(L,3) 

964  CONTINUE 
963  CONTINUE 
C 

DO  705  1=  1, PLIES 
K  =  MATN(I) 

SIGia.1)  =  Qll(K)*EPSl(I.l)  +  Q12(K)*EPS1(I,2) 

SIGia.2)  =  QnrtO^EPSia.l)  +  Q22(K)*EPS1(I,2) 

SIGia.3)  =  Q66(K)*£PSia.3) 

SIGlTa.1)  =  QllCQ*fc’PSlT(I.l)  +  Q12(K)*EPS1T(I.2) 

SIGlTa,2)  =  Q12(K)*EPS1T(I.1)  +  Q22(K)*EPSlTa,2) 

SIG1T(I,3I  =  Q66(K)*EPS1T(I,3) 

C 

SIGlBa.l)  =  011(K/EPS1B(I.1)  +  Q12(K)*EPS1B(I,2) 

SIG1B(I.2)  =  Q12(K)‘'EPS1B(I.1)  +  Q22(K)*EPS1B(I,2) 

SIGlBa,3)  =  Q66(K)*EPS1B(I,3) 

705  CONTINUE 
C 

c 

c 

C  . DISPLAY  ALL  DATA  OBTAINED . 

C 

PRINT  *,  ’THIS  IS  A  LAMINATE  COh.-IPOSITE  MATERIAL  PROGRAM’ 

PRINT  *■,  ’  ’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’THERE  IS  ’,  NOMAT,’  TYPE  OFMATEPJAL’ 

PRINT  *,  ’  ’ 

PRINT  ’THERE  ARE  ’.ANGLE,’  DIFFERENT  ANGLES  IN  THIS  LAMINATE’ 
PRINT  *,  ’  ’ 

DO  210 1  =  1. ANGLE 
PRINT  ’  ’ 

PRINT  *,  ’ROTATION  EQUAI.S  ’,THETrt(I),’  DEGREES’ 

210  CONTINUE 
PRINT  *,  ’  ’ 

PRINT  *.  ’  ' 

IF(ITEMP.EQ.  l)GOTO  HI! 

PRINT  *,  ’THERE  ARE  THERMAL  EFFECTS  IN  THIS  LAMINATE’ 

PRINT  ’ 

PRINT  *,  ’  ’ 
nil  CONTINUE 
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DO  201 1  =  1  J^OMAT 

PRINT  *,  ’THE  ENGINEERING  PARAMETERS  FOR  MATERIAL  ’,1,’  ARE:’ 
PRINT  *.  ’  • 

PRINT*’  Ell  E22  G12  NU12  ’ 

PRINT*  ’  ’ 

WRITE  (6,52)  Ell(I)£22(I),G12a)JW12(I) 

52  FORMAT  (3E11.5,1F5.2) 

201  CONTINUE 
PRINT  *.  ’  ’ 

PRINT  *,  ’  ’ 

DO  702 1  =  1 jsrOMAT 
PRINT  *,  ■  ’ 

PRINT  *,  ’ORTHOTROPIC  COMPLIANCE  MATRIX  FOR  MATERIAL 
PRINT  *,  ’  ’ 

DO10J=  1,6 

10  WRITE  (6,50)  (A(J.K,I).  K  =  1,6) 

50  FORMAT  (5X,  6E1 1.5) 

PRINT  *,  ’  ’ 

202  CONTINUE 
PRINT  *,  ’  ’ 

DO  203  I  =  1  J^OMAT 
PRINT  *,  ’  ’ 

PRINT  *,  ’ORTHOTROPIC  STIFFNESS  MATRIX  FOR  MATERIAL  ’,! 
PRINT  *,  ’  ’ 

DOll  J=  1,6 

1 1  WRITE  (640)  (B(J  is:, I),  K  =  1,6) 

PRINT  *,  •  ’ 

203  CONTINUE 
PRINT  *,  ’  ’ 

DO  204  I  =  1, NOMAT 
DO  205  K=l, ANGLE 
PRINT  *,  ’  ’ 

PRINT  *,  ’THE  REDUCED  STIFFNESS  MATRIX  FOR  MATERIAL  ’,  I 

nn  A  9  A  ^  I  f*N  Tk  A  T  •  -WW  •  • 

,1V,  lo; 

PRINT  *,  ’  ’ 

DO  12  J=  1,3 

12  WRrrE(641)  (E(J,M,I,K),  M=l,3) 

51  FORMAT  (5X,3E1 1.5) 

PRINT  *,  ’  ’ 

205  CONTINUE 

204  CONTINUE 
PRINT  *,  ’  ’ 

C 

DO  206  I  =  1, NOMAT 
DO  207  K=  l,.\NGLE 
PRINT  *,  ’  ’ 

PRINT  *,  ’THE  REDUCED  COMPLIANCE  MATRIX  FOR  MATCRIAL  ’,  I 
PRINl'  *,  ’AND ROTATION  ’,  K,  ’  IS:’ 

PRINT  *,  ’  ’ 

DO  13  J=  1.3 

13  WR1TE(641)  (D(J,M,1,K),  M  =  1,3) 
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PRINT  *, '  ’ 

207  CONTINUE 
206  CONTINUE 

PRINT  *,  ’  ’ 

PRINT  *,  ’  ’ 

1 

DO208I=  IJ^OMAT 
DO209K=  l^GLE 
PRINT  *,  ’  ’ 

PRINT  *.  *THE  LAMINA  ENGINEERING  PARAMETERS  FOR  MATERIAL  M 
PRINT  *,  ’AND  ROTATION  ’.K,’  ARE:* 

PRINT  *.  ’  ’ 

PRINT*  ’  EXX  EYY  NUXY  GXY  E^i'AX^’X  ETAXYY* 

PRINT  *,  ’  ’ 

WRITE(6^3)  EXXajC)JiYY(I,K)^XYa.K).GXya.K)^AXYXa,K), 

A  ETAXYYa,K) 

53  FORMAT  (2E1 1.5,F5.2,E1 1.5.2F6.3) 

PRINT  *.  ’  ’ 

209  CONTINUE 

208  CONTINUE 
PRINT  *,  •  * 

IF  aSYM  .EQ.  1)  GO  TO  800 
PRINT  *,  ’  ’ 

PRINT  *,  ’THE  LAMINATE  ENGINEERING  PARAMETERS  ARE:’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’EBARX  EBARY  GBARXY  NUBARXY  ETABXYX  ETABXYY 
AETABXXY  ETABYXY’ 

PRINT*  ’  ’ 

WRITE(644)  EBARX^BARY.GBARXY.NUBARXY^En'ABXYX, ETABXYY, ETABXXY, 
A  ETABYXY 

54  FORMAT  (3E9.3,F5.2,4E10.4) 

PRINT  *,  ’  ’ 

GO  TO  801 

800  CONTINUE 

PRINT  *.  ’THERE  ARE  NO  LAMINATE  ENGINEERING  PARAMETERS’ 

PRINT  *,  ’  ’ 

801  CONTINUE 

PRINT  *,  ’THE  LAMINATE  COEFFiaENTS  OF  THERMAL  EXPANSION  ARE:  ’ 

PRINT  *.  ’  ’ 

DO  881  =  1.3 
WRITE  (6,58)  ALAM(I) 

58  FORMAT  (5X,3E1 1.5) 

88  CONflNUE 
PRINT  *,  ’  ’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’THE  [A]  MATRIX  FOR  THE  LAMINATE  IS’ 

PRINF  *,  ’  ’ 

DO  211  1=  1.3 

21 1  WRITE  (6,51)  (AA(I,J),  1=1,3) 

PRINT  *,  ” 


PRINT  *,  ’  ’ 

PRINT  ♦.  *THE  [B]  MATRIX  FOR  THE  LAMINATE  IS’ 

PRINT  ’  ’ 

DO  212  I  1,3 

212  WRITE  (6,51)  (BB(IJ).  J  =  1,3) 

PRINT  ♦.  ’  ’ 

PRINT  •,  ’THE  [D]  MATRIX  FOR  THE  LAMINATE  IS’ 

PRINT  *,  ’  ’ 

DO  213 1=1,3 

213  WRITE  (6,51)  (DDHJ).  J=-T.3) 

PRINT  *,  ’  ’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’THE  MECHANICAL  STRESSES  IN  THE  X-Y  DIRECTION  ARE’ 
PRINT  *,  ’  ’ 

DO  214 1  =  13»LIES 
PRINT*,  ’IN PLY  ’,! 

WRITE  (6,55)  (SIGX(I.J),  J=l,3) 

214  CONTINUE 
PRINT  *,  ’  ’ 

PRINT  *,  ’THE  THERMAL  S'rRESSES  IN  THE  X-Y  DIRECTION  ARE’ 
PRINT  *,  ’  ’ 

DO  218  I  =  1  JULIES 
PRINT  *, ’D,’ PLY ’J,’:’ 

WRITE  (6,55)  (SIGXT(I,J),  J=l,3) 

218  CONITNUE 
PRINr  •.  ’  ’ 

PRINT  *,  ’THE  THERMAL  AND  MECHANICAL  STRESSES  ’ 

PRINT  *,  ’  IN  THE  X-Y  DIRECHON' 

PRINT  *,  ’  ’ 

DO  219 1  =  1,  PLIES 
PRINT  *,  ’IN  PLY ’,1,’:’ 

WRITE  (6,55)  (SIGXB(U),  J=l,3) 

219  CONTINUE 
PRINT  *,  ’  ’ 

55  FORMAT  (10X,3E11.5) 

PRINT  *.  ’  ’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’THE  MECHANICAL  STRESSES  IN  THE  1-2  DIRECTION  ARE’ 
PRINT  *:  ’  ’ 

DO  215  I  =  LPLIES 
PRINT  *,  ’IN  PLY 
V/RITE  (6,55)  (SIG1(I,J),J=1,3) 

215  CONTiNlTi 
PRINI’  *,  ’  ’ 

PRINT  *.  ’THE  'rHE;RMAL  STRESSES  IN  THE  1-2  DIRIZCTION  ARE’ 
PRINT  *,  ’  ' 

DO  220 1  =  LPLIES 
PRINT  *,  ’IN  PLY  ’.I,’:’ 

WRITE  (6,55)  ISiaiT(I,J),Jr:l,3) 

220  CONTINUE 
PRINT  *,  ’  ’ 


PRINT  *,  ’THE  THERMAL  AND  MECHANICAL  STRESSES  ’ 

PRINT  *.  ’  IN  THE  1-2  DIRECnON  ’ 

PRINT  *,  ’  ’ 

DO  221 1  =  1  J>LIES 
PRINT*  ’IN PLY ’.I,’:’ 

WRITE  (6^5)  (SIGlBaJ).J=1.3) 

221  CONTINUE 
PRINT  *,  ’  ’ 

PRINT  *.’  ’ 

IF  (ILOAD  .EQ.  0)  GO  TO  775 
PRINT  *.  ’  ’ 

PRINT  *.  ’THE  RESULTANT  FORCES  AND  MOMENTS  ARE’ 

PRINT  *,  ’  ’ 

PRINT  *,  ’  [N]  [M]’ 

DO  216 1  =  1,3 

216  WRITE  (6,57)  PN(D.PM(I) 

57  FORMAT  (F7.24X,F7.2) 

PRINT  *.  ’  ’ 

GO  TO  776 

775  CONTINUE 
PRINT  *,  ’  ’ 

PRINT  *.  ’THE  RESULTANT  MIDPLANE  STRAINS  AND  CURVATURES  ARE’ 
PRIf^T  *,  ’  ’ 

PRINT  ♦.  ’STRAINS  CURVATUP^S’ 

PRINT  *,  ’  ’ 

DO  217 1  =  1.3 

217  WRITE  (6,56)  EPS(r;,KAPa) 

56  FORMAT  (El  1.5.5XE1 1.5) 

776  CONTINUE 
C 

PRINT  *,  ’  ’ 

PRINT  *,  ’  ’ 

PPJNT  *,  ’THE  RESULTS  OF  THE  MAX  STRESS  FAILURE  CRITERION’ 
PRINT  *,  ’  ’ 

IF  (QUEST  .EQ.  1)  GO  TO  2000 
DO  2001 1  =  1,  PLIES 

PRINT  *,  ’THE  FORCE  RESULTANT  TO  PRODUCE  FIRST  PLY  FAIL’ 

PRINT  ♦.  ’  ’ 

PRINT  *,  'FOR  PLY  ’,  I,’;’ 

2001  WRITE  (6,55)  (NXX(IJ),  J=I,3) 

GOTO  2002 
2000  CONTINUE 

PRINT  *,  ’THE  CHANGE  IN  TEMPERATURE  TO  PRODUCE 
A  FIRST  PLY  FAILURE  IS’ 

PRINT  *,  ’  ’ 

DO  2003  I  =  14>LIES 
PRINl’  *,  ’FOR  PLY 
2003  WRI'rS  (6,55)  (DEL(I,J),  J=l,3) 

20i)2  CONT  INUE 
Pi’INT  *,’  ’ 

PRO'n'  *,  ’  ’ 
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WRITE  (645)  EPST1(1,1),EPSTH(1),KAPT(1) 

WRITE  (645)  EPSXT(1,1),EPST(1,1),ALPX(1.1) 

WRITE  (645)  T2(l,2.1).ALP2(l)J)ELTr 
WRITE  (645)  RNT(1)  JIMT(1),TEMPN 
PRINT*,” 

DO  3000 1  =  1,  PLIES 
WRITE  (645)  (ALPX(I4),J=1,3) 

WRITE  (645)  (EPSTa.J)J=1.3) 

WRITE  (645)  (ALPY(IJ),J=1.3) 

3000  CONTINUE 
PRINT  *,  ’  ’ 

C 

STOP 

END 

c 

c  . SUBROUTINE  TO  INVERT  A  3  B  Y  3  MATRIX . 

C 

SUBROUTINE  INVRS(U.V) 

C 

DIMENSION  U(3.3).  A(3.3).  V(3,3) 

DETl  =  U(2,2)*U(3,3)-U(2,3)*U(3.2) 

DET2  =  U(2.1)*U(3.2)-U(2,3)*U(3,1) 

DET3  =  U(2.1)*U(3,2)-U(2.2)*U(3,1> 

DET  =  U(1,1)*DET1-U(1,2)*DET2+U(1.3)*DET3 
C 

A(l,l)  =  U(2,2)*U(3.3j-U(2,3)*U(3A) 

A(l,2)  =  (U(1,2)*U(3,3)-U(1,3)*U(3,2))*(.1.0) 

A(l,3)  =  U(1,2)*(2.3)-U(1.3)*U(2.2) 

A(2,l)  =  (U(2.1)*U(3.3)-U(2.3)*U(3.1))*(-1.0) 

A(24)  =  U(1.1)*U(3.3)-U(1,3)*U(3,1) 

A(2,3)  =  (U(1.1)*U(2.3)-U(1,3)*U(2,1))*(-1.0) 

A(3.1)  =  U(2,1)*U(3.2)-U(2,2)*U(3,1) 

A(3,2)  =  (U(1,1)*U(3  2)-U(l,2)*U(3.1))*(-1.0) 

A(3.3)  =  U(1.1)*U(24)-U(1,2)*U(2.1) 

C 

DO  20,1=  1,3 
DO  10,  J=  1,3 
va  J)  =  A(U)/DET 
10  CONTINUE 
20  CONTINUE 
RETURN 
END 
C 

c 

SUBROUTINE  MPLY(U,V,T,M,K,N) 

C 

DIMENSION  U(M,K),V(K,N),T(M,N) 

DO  20,  1  =  14^ 

DO  20,  J=  141 
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Taj)=o.o 
DO20,L=  1,K 
T1  =  U(I,L)*V(L,J) 
TGJ)  =  TGJ)  +T1 
20  CONTINUE 
RETURN 
END 


C 

c 

SUBROUTINE  PLUS(U.V.T.K  J.) 

C 

DIMENSION  U(KW.V(KJL).T(KJL.) 
DOS.  I  =  1.K 
DOS,  ,T  =1,L 
Ta.J)  =  0.0 
5  CONTINUE 
DO  10,  I  =  l.K 
DO  10.  J  =  1.L 
TCIJ)  =  U('.J)  +  vaj) 

10  CONTINUE 
RETURN 
END 
C 

c 

SUBROUTINE  MINUS(U,V.T.K,L) 
C 

DIMENSION  U(K,L).V(KJL).T(Ki) 
DOS,  I  =  1.K 
DOS  J  =1,L 

Ta,J)  =  0.0  ^ 

5  CONTINUE 
DOlO,I=l,K 
DO10,J  =  l.L 
TaJ)  =  U(IJ)-Va.J) 

10  CONTINUE 
RETURN 
END 


Appendix  B 


Session  File 


B.l  Lo{;  In 


>PATRAN2 

>4107 

>GO 

>1 


B.2  Mode!  Generation 


>1 

>1 

>GR.i 

>GR.8.TR.0/9/0 

>VIEW 

>1 

>30,30,30 

>GR,4„0/3/D 

>OR,2,,0,.6,.8 

>GR,3„0,1. 1,1.3 

>GR,5,.0,1.6,1.6 

>GR,6„0,2.1,1.8 

>GR,7„0,4,1.9 

>GR,9,.0,6.6,1.1 

>LI,lJTr„14/3/5/6/4 

>LI,2,TR,80,1 

>LUJ^„4^/9/8 

>LI,4.TR,80,3 

>PA,UL„1,2 

>PA,22L.,3,4 

>PA3JvlI,Z,l 

>PA,4J^,Z,2 

>SETd.ABE,OFF 

>END 


B.3  Node  Generation 


>2 

>1 


>1 

>SET,CPLOT,ON 

>GF4>Al„4/25 

>GFJ?A2„4/25 
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>GF4»A3.;j5/4 

>GFJ»A4,^/4 

B.4  Material  Property  Definition 

r>PHAT,10,ORT^6.3E6^(l.lE6)..41..45..C0801.;2.2E6. 

>.5E6;<i.2E6r5E-6;2(6.5E-6) 

>PMAT.30,ISO..67E6„.36.^5E-6 

>DATA,15,.6,.4,1.1,2,2,1 

>PMAT.50.HALJ)15.10/30 

>PMAT.100/107J.AMS,4(,015).0/30/15/90.4(50) 

>END 


B.5  Element  Mesh  Generation 

>2 

>1 

>CFJ»A1T4,TRI.'3/53,T1 

>EMD 

B.6  Equivalencing 

1  1  ,Orr 

>3 

>N 

>2 

>1 

>N 

END 

B.7  Force  and  Constraint  Application 

>SET4’Hl,OFF 

>4 

>1 

>DFJ»A1.FORCE,100,1.ED2 
>DFJ>A2JFORCE,100,1  £D2 
>DFJ>A3JFORCE,100,1JED1 
>DFPA44^ORCF.,100.1JED  1 

>df,pai  ,d;oa)///.i  Wloo 
>END 

B.8  Model  Optimization 

>7 

>5 

>3 

>2 

>END 

B.9  Material  Property  Assignment 

>5 

>1 
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>PFJ>AlT4.'nU/3/53,1.106 

>END 

B.10  ANSYS  Interface 

>8 

>5 

>3 

>1 

>1 

>’ 

>BLADE1 

>Y 

>Y 

>2 

>2 

>5 

>6 


u  u  u 
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Appendix  C 


Equation  Solution  Program 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM  CALC 


*  PROGRAM  CALC 

* 

*  PATRICK  GRAHAM  FORRESTER 

*  UNIVERSITY  OF  VIRGINIA  1989 

*  THESIS  RESEARCH 

* 


THIS  PROGRAM  IS  DESIGNED  TO  READ  THE  DISPLACEMENT  INFORMATION 
PROUCED  BY  ANSYS.43A  FOR  24  NODES  OF  THE  COMPOSITE  BLADE 
MODEL  AND  USE  IT  TO  CAI.CULATE  THE  CHANGE  IN  DISPLACEMENT 
OVER  THE  LENGTH  OF  THE  BLADE  AND  THEN  SOLVE  THE  CLASSICAL 
BEAM  EQUATIONS. 


TIT?  AT  T  rvi  T  T  Txr*  T  T  yrwn  t 


REAL  UX(2,100).UY(2,100),UZ(2,1(X)) 

REAL  INT, SLOPE 

REAL  Y1,Y2,Z1.Z2,L1^2,T1.T2 

REAL  BErA(100),PHI(100),PSI(100).U(100) 

REAL  UXNOT(100),UYNOT(100),UZNOT(100) 

REAL  BETAP(1(X))T>HIP(100),PSIP(100),UP(100) 
DIMENSION  A(4,4,I2),B(4),X(4) 

INTEGER  BLADES, NODESiOADS 
REAL  H.C,ZETA(100) 

OPEN(UNtT=12,FILE=’NEWDATA’.STATUS=’OLD’) 
OPEN  (UNTr=13^FILE=’RES5^STATOS=’NEW’) 


ASSIGN  INPUT  DATA 


H  =  3.33333 
C  =  9.0 
NODES  =  2 


BLADES  = 1 
LOADS  =  4 
C 

C  . READ  IN  NODAL  DISPLACEMENTS  FROM  ANSYS. 

C 

C  1  =  LEADING  EDGE  NODE 
C  2  =  LEADING  EDGE  NODE 
C 

DO  1  K=  1, blades 
D0  9L-=1,L0ADS 
DO  2  J  =  1,3*N0DES 
C 

T1  =  9.0 
T2  =  0.0 
LI  =  0.0 
L2  =  0.0 
C 

READ(12.*)  UX1.UY1.UZ1.UX2.UY2,UZ2 
C 

UX(1.J)  =  UX1 
UY(1.J)  =  UYl 
1^1  J)  =  UZl 
UX(2,J)  =  UX2 
UY(2,J)  =  UY2 
UZ(2  J)  =  UZ2 
C 

C  . CALCULATE  ZETA  FOR  ISOTROPIC  MODEL . 

C 

IF(K.NE.  1)  GO  TO  2 
IF(L.NE.  l)GOT02 
C 

Y1  =L1  +UY1 
Y2  =  T1  +UY2 
Z1  =  L2  +  UZl 
Z2  =  T2  +  UZ2 
C 

SLOPE  =  (Y1  -  Y2)/(Z1  -  Z2) 

INT  =  (Yl-(Zl/SLOPE)) 

ZETA(J)  =  INT/C 
C 

2  CONTINUE 
C 

c 

C  . CALCULATE  DISPLACEMENTS  AND  ROTATIONS 

C 

C  . CALCULATE  ALL  BETA  FOR  MODEL . 

C 

D0  3  J  =  1,6 

BETA(J)  =■■  -(UZ(1,J)  -  UZ(2,.D)/C 

3  CONTINUE 
C 


on  n  noon  n  o  o  noon  noon  o  on 


■CALCULATE  ALL  PHI  FOR  MODEL, 


D04J=  1,6 

UZNOT(J)  =  UZ(1J)*(1  -  ZETA(J))  +  UZ(2J)*(2E'rA(J)) 
UYNOT(J)  =  UY(1.J)*(1  -  ZET  A(J))  h-  UY(24)*(ZETA(J)) 
UXNOT(J)  =  UX(1,J)*(1  -  ZETA(J))  +  UX(2J)''(ZETA(J)) 
4  CONTINUE 

PHI(l)  =  (UZNOT(3)  -  UZN0T(2))/H 
PHI(2)  =  (UZNOT(2)  -  U2NOT(l))/H 
PHI(3)  =  (UZNOT(6)  -  U2NOT(5))/H 
PHI(4)  =  (UZN0T(5)  -  UZN0T(4))/H 

. CALCULATE  ALL  PSI  FOR  MODEL . 


PSI(l)  =  (UYNOT(3)  -  UYNOT(2))/H 
PSI(2)  =  (UYNOT(2)  -  UYNOT(l))/H 
PSI(3)  =  (UYNOT(6)  -  UYNOT(5))/H 
PSI(4)  =  (UYNOT(5)  -  UYNOT(4))/H 


. CALCULATE  THE  FIRST  DERIVATIVE  WRT  X, 

BETAP(l)  =  (BETA(3)  -  BETA(1))/(2*H) 

BETAP(2)  =  (BETA(6)  -  BETA(4))/(2*H) 

PHIP(l)  =  (PHI(l)  -  PHI(2))/H 
PHIP(2)  =  (PHI(3)  -  Pin(4))/H 

PSIP(1)  =  (PSI(l)  -  PSI(2))/H 
PSIP(2)  =  (PSI<3)  -  PSI(4))/H 

Up/^n  ^  UXNOTG''  -  UXNOT'^UV'^2*K' 

UP(2)  =  (UXn6t(6)  -  UXNOT(4))/(2*H) 

DO30J=  IJ^ODES 

. IDENTIFY  LOADS  AND  DISPLACEMENTS.... 

X(l)=  1000 
X(2)  =  1000 
X(3)  =  1000 
X(4) =  1000 

B(l)  =  BETAP(J) 

B(2)  =  PHIP(J) 

B(3)  =  PSIP(J) 

B(4)=  UPu) 

. SOLVE  THE  EQUATION  MATRIX . . 


n  o  n  ^  n  ri  n 
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C 

DO  99 1  =  1,4 
Aa.LJ)  =  B(I)/X(L) 

99  CONTINUE 

. WRITE  THE  COLUMN  SOLUTION  FOR  FITI URE  USE . 

WRITE(13,*) 

WRrrE(13,*)  'NTJMBER  COLUMN  FOR  SECTION:  V,’,*, '  BLADE: 
A  K 

WRrrE(13.*) 

D06 1  =  1,4 
WRnE(13,666)  A(1,L,J) 

6  CONTINUE 

FORMAT  (15X,1E15,6) 

. BUILD  THE  M/.  IRIX  FOR  EACH  NODE . 

30CONTI>:UE 

9  CONTINUE 
DOll  J=  IJSTODES 
WRrrE(13,*) 

WRITE  (13,*)  ’MATRU  FOR  SECHON:  y,\’,  ’  BLADE:  ’,K 
WRITEdS,*) 

DO  101=  1,4  , 

WRITE  (13,667)  A(I,1,J),  Aa,2J),  A(I,3,J).  Aa,4,J) 

10  CONTINUE 

11  CONTINUE 

667  FORMAT  (12X,4E14.6) 

DO300J=  1,6 
WRITE  (13,669)  ZETA(J) 

669  FORMAT  (2X,1E14.6) 

300  CONTINUE 
1  CONTINUE 
STOP 
END 


I 


Appendix  D 


Sample  ANSYS  Input  File 


C***  PREP?  INPUT  PRODUCED  BY  "PATANS"  VERSION  1.7B 
C***  06-FEB-89  16:28:  9 


/ITTLE,  BLADE2MZ 
/NOPR 

C***  NODAL  COORDINATE  DEFINTHON 


N, 

1. 

0.0000000 

.  0.0000000 

O.OOOOOv... 

N. 

^  » 

0.0000000 

.  0.7652769 

r 

0.9910634 

N, 

3. 

0.0000000 

,  1.788908 

« 

1.688653 

N. 

4. 

0.0000000 

,  3.000000 

f 

2.000000 

N. 

5. 

3.333333 

.  0.0000000 

* 

0.0000000 

N, 

6. 

3.333333 

,  0.7652769 

9 

0.9910634 

N, 

7, 

3.333333 

,  1.788908 

9 

1.688653 

N. 

8, 

3.333333 

,  3.000000 

9 

2.000000 

N. 

9. 

6.666666 

.  0.0000000 

9 

0.0000000 

N. 

10, 

6.666666 

.  0.7652769 

9 

0.9910634 

N, 

11. 

6.666566 

,  1.788908 

9 

1.688653 

N, 

12. 

6.666666 

,  3.000000 

9 

2.000000 

N, 

15 , 

lu.uuuOO 

,  O.uuuuuuu 

9 

U.COOOOOO 

N. 

14, 

10.00000 

,  0.7652769 

9 

0.9910634 

N, 

15, 

10.00000 

,  1.788908 

9 

1.688653 

N, 

16. 

10.00000 

,  3.000000 

9 

2.000000 

N, 

17, 

13.33333 

,  0.0000000 

9 

O.IXKXXXX) 

N, 

18, 

13.33333 

,  0.7652769 

9 

0.9910634 

N. 

19. 

13.33333 

,  1.788908 

9 

1.688653 

N, 

20. 

13.33333 

,  3.000000 

9 

2.000000 

N. 

21. 

16.66666 

,  0.0000000 

9 

0.0000000 

N, 

22, 

16.66666 

,  0.7652769 

9 

0.9910634 

N. 

23. 

16.66666 

,  1.788908 

9 

1.688653 

N. 

24. 

16.66666 

,  3.000000 

9 

2.000000 

N, 

25, 

20.00000 

.  0.0000000 

9 

0.0000000 

N, 

26. 

20.00000 

,  0.7652769 

9 

0.9910634 

N. 

27. 

20.00000 

,  1.788908 

9 

1.688653 

N, 

28, 

20.00000 

.  3.000000 

9 

2.000000 

N. 

29, 

23.33333 

,  0.0000000 

9 

0.0000000 

N. 

30, 

23.33333 

,  0.7652769 

9 

0.9910634 
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N, 

31. 

23.33333 

,  1.788908 

,  1.688653 

N, 

32. 

23.33333 

,  3.000000 

.  2.000000 

N. 

33. 

26.66666 

,  0.0000000 

.  0.0000000 

N. 

34. 

26.66666 

,  0.7652769 

.  0.9910634 

N. 

35. 

26.66666 

.  1.788908 

.  1.688653 

N. 

36, 

26.66666 

,  3.000000 

,  2.000000 

N. 

37, 

30.00000 

.  0.0000000 

.  0.0000000 

N. 

38. 

30.00000 

,  0.7652769 

.  0.9910634 

N. 

39. 

30.00000 

,  1.788908 

.  1.688653 

-THRU- 


N, 

301  , 

79.99998 

1 

.7652769 

.  O.OOOOCOO 

N, 

302, 

79.99998 

• 

1.788908 

,  0.0000000 

N, 

303, 

79.99998 

« 

3.000000 

,  0.0000000 

N, 

304, 

79.99998 

• 

5.085196 

.  0.0000000 

N, 

305, 

79.99998 

• 

7.079720 

.  0.0000000 

N, 

306, 

0.000000 

1 

.7652769 

.  0.0000000 

N, 

307, 

0.000000 

• 

1.788908 

,  0.0000000 

N, 

308, 

0.000000 

t 

3.000000 

,  0,0000000 

N, 

309, 

0.000000 

t 

5.085196 

,  0.0000000 

N, 

310, 

0.000000 

$ 

7.079720 

,  0.0000000 

C***  ELEMENT  LIBRARY  DEMOTION 


ET,  1,  53,  0,  0,  8.  0,  0,  0 

C***  ELEMENT  CONNECTIVrrY  DEFINITION 

TYPE.  1  $MAT,  106  $REAL,  1 

E.  106,  103,  252,  252 

E.  i02,  103,  106,  106 

r? 

AJ  ,  IKAJt  iUU 

E,  254,  252.  251,  251 

E,  101,  102,  105,  105 

E,  106,  105,  102,  102 

E,  109,  106.  254,  254 

E,  105,  106,  109.  109 

E.  251.  253,  254,  254 

E.  253,  251,  178,  178 

E,  256.  254,  253,  253 

E,  254,  256,  109,  109 


E.  89,  93,  245,  245 

E,  245,  248,  249,  249 

E,  248,  245,  93.  93 

E,  99,  98,  94.  94 

E.  93,  94.  98,  98 


-THRU- 
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E.  98,  97,  93,  93 
E,  93,  97,  248,  248 

TYPE ,  1  $MAT  ,  106  $REAL ,  2 

E,  1,  2,  306,  306 

E,  2,  3,  307,  307 


E, 

2, 

307, 

306, 

306 

E, 

3, 

4,  307,  307 

E, 

4, 

308, 

307, 

307 

E, 

4. 

101, 

308, 

308 

E, 

101, 

309, 

308, 

308 

£, 

101, 

102, 

309, 

309 

E, 

102, 

310, 

309, 

309 

E, 

102, 

103, 

310, 

310 

E, 

1, 

306, 

176, 

176 

E, 

306, 

307, 

176, 

176 

E, 

307, 

177, 

176, 

176 

E, 

307, 

308, 

178, 

178 

E, 

307, 

178, 

177, 

177 

E, 

308, 

309, 

251, 

251 

E, 

308, 

251, 

178, 

178 

E, 

309, 

310, 

252, 

2.52 

E, 

309, 

252, 

251, 

251 

E. 

310, 

103, 

252, 

252 

E, 

97, 

98, 

301, 

301 

E, 

98, 

302, 

301, 

301 

E, 

98, 

99, 

302, 

302 

E, 

99. 

100, 

302, 

302 

E, 

100, 

303, 

302, 

302 

E, 

100, 

173, 

303, 

303 

E, 

173, 

304, 

303, 

303 

E, 

173, 

174, 

304, 

304 

17 

^  f 

*  f~^t 

anc 

on.i 

mJ' 

or%A 

E, 

174, 

175, 

305, 

305 

E, 

97, 

301, 

248, 

248 

E, 

301, 

302, 

248, 

248 

E, 

248, 

302, 

249, 

249 

E, 

302, 

250, 

249, 

249 

E. 

302, 

303, 

250, 

250 

E, 

303, 

299, 

250, 

250 

E. 

303, 

304, 

299, 

299 

E, 

304, 

305, 

300, 

300 

E, 

304, 

300, 

299, 

299 

E, 

305, 

175, 

300, 

300 

C**'>  RE.AL  CONSTANT  DEFINITION 

R,  1,. 015,50,0, .015, 50, 30 
RMORE,.015,50, 15, .015,50,90 
RMORE,.015,50,90,.015,50,15 
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RMORE..015^0,30,.015^0.0 
R,  2, .3, 106, 90, .3, 106,0 
RMORE,.3,106,90,.3,106,0 
RMORE,.3,106,90,.3, 106,0 
RMORE,.3,106,90,.3,106,0 

C***  MATERIAL  DEFINITION 

ALPX,  30, 0.2500000E-05 
EX  ,  30,  670000.0 
NUXY,  30,0.3600000 
GXY ,  30,  246323.5 
DENS,  50,  1.000000 

ALPX,  50,0.5157655E-06$ALPY,  50,-0.22 15747E-06$ALPZ.  50,-0.2215593E-06 
EX  ,  50.0.1920000E-K)8$EY  ,  50,  1560000.  $EZ  ,  50,  1560000. 

NUXY,  50,0.1950000E-01$NUYZ,  50.0.2100000  $NUXZ,  50, 0.801000E-02 
GXY,  50.  820000.0 
DENS,  106,  1.000000 

ALPX.  106.0.4391691E-06.«ALPY,  106, 0.4441 189E-06$ALPZ,  l06,-0.4642736E-06 
EX  .  106,  91963000.  $EY  ,  106,  61118000.  $EZ  ,  IfXS,  61118000. 

NUXY,  106,01.051700E-01$NUYZ.  106, 2.400697E-01$NUXZ,  106, 01.051700E-01 
GXY,  106,  8381000. 

C***  LOADS  AND  CONSTRAINTS  DEFINITION 
KTEMP,  0 

DDEL,ALL  $FDEL,ALL  $NTDEL,ALL 

TEDEL,  1,99999  $TDEL,1,99999  $EPDEL, 1.99999. 1,1  $RP6„„1 

F  ,  303MZ  .  1000.00 

D  ,  1,ALL  ,  0.0000000 

D  ,  2,ALL  .  0.0000000 

D  ,  3,ALL  ,  0.0000000 

D  ,  4,ALL ,  0.000000 

T*\  im  ATT  n  rwvwT 
f  x\f  f  V«\/W/\/v/v 

D  ,  102.ALL ,  0.000000 
D  ,  103.ALL  ,  0.000000 
D  ,  176,ALL  .  O.OOLOOO 
D  ,  177,ALL  ,  0.000000 
D  ,  178,ALL  ,  0.000000 
D  ,  251,ALL,  0.000000 
D  ,  252,ALL ,  0.000000 
LWRITE 
/GOPR 


Appendix  E 


Sample  ANSYS  Output  File 


GEOMETRY  STORED  FOR  310  NODES  616  ELEMENTS 
TnLE=  BLADE5MX 

*****  POSTl  NODAL  DISPLACEMENT  LISTING 

LOAD  STEP  1  ITERATION^  1  S£CTION=  I 
TIME=  0.00000  LOADCASE*  1 

T.IE  FOLLOWING  X.Y2  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

1  0.2283E-06-0.1450E-06  0.2482  -0.8246E-OI-0.2335E-06  0.1792E-06 

2  -0.i302E-06  0.8l99E-0i  0.1849  -0.8043E-01-0,2723E-06  0.6654E-07 

3  ■0.2755E-06  0.1397  0.1002  -0.8131E-01'0.1494E-06-0.1079E-06 

4  ^.6522E-07  0.1654  -0.3092E-07-0.6561E-01-0.il  18E-07-0.2450E-06 

5  0.7092E-06  0.254  lE-04  0.2487  -0.8299E-01-0.2331E-03  0.7063E-04 

6  -0.6694E-04  0.8227E-01  C.1832  -0.8298E-0l-0.7047E.03  0.1428E-02 
7 -0.9195E-04  0.1400  0.1005  -0.8255E-01-0.3136E-03  0.il29E-02 

8 .0.8979E-04  0.1658  0.3806E-04-0.8374E-01  0.3025E-02  0.1229E-01 

9  -0.71 12E-06  0.1869E-04  0.2495  -0.831 1E-01-0.2254E-03-0.1435E-04 

10  -O.I  lOOE-03  0.8234E-01  0.1859  -0.8304E-01-0.3756E-03  0.1526E-03 
11 -0.1601E-03  0.1402  0.1009  -0.8306E-01-0.4591E-03  0.4760E-03 

12  -0.1681E-03  0.1661  0.2943E-03-0.8301E-01-0.4955E-03  0.1485E-O2 

*****  POSTl  NODAi  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 

THE  FOLLOWING  X.Y^  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

13  0.6165E-06  0.3206E-07  0.2503  -0,8328E-0l-0.278lE-03-0.2460E4)5 

14  -0.1491E-03  0.8250E-01  0.1866  -0.8324E-01-0.2747E-03  0.9224E-04 

15 -0.23 19E-03  0.1405  0.1014  .0.8320E-01-0.2429E-03  0.2577E-03 

16  -0.25  lOE-03  0. 1664  0.6454E-03-O.83 17E-01-0.7682E-04  0.6402B-03 

17  0.2653E-05-0.9957E-05  0.2513  -0.8344E-01-0.3180E-03-0.1159E-05 

18  -0.1889E-03  0.8265E-01  0.1875  -0.8340E-01-0.2912E-03  0.5352E-04 

19  -0.301  lE-03  0.1408  0.1021  -0.8336E-01-0.2565E-03  0.1227E-03 

20  -0.3316E-03  0.1667  0.1 171E-02-0.8332E-OI-0.1905E-03  0.2422E-03 

21  0.4029E-05-0.1 126E-04  0.2525  -0.8359E-01-0.3570E-03  0.14<>9E-05 

22  -0.2279E-03  0.8280E-01  0.1885  -0.8356E-01-0.3235E-03  0.4998E-04 

23  -0.3685E-03  0,141 1  0.1030  -0.8352E-01-0.2786E-03  0.9272E-04 

24  -0.4niE-03  0.1670  0.1821E-02-0.8349E-OI-0.2149E-03  0.I447E-03 


•**■**  POSTl  NODAL  DISPLACEMENT  LISTING  *'*'••* 

LOADS'IEP  1  ITERATION*  I  SECTION*  1 
TIME*  O.OOOX)  LOADCASE*  1 

THE  FOLLOWING  X,YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

25  0.4624E-05-0.7447E-05  0^1537  -0.8373E-01-0.3962E-03  0^56E-05 

26  -0.2664E-03  0.8295E-01  0.1896  -0.837 1E-01-0.3598E-03  0.4692E-04 

27 -0.4349E-03  0.1413  0.1039  -0.8368E-01-0.3133E-03  0.8112E-04 

28  -0.4901E-03  0.1673  0.2605E-02-0.8365E-01-0.2551E-03  0.1086E-03 

29  0.4668E-05-0.1461E-05  0.2551  -0.8388E-01-e.4353E-03  0.2056E-05 

30  -0.3048E-03  0.831  lE-01  0.1909  -0.8387E-01-d.3975E-03  0.4624E-04 

31  -0.5009E-03  0.1416  0.1051  -0.8384E-01-0.3501E-03  0.7815E-04 

32  -0.5689E-03  0.1677  0.3519E-02-0.8382E-01-0.2934E-03  0.9716E-04 

33  0.4400E-05  0.4874E-05  0.2566  -0.8404E-01-0.4744E-03  0.1770E-05 

34  -0.3435E-03  0.8327E-01  0.1923  -0.8402E-01-0.4360E-03  0.4609E-04 

35  -0.5669E-03  0.1419  0.1063  -0.8400E-01-0.3882E-03  0.7698E-04 

36  -0.6475E-  03  0. 1680  0.4564E-02-0.8398E-01-0.3325E-03  0.9279E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION*  I  SECTION-  1 
TIME*  0.00000  LOADCASE*  1 

THE  FOLLOWING  X.Y^:  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

37  0.3976E-05  0.1064E-04  0.2583  -0.84 19E-01 -0.5 135E-03  0.1403E-05 

38  -0.3824E-03  0.8343E-01  0.1938  -0.8417E-01-0.4750E-03  0.4610E-04 

39 -0.6330E-03  0.1421  0.1077  -0.8416E-01-0.4269E-03  0.7657E-04 

40  -0.7260E-03  0.1683  0.5740E-02-0.8414E-01-0.3716E-03  0.91 14E-04 

41  0.3478E-05  0.1542E-04  02600  -0.8435E-01-0.5528E-03  0.1028E4)5 

42  -0.4215E-03  0.8359E-01  0.1955  -0.8433E-01-0.5142E-03  0.4609E-04 
43 -0.6993E-03  0.1424  0.1092  -0.8431E-01-0.4659E-03  0.7641E-04 
44  -0.8046E-03  0. 1686  0.7045E-02-0.8430E-01-0.4 108E-03  0.9045E-04 

tj  w.iyum-vA  -v/.«wjuc-vi-u.oyiii3-0j  U.OOZ/t-UO 

46  -0.4607E-03  0.8374E-01  0.1973  -0.8449E-01-0.5535E-03  0.4603E-04 

47  ■0.7658E-03  0.1427  0.1108  -0.8447E-01-0.5051E-03  0.7633E-04 

48  -0.8832E-03  0.1689  0.8481E-02-0.8445E-01-0.4499E-03  0.9017E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION*  1  SECTION*  1 
TIME*  0.00000  LOADCASE*  1 

THE  FOLLOWING  X,Y2  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

49  0.2374E-05  0.2I32E-04  0.2640  -0.8466E-01-0.63I6E-03  0.3082E-06 

50  -0.5001E-03  0.8390E-OI  0.1992  .0.8464E-01-0.5929E-03  0.4592E-04 
51 -0.8323E-03  0.1429  0.1125  -0.8463E-01-0.5443E-03  0.7631E-04 

52  -0.9618E-03  0. 1692  0.1005E-01-0.8461E-OI-0.4890E-03  0.9007E-04 

53  0. 1777E-05  0.2224E-04  0.2662  -0.8482E-01  -0.671 1E-03-0.4259E-07 

54  -0.5395E-03  0.8406E-01  0.2012  -0.8480E-01-0.6324E-03  0.4579E-04 

55  -0.8990E-03  0. 1432  0. 1 144  -0.8479E-01  -0.5837E-03  0.7634E-04 
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56  -0.1040E-02  0.1696  0.1 174E-01-0.8477E-01-0.5280E-03  0.9007E-04 

57  0.1 153E-05  0.2168E-04  02685  -0.8498E-01-0.7106E-03-0.3972E-06 

58  -0.5790E-03  0.8422E-01  0.2034  -0.8496E-01-0.6721E-03  0.4571E-04 

59 -0.9656E-03  0.1435  0.1164  -0.8494E-01-0.6231E-03  0.7649E-04 

60  ~0.1119E-02  0.1699  0.1357E-01-0.8493E-01-0.5670E-03  0.9021E-04 

f  •***  POSTl  NODAL  DISPLACEMENT  USTING  ***** 


LOAD  STEP  1  rreRATION=  1  SECTI0N=  1 
TIME=  0.00000  LOADCASE*  1 

THE  FOLLOWING  X.Y^i  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  U2  ROTX  ROTY  ROTZ 

61  0.5119E-06  0.1953E-04  0.2709  -0.8514E-01-0.7504E-03-0.8223E-06 

62  -0.6186E-03  0.8437E-01  0.2057  -0.85 12E-01-0.7120E-03  0.4583E-04 

63 -0.1032E-02  0.1437  0.1186  -0.8510E-01-0.6628E-03  0.7702E-04 

64  -0.1198E-02  0.1702  0.1552E-01-0.8508E-01-0.6060E-03  0.9064E-04 

65  -0.1224E-06  0.1570E-04  0.2735  -0.8530E-01-0.7898E-03-0.8985E-06 

66  -0.6581E-03  0.8453E-01  0.2082  -0.8528E-01-0.7525E-03  0.4681E-04 
67 -0.1099E-02  0.1440  0.1209  -0.8526E-01  •0.7031E-03  0.7868E-04 

68  -0.1276E-02  0.1705  0.1761E-01-0.8525E-01-0.6448E-03  0.9202E-04 

69  -0.6993E-06  0. 1005E-04  0.2762  -0.8545E-01-0.8353E-03-0.53 14E-05 

70  -0.6975E-03  0.8468E-01  0.2107  -0.8544E-01-0.7974E-03  0.4924E-04 

71 -0,1 166E-02  0.1443  0.1233  -0.8543E-01-0.7456E-03  0.8369E-O4 

72  -0.1355E-02  0.1708  0.1982E-01-0.8541E-01-0.6838E-03  0.96i9E-C4 

POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION*  1  SECHON*  1 
TIME*  0.00000  LOADCASE*  I 


THE  FOLLOWING  X.YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 


NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

73  -0. 1 168E-05  0.2545E-05  0.2790  -0.8561E-01-0.8485E-03  0. 1495E-04 

74  -0.7365E-03  0.8483E-01  0.2134  -0.8560E-01-0.8374E-03  0.6572E-04 


A  0€et\n  At  A  ^At  A  t  AAt 


76  -0.1434E-02  0.171 1  0.2216E-01-0.8557E-01-0.7233E-03  0.1095E-03 

77  -0.1 141E-05-0.7248E-05  0.2819  -0.8577E-01-0.9971E-03-0.6374E-04 

78  -0.  /748E-03  0.8498E-01  0.2163  -0.8576E-01-0.9716E-03  0.7177E-04 
79 -0.1298E-02  0.1448  0.1285  -0.8575E-01-0.8764E-03  0.1399E-03 

80  -0.1512E-02  0.1715  0.2464E-01-0.8573E-01-0.7690E-03  0.1493E-03 

81  -0.1250E-07-0. 1645E-04  0.2850  -0.8594E-01-0.8935E-03  0.21 13E-04 

82  -0.8125E-03  0.8514E-01  C.2192  -0.8593E  01-0.8052E-03  0.1593E-03 

83 -0.1363E-02  0.1451  0.1313  -0.859 1  E-01-0.83 18E-03  0.2522E-03 

84  -0.1591E-02  0.1718  0.2724E-01-0.8589E-01-0.7812E-03  0.2726E-03 


*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 


LOAD  STEP  1  ITERATION*  1  SECTION*  1 
TIME*  0.00000  LOADCASE*  1 


THE  FOLLOWING  X.Y^  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 
85  0.3426E-05-0.1887E-04  0.2882  -0.8610E-0I-0.8927E-03  0.1602E-03 


86  -0.8497E-03  0.8529E-01  0.2223  -0.8609E-01-0.1473E-02  0.3925E-03 

87  -0.1428E-02  0.1453  0.1342  -0.8609E-01-0.1834E-02  0.9139E-03 

88  4).1663E-02  0. 1721  0J2998E-01-0.8601E-01-0.1229E-02  0.8389E-03 

89  0.90J9E4)5-0.1026E-05  0.2916  -0.8623E-01-0.9567E-03-0.6209E-04 

90  -0.8862E-03  0.8544E-01  0.2256  -0.8629E-01-0. 1727E-02  0.2659E-03 
91 -0.1491E-02  0.1456  0.1373  -0.8616E-01-0.9865E-03  0.7199E-03 

92  -0.1744E4)2  0.1724  0.3280E-01-0.8626E-01  0.5281E-03  0.1162E-03 

93  0.8251E4)5  0.3780E-04  0.2951  -0.8635E-01-0.1 197E-02  0.4802E-03 

94  -0.9152E-03  0.85520-01  0.2290  -0.8624E-01-0.2103E-02  0.6081E-03 

95 -0,1555E-02  0.1459  0.1404  -0.8695E-01 0.1485E-02-0.1051E-02 

96  -0.1833E-02  0.1728  0.3575E-01-0,8549E-01-0.4080E-02  0.9093E-02 

*****  posTl  NODAL  DISPLACEMENT  USTiNG  ***** 


LOAD  STEP  1  ITERATIONS  1  SECTION=  1 
TIME-  0.00000  LOADCASE*  1 

THE  FOLLOWING  X.Y^  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

97  0.1600E-05  0.4343E-04  0.2986  -0.8676E-01-0.9476E-03  0.5884E-06 

98  -0.9377E-03  0.8575E-01  0.2324  4).8880E-01-0.9475E-03  0.4544E-06 

99 - 0.1599E-02  0.1461  0.1439  -0.8790E-01-0.9473E-03-0.2643E-07 

100 - 0.1893E-02  0.1730  0.3913E-01-0.1040  -0.9469E-03-0.9950E-07 
101  0.4105E-06  0.1363  -0.1725  -0.5966E-01  0.2545E-06-0.1400E-06 
102 .0.6052E-07  0.7385E-01-0.3375  -0.7488E-01 0.4199E-06  0.31 17E-06 

103  -0. 1633E-05  0.1 145E-05-0.4963  -0.7847E-01 0.2525E-06  0.8813E-06 

104  4).5301E-04  0.1365  -0.1727  -0.8408E-01  0.6142E-02  0.7814E-02 

105  -0.1845E-O4  0.7431E-01-0.3371  -0.8262E-01  C.2739E-02  0.2078E-02 

106  -0.1708E-04  0.808  lE-04-0.4967  -0.8438E-01 0.1691E-02-0.3974E-03 
107 -0.1082E-03  0.1368  -0.1729  -0.8296E-01  0.3590E-04  0.3608E-02 
108  -0.4432E-04  0.7400E-01-0.3386  -0.8302E-01-0.4169E-04  0.2629E-02 


*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 


LOADS'IEP  1  ITERATIONS  1  SECTION=  I 
TIMES  0.00000  LOADCASES  1 

THE  FOLLOWING  X,Y.Z  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 
109  0. 121  lE-05  0.693  lE-04.0.4975  -0.8247E-01-0.6907E-03  0.6106E-03 
110-0.1729E-03C.1371  -0.1727  -0.8311E-01  0.5550E-03  0.1316E-02 

111  -0.7552E-04  0.7434E-01-0.3384  -0.8311E-01 0.9901E-03  0.1157E-02 

112  0.4388E-05  0.4588E-04-0.4981  -0.8328E-01  0.6848E-03-0.4932E-03 

113  -0.2345E-03  0.1373  -0.1725  -0.8329E-01  0,2957E-04  0.5057E-03 

114  -0.1083E-03  0.7442E-01-0.3386  -0.8328E-()1  0.1907E-03  0.5795E-03 

115  0.2945E-05  0.3405E-04 -0.4985  -0.8328E-01-0.5128E-04  0.2282E-03 
116 -0.2987E-03  0.1376  -0.1722  -0.8345E-01-0.6216E-04  0.2193E-03 

117  -0.1429E-03  0.7456E-01-0.3386  -0.8344E-01  0.8841E-04  0.2242E-03 

118  0.1851E-05  0.2615E-04-0.4988  -0.8344E-01  0.1137E-03-0.i347E-04 
119 -0.3638E-03  0.1378  -0.1718  -0.8362E-01-0.1374E-03  0.1229E-03 
120  -0.1784E-03  0,7469E-01  0.3385  -0.8360E-01-0.2461E-04  0.1071E-03 

*****  posTl  NODAL  DISPLACEMENT  LISTING  ***** 


LOAD  STEP  1  riERATIONs  1  SECTION*  1 


TIME-  0.00000  LOADCASE-  1 


THE  FOLLOWING  X.Y^  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

121  0.6841E-06  0J2160E-04-0.4990  -0.8359E-01  0.2737E-04  0.1387E-04 

122  -0.4292E-03  0.1381  -0.1712  -0.8378E-01-0.1878E-03  0.9022E-04 

123  -0.2142E-03  0.7482E-01-0.3383  -0.8376E-01-0.8478E-04  0.6258E-04 

124  -0.1328E-06  0.1914E-04-0.4991  -0.8373E-01-0.5041E-05  0.3202E-05 

125  -0.4946E-03  0.1383  -0.1705  -0.8394E-01-0.2321E-03  0.7927E-04 

126  -0.2500E-03  0.7496E-01-0.3379  -0.8391E-01-0.1341E-03  0.4766E-04 

127  -0.7424E-06  0.1797E-04-0.4990  -0.8389E-01-0.4690E-04  0. 1786E-05 

128  -0.5599E-03  0.1386  -0.1697  -0.8410E-01-0.2735E-03  0.7555E-04 

129  -0.2857E-03  0.7510E-01-0.3374  -0.8407E-01-0.1775E-03  0.4253E-04 

130  -0.1231E-05  0.1769E-04-0.4988  -0.8404E-01-0.8667E-04  0.7756E-06 

131  -0.6250E-03  0.1389  -0.1687  -0.8426E-01-0.3137E-03  0.7419E-04 

132  -0.3213E-03  0.7524E-01-0.3367  -0,8423E-01-0.2185E-03  0.4078E-04 


*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 


LOAD  STEP  1  ITERATION*  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 

THE  FOLLOWING  X.Y 2  DISPLACEMENTS  APJE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

133  -0.1659E-05  0.1812E-04-0.4984  -0.8420E-01-0.1.264E-03  0.3820E-06 

134  -0.6900E.03  0.1391  -0.1676  -0.8442E-01-0.3532E-03  0.7363E.04 

135  ■0.3568E-03  0.7538E.01-0.33S9  -0.8439E-01-0.2584E-03  0.4014E-04 

136  -0.2055E-05  0.1918E-04-0.4979  -0.8436E4)l-0.1660E-03  02365E-06 

137  -0.7549E-03  0.1394  -0.1663  -0.8458E-01-0.3923E-03  0.7336E-O4 

138  -0.392  lE-03  0.7553E-01-0.3350  -0.8455E-0i-0.2978E-03  0.3989E-04 

139  -0.2435E-05  0.2080E-044).4973  -0.8452E-01-C.2054E-03  0.2574E-06 

140  -0.8196E-03  0.1397  -0.1650  -0.8474E-01-0.4312E-03  0.7326E-04 

141  -0.4274E-03  0.7567E-01-0.3339  -0.8470E-01-0.3369E-03  0.3982E-04 

142  -0.2810E-05  0.2298E-044).4963  -0.8467E-01-0.2448E-03  0.3855E-06 

143  -0.8843E-03  0.1399  -0.1635  -0.8489E-01-0.4698E-03  0.7335E-04 

144  -0.4626E-03  0.7581E-01-0.3328  -0.8486E-01-0.3758E-03  0.3987E-04 

*****  POSTl  NODAL  DISPLACEMENT  USTING  ***** 

LOAD  STEP  1  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 

THE  FOLLOWING  X.YJZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

145  -0.3190E-05  0.2570E-04-0.4957  -0.8483E-01-0.2841E-03  0.5941EO6 

146  -0.9490E-03  0.1402  -0.1619  -0.8505E-01-0.5082E-03  0.7379E-04 

147  -0.4978E-03  0.7595E-01-0.3314  -0.8502E-01-0.4144E-03  0.4012E-04 

148  -0.3577E-05  0.2898E-04-0.4947  -0.8499E-01  0.3234E-03  0.8722E-06 

149  -0.1014E-02  0.1404  -0.1601  -0.852 lE-Ol -0.546 lE-03  0.7522E-04 

150  -0.5330E-03  0.76  lOE-0 1-0.3300  -0.8517E-01-0.4528E-03  0.4086E-04 

151  -0.3965E-05  0.3286E-04-0.4935  -0.8514F.-0I-O.3628E-03  0.1226E-05 
152 -0.1078E-02  0.1407  -0.1582  -0.8537E-01-0.5833E-03  0.7951E-04 

153  ■0.5683E-03  0.7624E-01-0.3284  -0.8533E-01-0.4906E-03  0.4298E-04 

154  4).4324E-05  0.3733E-04-0.4922  -0,8530E-01-0.4021E-03  0.1685E-05 


155  ^.1143E-02  0.1410  -0.1562  -0.8552E-01-0.6186E-03  0.9233E-04 

156  -0.6036E-03  0.7638E.01-0.3267  -0.8548E-01-0.5270E-03  0.491  lE-04 

*****  POSTl  NODAL  DISPLACHvIENT  USTING  ***** 

LOAD  STEP  1  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 


THE  FOLLOWING  X,Y2  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  DZ  ROTX  AOTY  ROTZ 
157  -0.4588E-05  0.423  lE-04-0.4908  -0.8545E-0i-0.4415E-03  0.2323E-05 
158 -0.1209E-02  0.1412  -0.1541  -0.8568E-01-0.6475E-03  0.1286E-03 

159  -0.63y2E-03  0.7652E-01-0.3249  -0.8563E-01-0.5590E-03  0.6612E-04 

160  -0.4638E-05  0.476 lE-04-0.4893  -0.8560E-01-0.4809E-03  0.3361E-05 
161 -0.1275E-02  0.1415  -0.1518  -0.8584E-01-0.660CE-03  C.23d2E-03 

162  -0.6754E-03  0.7666E-01-0.3229  -0.8579E-01-0.5799E-03  0.1 185E-03 

163  -0.4264E-05  0.5270E-04-0.4876  -0.8575E-01-0.5221E-03  0.3852E-05 

164  -0.1343E-02  0.1417  -0.1494  -0.8598E-01-0.6383E-C3  0.5498E-03 

165  -0.7129E-03  0.7681E-01-0.3208  -0.8596E-01-0.5539E-03  0.2543E-03 

166  -0.3085E-05  0.5658E-04-0.4858  -0.859 1E-01-0.5570H-03  0.1756E-04 

167  -0.1412E4)2  0.1420  -0.1468  -0.8627E-01  0.1229E-03  0.1 1 13E-02 

168  -0.7539E-03  0.7698E-01-0.3185  -0.8612E-01-0.5036E-03  0.7403E-03 

*****  POSTl  NODAL  DISPLACEMENT  USTING  ***** 

LOAD  STEP  1  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 


THE  FOLLOWING  X.YZ  DISPLACEMF.NTS  ARE  IN  NODAL  COORDINATES 


NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 
169  -0.4994E-06  0.564  lE-04-0.4839  -0.8607E-01-0.6280E-03-0.3022E-04 
170 -0.1483E-02  0.1424  -0.1438  -0.8449E-01-0.2202E-02  0.8340E-02 

171  -0.7818E-03  0.7748E-01-0.3153  -0.8615E-01-0.7363E-03  0.2577E-02 

172  0.1782E-05  0.3875E-04-0.4817  -0.8660E-01-0.4814E-03-0.1184E-03 
173 -0.i559E-u^  0.1425  -0.1412  -0.1097  -0.9461E-03  0.i274E-07 

174  -0.8446E-03  0.7725E4) 1-0.3 137  -0.9442E-01-0.9461E-03  0.3650E-06 

175  -0. 1 149E-06  0.438 lE-04-0.4797  ■0.9037E-01-0.9464E-03  0.2743E-06 

176  0.3240E-06-0.8199E-01  0.1849  -0.7589E-01-O.I797E.06  0.1829E-06 

177  0.1916E-06  0..);»7  0.1002  -0.7287E-01-e.5586E-07  0.2508E-06 

178  -0.1525E-06-0.l6.54  0.2653E-%-0.4791E-01  0.1620E-06  0.2933E-06 

j79  0.6537E-04-0  S227E-01  0.1852  -0.8329E-01-0.8238E-03-0.8916E-03 
180  0.9284E-04-0.i400  0.1004  -0.8350E-01-0.9032E-03-0.2492E  02 


*****  POST  1  NODAL  DISPLACEMENT  LISTING  ***** 


LOAD  STEP  I  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 

THE  PC.  OWING  X,YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROT'Y  ROTZ 

181  0.I0I8E-03-O.1659  -0.3114E-03-0.8462E-01  0.1070E-02-0.3890E-02 

182  0.1 133E-03-0.g239E-01  0.1858  -0.8314E-01-0.2151E-03-0.4507E-04 

183  0.169gE-03-0.l403  0.1007  -O.S300E-01-0.1638E-03-0.2147E-04 

184  U.1834E-03-0.1661  0.2506E-03-0.8280E-01-0.3647E-03-0.1335E-03 


185  0.1559E-03-0.8255E-01  0.1866  -0.8329E-01-0.2526E-03-0.5205E-04 

186  0.2425E-03-0.1406  0.1013  -0.8317E-01-0.2132E-03-0.1017E-03 

187  0.2593E-03-0.1664  0.7082E-03-0.8301E-01-0.1247E-03-0.1559E-03 

188  0.1965E-03-0.8269E-01  0.1875  -0.8342E-01-0.2907E-03-0.4630E-04 

189  0.3095E-03-0.1408  0.1021  -0.8332E-01-0.2480E-03-0.7179E-04 

190  0.3349E-03-0.1667  0.1248E-02-0.8324E-01-0.1838F.-03-0.8376E-04 

191  0.2354E-03-0.8282E-01  0.1885  -0.8355E-01-0.3271E-03-0.4274E-04 

192  0.3743E-03  0.14U  0.1030  -0.8349E-01-0.2795E-03-0.727.5E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING 

LOAD  STEP  1  ITERATION-  1  SECTION-  I 
TIME-  0.00000  LOADCASE-  I 

THE  FOLLOWING  X.YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

193  0.4120E-03-0.1670  0.1892E-02-0.8344E-01-0.2122E-03-0.8890E-04 

194  0.2735E-03-0.8296E-01  0.1896  -0.8370E-01-0.3628E-03-0.4235E-04 

195  0.4390E-03-0.1413  0.1040  -0.8366E-01-0.3138E-03-0.7360E-04 

196  0.4900E-03-0.1673  0.2660E-02-0.8363E-01-0.2504E-03-0.8878E-04 

197  0.31 15E-03-0.8310E-01  0.1909  -0.8385E-01-0.3995E-03-0.4284E-04 

198  0.5042E-03-0.1416  0.1051  -0.8383E-01-0.3502E-03-0.7496E-04 

199  0.5684E-03-0.1676  0.3559E-02-0.8381E-01-0.2894E-OJ-0.9007E-04 

200  0.3497E-03-0.8325E-01  0.1923  -0.8401E-01-0.4373E-03-0.4361E-04 

201  0.5698E-03-0.1418  0.1063  -0.8399E-01-0.3880E-03-0.7606E-04 

202  0.6472E-03-0.1680  0.4590E-02-0.8397E-OI-0.3291E-03-0.9094E-04 

203  0.3881E-03-0.8340E4)1  0.1938  -0.8417E-01-0.4758E-03-0.4433E-04 

204  0.6359E-03-0.1421  0.1077  -0.8415E-01-0.4267E-03-0.7684E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***•* 

LOAD  STEP  I  ITERATION-  1  SECTION-  i 
TIME-  0.00000  LOADCASE-  1 

THE  FOLLOWING  X,Y.Z  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

205  0.7262E-03-0.1683  0.5754E-02-0.8414E-01-0.3690E-03-0.9156E-04 

206  0.4267E-03-0.8355E-01  0.1955  -0.8433E-OI-0.5149E-03-0.4491E-04 

207  0.7023E-03-0.1424  0.1092  -0.8431E-01-0.4660E-03-0.7735E-04 

208  0.8053E-03-0.1686  0.705 1E-02-0.8430E-01-0.4090E-03-0.9197E-04 

209  0.4655E-03-0.8371E-01  0.1973  -0.8449E-01-0.5543E-03-0.4538E-04 

210  0.7689E-03-0.1426  0.1108  -0.8447E-01-0.5057E-03-0.7767E-04 

211  0.8845E-03-0.1689  0.8481E-02-0.8445E-01-0.4491E-03-0,9222E-04 

212  0.5044E-03-0.8386E-01  0.1992  -0.8464E-01-0.5939E-03-0.4576E-04 

213  0.8357E-03-0.1429  0.1125  -0.8463E-01-0.5455E-03-0,7786E-04 

214  0.9638E-03-0.1692  0.1005E-01-0.8461E-01-0.4892E-03-0.9234E-04 

215  0.5433E-03-0.8402E-01  0.2012  -0.8480E-01-0.6337E-03-0.4606E-04 

216  0.9025E-03-0.1431  0.1144  -0.8478E-01-0.5855E-03-0.7793E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION-  1  SECTION-  1 
TIME-  0.00000  LOADCASE-  1 

THE  FOLLOWING  X.YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 


NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

217  0.1043E-02-0.1695  0.1174E-01-0.8477E-01-0.5293E-03-0.9234E-04 

218  0.5822E-03-0.8417E-01  02034  -0.8496E-01-0.6735E-03-0.4628E-04 

219  0.9694E-93-0.1434  0.1164  -0.8494E-01-0.6257E-03-0.7790E-04 

220  0.U22E-02-0.1698  0.1357E-01-0.8493E-01-0.5695E-03-0.9223E-04 

221  0.6212E-03-0.8433E-01  02057  -0.851 1E-01-0.7135E-03-0.4641E-04 

222  0.1036E-02-0.1437  0.1186  -0.8510E-01-0.6659E-03-0.7776E-04 

223  0.1?X)2E-02-0.1701  0.15.54E-01-0.3508E-01-0.6096E-03-0.9199E»04 

224  0.6602E-03-0.8449E-01  0.2082  -0.8527E-01-0.7534E-03-0.4641E-04 

225  0.1103E-02-0.1439  0.1209  -0.8525E-01-0.7061E-03-0.7742E-04 

226  0.1281E-02-0.1705  0.1764E-01-0.8524E-01-0.6496E-03-0.9162E-04 

227  0.6992E-03-0.8464E-01  0.2108  -0.8542E-01-0.7944E-03-0.4662F-04 

228  0.1 170E-02-0. 1442  0.1233  -0.8541E-01-0.7468E-03-0.7718E-04 

*****  POSTl  NOD/iL  DISPLACEMENT  LISTING  *•*** 

LOAD  STEP  1  ITERATIONS  1  SECnON=  1 
TIMES  0.00000  LOADCASE®  1 

THE  FOLLOWING  X,Y7-  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

229  0.1361E-02-0.1708  0.1987E-01-0.8540E-01-0.6896E-03-0.9133E-04 

230  0.738m-03-O.8480E-0l  0.2135  -0.8557E-01-0.8305E-03-0.4507E-04 

231  0.123/11-02-0.1445  0.1259  -0.8556E-01-0.7847E-03-0.7546E-04 

232  0.1440E-02-0.1711  0.2224E-01-0.8556E-01-0.7284E-03-0.9046E-04 

233  0.7765E-03-0.8496E-01  0.2163  -0.8571E-01-0.8870E-03-0.5491E-04 

234  0.1303E-02-0.1447  0.1286  ■0.8572E-01-0.8360E-03-0.8137E-04 

235  0.1519E-02.0.1714  0.2473E-01-0.8573E-01-0.7731E-03-0.9701E-04 

236  0.8141E-03-0.8512E-01  0.2193  -0.8586E-01-0.8836E-03-0.5124E-04 

237  0.1368E-02-0.1450  0.1314  -0.8587E-01-0.8260E-03-0.6253E-04 

238  0.1598E-02-0.1717  0.2737E-01-0.8591E-01-0.7836E-03-0.6628E-04 

239  0.85 12E-03-0.8530E-01  0.2224  .0.8601E-01-0.1085E-02-0.1141E-03 

240  0.143  lE-02-0. 1453  0.1343  -0.8597E-01-0.1064E-02-0.1831E-03 

POSTl  NODAL  DISPLACEMENT  LISTING  *•••• 

LOAD  STEP  1  ITERATIONS  1  SECTlON=  1 
TIMES  0.00000  LOAD  CASE®  1 

THE  FOIXOWING  X.YZ  DISPLACEMENTS  AFE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

241  0.1677E-02-0.1720  0.3014E-01-0.8593E-01-0.9888E-03.0.3213R-03 

242  0.8882E-03-0.8548E-01  0.2255  -0.8630E-01-0.9774E-03-0.2885E  03 

243  0.1486E-02-0.1456  0.1373  -0.8614E-01-057.ME-03-0.2268E-03 

244  0.1757E-02-0.1723  0.3306E-01-0.8710E-01-0.3391E-03  0.9434E-OJ 

245  0.9283E-03-0.8566E-01  0.2289  -0.8598E-01-0.2814E-02-0.1298E-02 

246  0.1563E-02-0.1460  0.1403  -0.8540E-01-0.4198E-02-0.4295E-02 

247  0.1827E-02-0.1727  0.3641E-01-0.7997E-01-0.2706E-02-0.6971E-02 

248  0.9401E-03-0.8566E-01  0.2324  -0.9340E-01-0.9474E-03  0.6722E-06 

249  0.1600E-02-0. 1460  0.1439  -0.9656E-01-0.9466E-03  0.9592E-06 

250  0.1894E-02-0.1729  0.3913E-01-0.1226  -0.9467E-03  0.4926E-06 

251  -0.3625E-06-0. 1363  -0.1725  -0.4332E-01  0.2258E-06  0.1503E-06 

252  -0.7122E-06-0.7385E-01-0.3375  -0.6374E-01  0.3074E-06  0.3940E-06 


**«■**  POSTl  NODAL  DISPLACEMENT  LISTING  •*••• 


LOAD  STEP  1  rTERATION=  1  SECTION*  1 
TIME-  0.00000  LOADCASE*  I 


THE  FOLLOWING  X,YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINa:  5S 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

253  0.6955E-04-0.1370  -0.1714  -0.8863E-01 0.2908E-02-0.4036E-02 

254  0.9713E-05-0.7479E-01-0.3357  -0.8870E-01 0.2344E-02-0.1802E-02 

255  0.9650E-04-0.1370  -0.1722  -0.8215E-01-0.2142E-03-0.7699E-03 

256  0.3763E-04-0.7413E-01-0.3380  -0.8196E-01 0.3014E-03-0.8906E-03 

257  0.1603E-03-0.1372  -0.1723  -0.8305E-01 0.1546E-03-0.1641E-03 

258  0.6870E-04-0.7437E-01-0.3381  -0.835 lE-Ol  0.2934E-03-0.4028E-03 

259  02250E-03-0.1374  -0.1723  -0.8324E-01-0.4686E-04-0.7246E-04 

260  0.1034E-03-0.7444E-01-0.3384  -0.8335E-O1 0.8070E-04-0.5995E-04 

261  0.2919E-03-0.1376  -0.1721  -0.8344E-C1-0.8589E-04-0.6043E-04 

262  0.1389E-03-0.7455E-01-0.3385  -0  ?350E-0l  0.1094E-04-0.4276E-04 

263  0.3589E-03-0.1378  -0.1717  -0.8.  62E-01-0.1380E-03-0.6539E-04 

264  0.1752E-03-0.7467E-01-0.3385  -O.c  jtj3E-0l-0.434lE-04-0.3l66E-04 


POSTl  NODAL  DISPLACEMENT  LISTING  ***** 


LOAD  STEP  1  ITERATION*  1  SECTION*  1 
TIME*  0.00000  LOADCAS&r  1 

THE  FOLLOWING  X.Y^Z  DISPLACEMELrTS  APxE  IN  NODAL  COORDINATES 

NODE  UX  UY  IjZ  ROTX  ROTY  ROTZ 

265  0.4256E-03-O.I381  -0.1712  -0.8379E.OI-0.1837E-03-0.6974E-04 

266  0.2113E-0343.7479E-01-0.3382  -0.8377E-01-0.9073E-04-0.3551E-04 

267  0.4919E-03-0.1383  -0.1705  -0.8395E-01-0.2275E-03-0.7280E-04 

268  0J2473E-O3-O.7493E-01-O.3379  -0.8392E-OI-0.1344E-03-0.3786E-04 

269  0.5578E-03-0.1386  -0.1697  -0.8411E-01-0.2698E-03-0.7448E-04 

270  0.283 1E-03-0.7507E-01-0.3374  -0.8407E-01-0.1765E-03-0.3951E-04 

271  0.6235E-03-0.1388  -0.1687  -0.8427E-01-0.3110E-03-0.7530E-04 

272  0.3186E-03-0.7520E-01-0.3367  -0.8423E-01-0.2174E-03-0.4032E-04 

273  0.6890E-03-0.1391  -0.1676  -0.8442E-01-0.3515E-03-0.7559E-04 

274  0.3540E-03-0.7534E-01-0.3359  -0.8439E-01-0.2575E-03-0.4060E-04 

275  0.7543F.-O3-O.1394  -0.1663  -0.8458E-01-0.3916E-03-0.7557E-04 

276  0.3892E-O3-O.7548E-OI -0.3350  -0.8455E-01-0.2972E-03-0.4056E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATION*  1  SECTION*  1 
TIME*  0.00000  LOADCASE*  1 

THE  FOLLOWING  X,Y;z;  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  ROTY  ROTZ 

277  0.8 196E-03-0. 1396  -0.1650  -0.8474E-01-0.4314E-03-0.7539E-04 

278  0.42/UE-03-0.7562E-01-0.3339  -0.8470E-01-0.33<36E-03-0.4032E-04 

279  0.8847E-03  0.1399  -0.1635  -0.8490E-01-0.4709E-03-0.7510E-04 
,280  0.4595E-03-0.7576F.-01-0.3328  -0.8486E-01-0.3758E-03-0.3996E-04 

281  0.9499E-03-0.1401  -0.1618  -0.8506E-01-0.5103E-03-0.7477E-04 

282  0.4946E-03-0.7590E-01'0.3314  -0.8502E-01-0.4147E-03-0.3951E-04 

283  0.1015E-02-0.1404  -0.1601  -0.8522F,-01-0.5494E-03-0.7443E-04 

284  0.5297E-03-0.7603F  01-0.3300  -0.8518E-01-0.4536E-03-0.3902E-04 

285  0.1080E-02-0.1406  -0.1582  -0.8538E-01-0.5884E-03-0.7423E-04 


286  0.5647E-03-0.7617E-01-0.3284  -0.8534E-01-0.4924E-03-C.3857E-04 

287  0.1 145E-02-0. 1409  -0.1562  -0.8554E-01-0.6273E-03-0.7438E-04 

288  0.6000E-03-0.7630E-01-0.3267  -0.8550E-0 1-0.53 12E-03-0.3843E-04 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATIONS  1  SECTION=  1 
TIMES  0.00000  LOADCASE-s  1 

THE  FOLLOWING  X.YZ  DISPLACEMENTS  ARE  IN  NODAL  COORDINATES 

NODE  UX  UY  UZ  ROTX  R07Y  ROTZ 

289  0.1210E-02-0.1411  -0.1540  -0.8571E-01-0.6650E-03-0.7816E-04 

290  0.6354E-03-0.7644E-01-0.3249  -0.8565E-01-0.5700E-03-0.3984E-04 

291  0.1276E-02-0.1414  -0.1517  -0,8587E-01-0.7006E-03-<).8176E-04 

292  0.6713E-03-0.7657E-01-0.3229  -0.8581E-01-0.60SOE-03-0.4673E-04 

293  0.1341E-02-0.1417  -0.1493  -0.8600E-01-0.739.5E-03-0.1918E-03 

294  0.7079E-03-0.7670E-01-0.3208  -0.8598E-01-0.6324E-03-0.8570E-04 

295  0.1409E-02-0.1419  -0.1468  -0.8665E-01-0.5437E-03  0.2074E-04 

296  0.7463E-03-0.7686E-01-0.3185  -0.8612E-01-0.6250E-03-0.2036E-03 

297  0.1474E-02-0.1422  -0.1440  -0.8189E-01-0.1004E-02-0.3539E-02 

298  0.7778E-03-0.7702E-01-0.3161  -0.8481E-01-0.6506E-03-0.9445E-03 

299  0.1559E-02-0.1424  -0.1412  -0.1262  -0.9464E-03-0.1304E-06 

300  0.8449E-03-0.7716E-01-0.3137  -0.1056  -0.9461E-03-0.1230E-06 

*****  POSTl  NODAL  DISPLACEMENT  LISTING  ***** 

LOAD  STEP  1  ITERATiONs  1  SECTIONS  1 
TD/IEs  O.OCOOO  LOAD  CASES  1 


THE  FOLLOWING  X.Y.Z  DISPLACEMENTS  ARE  IN  NOD/  L  COORDINATES 


NODE  TJX  UY  UZ  ROTX  ROTY  ROTZ 

301  0.1 193E-05  C.4347E-04  0.2324  -0.9610E-01-0.9475E-03  0.5418E-06 

302  0.7784E-060.4367E-04  0.1439  0.1118  -0.9471E-03  0.4921E-06 

303  0.3327E-06  0.4374E-04  0.3913E- 11-0.1697  -0.9470E-03  0.355(E-06 

'2fil  n  lOOi^C.A^  A  i4170C  A/l  n  1 

305  0.1605E-060.4377E  04-0.3137  -0.1031  •0.9462EO3  0.2380E-06 

306  0.1 155E-06O.1835E-06  0.1849  -0.7325E-01-0.2278E-06  0.1 1  lOE-06 

307  0.3081E-070.8591E-07  0.1002  -0.5778E-01-0.1545B-06  0.54%E-07 

308  0.0000  O.OCOO  0.0000  0.0000  0.0000  0.0000 

309  0.2274E-('3  0.3168E-06-0.1725  -0.47930-01  0.2268E-06  0.7406E-07 

310  -0.4158E-06  0.7867E-06-0.3375  -0.661  lE-Ol  0.3620E-06  0.3978E-06 


MAXIMUM^ 

NODE  250  100  124  303  104  8 

VALUE  0,1894E-02  0.1730  -0.4991  -0.1697  0.6142E-02  0.1229E-01 


